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TASK 2 

TEST OF SPECTRAL/SPATIAL CLASSIFIER 


INTRODUCTION 

Contemporary classifiers for analysis of remotely sensed data compare 
spectral measurements from each feature of each pixel to class statistics, 
computing a likelihood discriminant function associated with each class, 
and categorizing the point according to the class with the largest discri- 
minant function value. Each point is classified on the basis of its spectral 
measurements alone. One premise of this technique is that the objects of 
interest are large in comparison to the size of a point. If this were not 
so, a large proportion of points would be composites of several classes, 
making statistical pattern classification unreliable since pre-specif ied 
categories would be inadequate to describe actual states of nature. From 
this premise it follows that objects are represented by arrays of points, 
and that a statistical dependence exists between consecutive points. 
Contemporary classifiers fail to exploit the statistical dependence between 
adjacent points when assigning classes. 

The ECHO processors benefit from spatial information by first aggregating 
into groups points whose spectral responses are not significantly different 
in a statistical sense, and then applying a maximum likelihood classification 
rule to these homogeneous groups. Homogeneous objects are identified in 
a three step process. First, cells are formed by systematically partitioning 
the data into N by N sized blocks of pixels. The statistics of each cell are 
then compared to a homogeneity criterion. Points which do not comprise 
homogeneous groups are classified on a point-by-point basis, just as con- 
temporary classifiers classify all points. Statistics of adjoining homogeneous 
cells are then compared to annexation thresholds. Adjoining cells which 
appear to belong to the same statistical population are combined into a 
single object. 

Two separate ECHO algorithms have been developed. The first, Supervised 
ECHO, makes use of pre-specif ied class statistics to identify homogeneous 
objects. The second, Nonsupervised ECHO, identifies homogeneous objects 
without the use of class statistics. Consequently, those objects identified 
by the Nonsupervised algorithm may be used to aid in the training process. 

Past Work at LARS 


Much of the background research on the ECHO concept was performed at 
LARS during 1975 and is documented in the Final Report for 1975 [Y] , R. L. 
Kettig's doctoral thesis[2j, a LARS Information Note [3] , and in symposium 
proceedings [4] . 


*ECHO stands for Extraction and Classification of Homogeneous Objects. 
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The work related to the development, testing, and documentation of the 
Supervised ECHO algorithms is documented in Volume I of the Final Technical 
Report on NASA Contract NAS9-14970, June 1, 1976 - May 31, 1977 [5] . 

Included in that report are a more detailed background on the past develop- 
ment work on ECHO at LARS, an appendix containing program listings and 
documentation for the Supervised ECHO algorithms, the results of systematic 
tests of the Supervised processor on MSS data for agricultural regions as 
observed by the Landsat satellites, aircraft scanners, and simulated Thematic 
mapper data, and an example product of an object map enabling the determi- 
nation of the utility of object maps to a LACIE Analyst Interpreter in the 
selection and labeling of training fields. 
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Objectives 

The objectives for the contract extension were to: 

* Complete the Fortran implementation and documentation of the Non- 
supervised ECHO (Extraction and Classification of Homogeneous 
Objects), making the algorithms available for use by JSC. 

* Classify, at different parameter settings, Landsat, simulated 
Thematic Mapper and aircraft data sets. Compare the resulting 
CPU time required, full field performance, field center pixel per- 
formance and classification variability of results with those 
achieved by the perpoint and Supervised ECHO processors. 

* Evaluate the objects identified by the Nonsupervised algorithms in 
terms of two types of errors: 

1. more than one field on the ground being identified by ECHO 
as a single object, and 

2. a single field on the ground being subdivided into distinct 
objects. 
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DESCRIPTION OF WORK 


Produce Documented Fortran Programs 


Fortran listings and program abstracts for the Nonsupervised ECHO 
processor are presented in Appendix A. This processor performs field 
extraction without the benefit of class statistics. Statistics are 
necessary for classification of the objects identified, however. The 
software is designed to function in a two phase mode. The first phase is 
the Nonsupervised Field Extraction Algorithm. This phase of the processor 
partitions the data to be classified into a set of points and homogeneous 
fields and calculates the channel means and covariance matrix for each 
homogeneous field which is identified. 

The Nonsupervised processor utilizes a homogeneity test which compares 
a ij/^ij t0 a threshold t^; where i refers to feature or channel and j refers 
to thecell. If the standard deviation divided by the mean for feature i 
exceeds the user-specified threshold t . , the cell is "singular" and elements 
of the cell will be classified on a point-by-point basis. Adjoining, homo- 
geneous cells are annexed to fields on the basis of a two-step test, 1) that 
the channel variances of the field are equivalent to the channel variances 
of the cell and 2) that the channel means of the field are equivalent to 
the channel means of the cell. Should either of these criteria not be met 
for any channel, annexation will not take place. 

The research software was designed to run in two phases. The first 
phase merely calculated the cell mean and covariance matrices and wrote them 
on tape; the second phase proceeded to perform field extraction followed by 
classification. This process has two disadvantages. First, although an 
object map could be produced, it had to be produced in phase two, the same 
phase that required a class statistics deck for input in order to perform 
the classification. Second, since only cell mean and covariance matrices 
were written on the intermediate tape, when singular cells were encountered 
in phase two, they had to be classified as small samples, a cell at a time, 
rather than on a point-by-point basis, because data values for invididual 
points were not available to the phase two classifier. 

The Nonsupervised ECHO software was restructured to move the field 
extraction algorithm into phase one of the processor and produce an inter- 
mediate tape containing : 

* class means and covariance matrices for each homogeneous object 
identified. 

* an object map containing the mean for channel i of object j in 
every pixel of object j and the original data values for those 
pixels belonging to singular cells. 

* a tag array identifying whether a pixel falls in a singular cell 
(and should be classified individually), or inaa homogeneous object 
and hence should receive a class assignment based on the sample 
classification of the object to which it belongs). 
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The second phase of the restructured Nonsupervised processor reads 
the intermediate tape and the class statistics deck, performs a maximum 
likelihood point-by-point classification on points falling in singular cells, 
and a maximum likelihood sample classification of the homogeneous objects 
identified. Figure 2a-l presents the general processing flow for phase 
one (field extraction) of the Nonsupervised processor. The general processing 
flow for phase two (classification) is presented in Figure 2a-2. Inputs 
to outputs from the Nonsupervised ECHO processor's two phases are described 
in Table 2a- 1. 

The formats of the disk and tape files which support the Nonsupervised 
ECHO processors are presented in Appendix B. User documentation is available 
in the forms of an ECHO User's Guide [6] and an Echo Case Study [7] . 
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Figure 2a-l 
GENERAL FLOW OF 
NONSUPERVISED ECHO 
PHASE 1 























2a-8 


Table 2a-l 
NONSUPERVISED ECHO 
PROCESSOR 

INPUT OUTPUT 

Phase One : 

Channel Selection Intermediate Tape containing 

Cell Width object map, object statistics 

Cell Homogeneity Thresholds and pixel tag array. 

<°i/^i ^ or c * ianne l *■) 

Annexation Thresholds 

(mean and variance) 

Multispectral Image Storage 
Tape 

Intermediate Results Tape 
and File Specification 

Phase Two : 

Class Statistics Results File 

Intermediate Tape from 
Phase One 

Specification of Results 


Tape and File 
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Test the Nonsupervised ECHO Algorithms 
1. Data Sets 

The second objective of the ECHO Extension work is to test the Non- 
supervised ECHO algorithms on MSS data for agricultural regions. Data 
sets are to include Landsat, aircraft, and simulated Thematic Mapper data. 
Nine Landsat, one aircraft and six simulated Thematic Mapper data sets 
(form resolutions from one site and two resolutions from a second site) 
were selected for analysis. The data sets are summarized in Table 2a 
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Table 2a~2 




Data Sets 




(LANDSAT) 



CHANNEL 

NUMBER 


WAVE BAND 
UM 


1,5,9 


.5-. 6 


2,6,10 


.6-. 7 


3,7,11 


.7-. 8 


4,8,12 


.8-1.1 


Area 

Channles Used 

Data 

Collected 

Data Set 

Graham County, Kansas 

9, 10, 11, 12 

5/26/74 

LACIE/SRS 

Grant County, Kansas 

5, 6, 7, 8 

5/19/74 

LACIE/SRS 

Haskell County, Kansas 

9, 10, 11, 12 

5/27/74 

LACIE/SRS 

Kearny County, Kansas 

9, 10, 11, 12 

5/27/74 

LACIE/SRS 

Huntington County, Indiana 

1, 2, 3, 4 

7/7/73 

CITARS 

Shelby County, Indiana 

1, 2, 3, 4 

9/7/73 

CHARS 

White County, Indiana 

1, 2, 3, 4 

8/21/73 

CITARS 

Fayette County, Illinois 

1, 2, 3, 4 

8/21/73 

CITARS 

Lee County, Illinois 

1, 2, 3, 4 

8/5/73 

CITARS 
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Table 2a-2 (Continued) 


Simulated 

Thematic Mapper 

I CHANNEL 

WAVE BAND 

NUMBER 

pM 

1 

.45-. 52 

2 

.52-. 60 

3 

.63-. 69 

4 

.74-. 80 

5 

.80-. 91 

6 

1.55-1.75 

7 

10.4-12.5 

8 

.74-. 91 


Area 

Channels 

Used 

Date 

Resolution 

Williams County, ND 

1, 

2, 

3 , 

6, 8 

8/15/75 

30m 

Williams County, ND 

1, 

2, 

3 , 

00 

vO 

8/15/75 

40m 

Finney County, Kansas 

2, 

3 , 

5, 

6, 7 

7/6/75 

30m 

Finney County, Kansas 

2, 

3 , 

5, 

6, 7 

7/6/75 

40m 

Finney County, Kansas 

2, 

3 , 

5, 

6, 7 

7/6/75 

50m 


Finney County, Kansas 


2, 3, 5, 6, 7 


7/6/75 


60m 
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Area 


Table 2a-2 (Continued) 

Aircraft 

Wavebands Used Data Collected 

(yM) 


.4-. 49, .59-. 64, .65-. 69, 


Finney County, KS 


82-. 88, 1.53-1.62, 10.1-11.0 


7/6/75 
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2. Training Procedures 

a. Landsat Data Sets 

The LACIE training sets were created using ground truth information 
provided by JSC for both test and training fields. Odd-numbered fields 
appearing in each subclass were used for training; even-numbered fields 
were used for test. For Graham and Grant counties, statistics were generated 
by using the STATISTICS processor. For Haskell and Kearny counties sta- 
tistics were secured by clustering the training fields of each class into 
four subclasses which were then pooled on the basis of output from the 
SEPARABILITY processor. It should be noted that both training and test 
fields for the LACIE/SRS data are large enough to inset field boundaries 
approximately two pixels inside estimated field boundaries. This inset 
allows for any image misregistration which may occur between any two dates 
on the multitemporal input runs. Proportion estimate for the 1974 LACIE/SRS 
segments were available in ground truth packets provided by JSC. 

The CITARS training sets were originally created by a supervised 
procedure using ground truth information provided by the Agricultural 
Stabilization and Conservation Services (ASCS) to select both the training 
and test fields|8|. Those training sets were used without change in the 
ECHO tests. 

Five counties in the CITARS experiment were used as test sites for the 
Nonsupervised processor evaluation. A data set free of clouds which occured 
late in the growing season was required for each of the five counties. Dates 
from July, August or early September were selected. The classification 
results for all of the CITARS experiments are catalogued on a series of LARS 
tapes. The desired data sets were located on the catalogued CITARS tapes 
and the statistics which had been used for the CITARS experiment were obtained 
by using the LARSYS PUNCHSTATISTICS processor. The pooling of classes was 
determined by running the PRINTRESULTS processor and requesting training 
field results. This combination of requests produced a table of available 
classes in the statistics deck and also the informational names under 
which they were classified. By using this list it was possible to reconstruct 
the combinations of spectral classes and pooled spectral classes which were 
needed to reproduce the original classification. 

Training and test fields for each CITARS classification were secured 
and appropriate control cards added to evaluate the ECHO classifications 
of CITARS data sets. The proportion estimates used for the CITARS evalua- 
tion are estimates of the proportions of the various classes for the entire 
county made by the SRS. 

b. Simulated Thematic Mapper Data Sets 

Training sets from the simulated Thematic Mapper tests performed at 
LARS in 1976 were used for ECHO analysis. These training sets were generated 
by selecting fields of known cover types and clustering each informational 
class separately to define subclasses. 
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Color infrared photographic mosaic prints were made from photographic 
data collected concurrently with the scanner data. Informational class 
information provided by ground observations was transferred to clear plastic 
overlays on the mosaic print. The analyst could then easily locate the 
corresponding fields in his cluster maps and assign the field coordinates 
to the informational classes. 

Statistics were calculated for each training area and compared using 
the SEPARABILITY processor. Similar classes were combined, where indicated 
and the data set was used to classify the flightlline. Training areas were 
not excluded from the test fields since the test fields had been pre-selected 
for the entire flightline. 

Two to four subclasses were found in each informational class. The 
Kansas flight was an . exception. Because of severe line-to-line changes 
in signal level in the original 6 meter scanner data, the analyst was 
forced to create more spectral classes to account for the within-class 
variations due to excessive noise. This was most apparent in the 30 and 
40 meter resolution data. The effect was reduced but not eliminated in 
the 50 and 60 meter resolution data. Alternate fields were used for 
training and test decks, respectively. The procedure was repeated for each 
of the four resolution sizes. As resolution size increases, the number 
of spectral classes decreases. 

The entire training set selection procedure was repeated for each 
resolution size so that any effects on training set selection which might 
be caused by data resolution would be included in the analysis results. 

An example is the increasing difficulty and eventual impossibility of 
selecting samples from small, or narrow, fields as the resolution size 
increases. 

c. Aircraft 

A 6-meter aircraft scanner data set used to generate the Finney County, 
Kansas simulated Thematic Mapper data set was used. The same training and 
test fields used for the Thematic Mapper were available in the six meter 
data set. Because of the very large number of data points in the full 
data set, only the first two miles (one third of the total flight line) 
were used for the Nonsupervised ECHO evaluation. 

The data set was not corrected for sun and scanner angle effects. To 
compensate for these angle effects, training fields were distributed across 
the width of the flight line. At least six fields in each informational 
class were used in the training set. Fields from informational classes were 
clustered together into five spectral classes. All spectral classes from 
all informational classes were combined into a statistics deck and appro- 
priate pooling was done based on SEPARABILITY results. The classification 
results were strongly influenced by angle effects. 
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3. Dependent Variables 

There are six variables which were monitored to evaluate the ECHO 
algorithm. 

' CPU time 

Field center pixel classification performance 

* Training field classification performance 

‘ Full field classification performance 
RMS proportion estimate error 
Classification variability 

These variables are related to reasons for adopting a new classification 
technique: cost, accuracy, and usability of results. The CPU time required 

to perform a classification is one way to measure the cost of classification. 
Field center pixel, full field, and training field performances and RMS 
proportion estimate error are all wasy to evaluate the accuracy of the 
classifier. Classification variability is a measure of "salt and pepper 
effect" in classification results. 

The CPU time required to execute each of the ECHO classifications has 
been recorded so that the effects of varying the cell homogeneity and 
annexation thresholds may be monitored. The CPU time required to perform 
the perpoint classifications have been adjusted to reflect the increased 
efficiency of the LARSYS perpoint classifier which is coded in assembly 
language. Thus, the CPU time recorded for a perpoint classification is 
what a FORTRAN classifier would have required to perform the classification. 

The indices of classification performance were applies in several 
ways. Classification accuracy (identification) was evaluated utilizing 
field center pixel, "full field" and test field sample performances for 
all data sets. Proportion estimation was carried out for the Landsat and 
Simulated Thematic Mapper data sets. 

The training performance is the overall classification accuracy (number 
of training pixels correctly classified divided by the total number of 
training pixels) of the pixels used to calculate the class statistics. 

Field center pixel performance is the overall classification accuracy of 
pixels inset at least one pixel from the field boundary. For the registered 
LACIE/SRS data the field center pixels are inset at least two pixels from 
the field boundary. Although this procedure insures that the pixels exa- 
mined are not mixture pixels, it has the unfortunate effect of eliminating 
smaller fields from consideration. The third measure of classification 
accuracy, "full field" performance, includes those pixels on the boundaries 
of the fields in the classification performance. The "full field" pixels 
were generated by expanding the field center pixel boundaries one pixel in 
all directions. 

The RMS error of informational class proportion estimates for each 
flightline was found by calculating the percent of the flightline classi- 
fied as a particular class and comparing it with the ground-collected 
estimage using equation (1). 



RMS Error = J 
v 


( 1 ) 


l 

i=l 


(C^C'p 

“n 


where, N = number of informational classes, 

Ci= percent classified as informational class i, and 

C'i= percent of class i estimated from ground-collected data. 

RMS error is calculated for the Landsat and Thematic Mapper data 
runs. The Agricultural Stabilization and Conservation Service (ASCS) 
provided the ground truth proportion estimates for the simulated Thematic 
Mapper data set. Proportion estimates for the 1974 LACIE/SRS segments 
were provided in ground truth packets received from JSC. The SRS county 
proportion estimates were used to calculate RMS proportion error for the 
CITARS data set. 

Average variability is a measure of the rate of change from one 
information class to another. It should reflect the degree to which ECHO 
reduces the "salt and pepper effect" which is sometimes present in per- 
point classifications. Variability is calculated by systematically 
selecting 50 lines of the classified area, counting the number of informa- 
tion class changes, and dividing by the number of opportunities for class 
changes . 

Variability - NCC/ (50*(NS-1)) (2) 

Where: 

NCC = the number of class changes over the 50 selected lines, and 

NS = the number of classified pixels/lines. 
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4 . Results 

This section outlines the results of the tests performed on the 
Nonsupervised ECHO classifications and the comparison of the Nonsuper- 
vised ECHO processor with the Supervised ECHO and the perpoint 
classifiers. Results are discussed separately for each scanner type, 
Landsat, simulated Thematic Mapper, and aircraft. 

For all data sets considered, the Nonsupervised ECHO processor was 
run with the following parameter settings: 

1. Cell width of 2. 

2. Cell homogeneity thresholds of 0.05, 0.10, and 0.25. 

3. Both the mean and the covariance annexation thresholds 
set at 0.001, 0.01, and 0.1. 

a. Landsat Results 

Training and test information for the nine Landsat test data sets 
were drawn from the 1974 Kansas LACIE/SRS data sets and the 1973 CITARS 
data sets. Four LACIE/SRS and five CITARS data sets were considered. 

i LACIE/SRS Results 

Figures 2a-3 through 2a-8 present the average results for the four 
LACIE/SRS sites examined at a cell width of two for three homogeneity 
and three annexation parameters. Figure 2a-3 plots the average CPU time 
in seconds required by the perpoint classifier (represented by the line 
of 'P's) versus the average CPU time required by the Nonsupervised ECHO 
routine to classify the four LACIE/SRS data sets at each of nine Non- 
supervised ECHO parameter settings for 2 by 2 pixel cells. The cell 
homogeneity threshold is plotted along the horizontal axis. As this 
threshold increases, it becomes more likely a cell will be classified as 
a unit, less likely that a cell will be split and its constituent pixels 
classified individually. The dependent variable, CPU time, appears on 
the vertical axis. The cell-to-f ield annexation parameter for each cell 
homogeneity threshold is represented on the plot. A '1' appears in the 
position for the results achieved when both the mean and covariance 
annexation thresholds are set to lO - -^; a '2' when they are set to 10“^; 
and a '3' when they are set to 10 - ^. As the annexation thresholds 
become smaller, the likelihood that adjoining homogeneous cells will be 
annexed into a single field increases. When two or more annexation 
thresholds achieve the same performance, a star appears in the position 
on the plot. 

The statistical significance of the effects of the Nonsupervised 
ECHO parameters on the LACIE/SRS results are presented in Table 2a-3. 
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Figure 2a-3 


Nonsupervised ECHO CPU Requirements 
for the LACIE/SRS Data Sets 
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Figure 2a-4 


Nonsupervised ECHO Field Center Pixel Performance 
for the LACIE/SRS Data Sets 
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Figure 2a-5 

Nonsupervised ECHO Full Field Performance 
for the LACIE/SRS Data Sets 
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Figure 2a-6 

Nonsupervised ECHO Training Field Performance 
for the LACIE/SRS Data Sets 
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Figure 2a-7 

Nonsupervised ECHO RMS Proportion Estimate Error 
for the LACIE/SRS Data Sets 
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Figure 2a-8 

Nonsupervised ECHO Classification Variability 
for the LACIE/SRS Data Sets 
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Table 2a-3 

Overall Landsat Comparisons 
Effects of Parameters on LACIE Data 


Variable 

Homogeneity 

Thresholds 

Annexation 

Thresholds 

Homogeneity X 
Annexation 

Degrees of Freedom 

2,6 

2,6 

4,12 

CPU Time 

1.86 

6.69* 

21. 37*** 

Field Center Pixel 
Performance 

1.27 

7.02* 

1.90 

Full Field 
Performance 

1.31 

11.09** 

3.74* 

Training Field 
Performance 

.70 

3.60 + 

5.89** 

RMS Proportion 
Error 

. 86 

4 . 0 3 + 

1.18 

Classification 

Variability 

13.64** 

14.39** 

12.59*** 


Significance 

Levels 

Table entries are F-values 

+ 

10% 


* 

5% 


* * 

1% 


*** 

.1% 
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It can be seen from Figure 2a-3 that, as the cell homogeneity 
parameter increases from .05 to .25, the CPU time required is effected 
by the interaction between the homogeneity and the annexation threshold. 
This interaction is probably the result of two separate influences. 

The LACIE/SRS, CITARS, and aircraft data runs of the Nonsupervised 
processor were produced utilizing the "MAP" option. When this option 
is specified, the data value for each point which falls within a 
field is replaced by the mean value for the field. This replacement 
requires a significant amount of computer time. Therefore, as more 
fields are identified, more CPU time will be used. As the homogeneity 
parameter increases, there is a tendency for more cells to be identified 
as homogeneous, and therefore, it is possible for more fields to be 
identified. The effect of this is reflected by the increase in CPU 
time of the runs with a 0.1 annexation parameter. However, with 
annexation parameters of 0.01 and 0.001, adjoining homogeneous cells 
are frequently annexed into fields reducing the number of passes 
through the classification equation and reducing the amount of computer 
time required. 

Figure 2a~4 graphs the field center pixel performance for the 
Nonsupervised ECHO classifier. The only significant effect is that 
as the annexation parameter goes from 10 - 1 to 10 ~ 3 , the field center 
pixel performance declines. 

For training field and full field performance measures presented 
in Figures 2a-5 and 2a-6, as the annexation threshold becomes smaller 
the performance decreases and as the homogeneity value becomes 
larger the effects of the annexation parameters are increased. 

Parameter selection seems to have an only slightly significant 
effect on the RMS proportion estimate error (Figure 2a-7) . As the 
annexation parameter becomes smaller, the RMS error tends to become 
larger. 

As Figure 2a-8 indicates, classification variability is strongly 
influenced by both the annexation and the homogeneity parameters and 
by their interaction. 
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ii CITARS Results 


Figures 2a-9 through 2a-14 present plots for the 2 by 2 cell size 
Nonsupervised ECHO results achieved over the five CITARS data sets. 
Characteristic differences between the CITARS and the LACIE/SRS 
data sets include: 

* The CITARS data set has a much smaller field size than the 
LACIE/SRS set. 

* The information classes are different. CITARS information 
classes are corn, soybeans, and other; LACIE/SRS classes 
are wheat and other. 

* the ground truth proportion estimates for the LACIE/SRS 
sites were for the area of the LACIE segment. The ground 
truth proportion estimates for the CITARS sites are for 
the whole county in which the data set lies, not for the 
area of the county which was actually sampled. Analysis 
of various results for the five CITARS data sets are 
presented in Table 2a-4. 

Figure 2a-9 indicates that the CITARS data has the same interaction 
of annexation and homogeneity parameters as the LACIE/SRS data sets 
(Figure 2a-3) had with respect to CPU time required. This interaction 
is due to: 


* the additional CPU time required to replace pixel values 
with field means as additional fields are identified, and 


* the reduction of CPU time required by the classifier as 
larger fields are identified, resulting in fewer passes 
through the classification equation. 

As the annexation thresholds go from 10"1 to lO - ^, the CPU time 
required to perform the Nonsupervised ECHO classification is reduced. 
For the CITARS data, as the homogeneity parameter increases, the 
decrease in average CPU time required is statistically significant. 


There are no statistically significant parameter effects on the 
CITARS field center pixel or full field performance measures. Training 
field performance decreases as the annexation thresholds are decreased 
from 10 - ^ to 10” 3 (Figure 2a-12) . As the homogeneity parameter 
increases, the effect of the annexation parameters is increased for 
training field performance, RMS proportion error and classification 
variability measures. 


As the homogeneity parameter increases, RMS proportion estimate 
error increases (Figure 2a-13) and classification variability 
decreases (Figure 2a-14) . As the annexation thresholds decrease, the 
training field performance decreases, the RMS proportion error 
increases, and the classification variability decreases. 
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Figure 2a-9 

Nonsupervised ECHO CPU Requirements 
for the CITARS Data Sets 
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Figure 2a-10 

Nonsupervised ECHO Field Center Pixel Performance 
for the CITARS Data Sets 
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Figure 2a-ll 

Nonsupervised ECHO Full Field Performance 
for the CITARS Data Sets 
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Figure 2a-12 

Nonsupervised ECHO Training Field Performance 
for the CITARS Data Sets 



SM-V1 


graph of kms erh fqh 9 areaisj 
FA-V HU-lll LE-1V 

CELL SUE IS 2X2 


KKS LR» 
40.29 

If. 49 
14.69 
11 . 68 
29.0ft ■ 
26.2 8 
21.4/ 
20.61 • 
I T.HT ■ 
IS. 06 • 
12.26 ■ 


O.OSO 


0. too 


PPPPPPPPPPPPPPPPPPP^PPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPMHPP 


Figure 2a-13 

Nonsupervised ECHO RMS Proportion Estimate Error 
for the CITARS Data Sets 
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Figure 2a-14 


Nonsupervised ECHO Classification Variability 
for the CITARS Data Sets 



2a-27 


Table 2a-4 

Overall Landsat Comparisons 
Effects of Parameters on CITARS Data 


Variable 

Homogeneity 

Thresholds 

Annexation 

Thresholds 

Homogeneity X 
Annexation 

Degrees of Freedom 

2,8 

2,8 

4,16 

CPU Time 

7.57* 

13.73** 

31. 73*** 

Field Center Pixel 
Performance 

.18 

1.48 

. 16 

Full Field 
Performance 

.63 

.84 

.31 

Training Field 
Performance 

4 . 18 + 

8.62* 

3.00* 

RMS Proportion 
Error 

10.29** 

4.81* 

4.22* 

Classification 

Variability 

158.47** 

3. 74 + 

7. 81*** 


Significance 

Levels 

Table entries are F-values 

+ 

10% 


* 

5% 


* * 

1% 


* * * 

. 1% 
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In both the CITARS and the LACIE/SRS data sets, the classifi- 
cation performance measurements are markedly inferior to the perpoint 
performance over the same areas. This circumstance probably results 
from errors in object identification. It Is of interest that while 
Nonsupervised ECHO performance measures for the LACIE/SRS data sets 
are 4 to 10 percentage points lower than the perpoint results, the 
average performance measures for the CITARS data sets are 10 to 30 
percentage points inferior. This difference may be due to the effects 
of the smaller average field size for the CITARS data. The RMS 
proportion error for the Nonsupervised algorithm's results is superior 
to the perpoint algorithm results in the LACIE/SRS data, but inferior 
in CITARS case. 

The Nonsupervised ECHO field center pixel performances are 
superior to the perpoint results in three of the five data sets 
(LE-IV, Hu-III, and Fa-V) by a few percentage points, but on the 
order of 17 to 20 points worse in the other two data sets (Sh-V 
and Wh-V) . This indicates that the Nonsupervised processor has the 
potential to improve the classification accuracy achieved, but that 
it also may substantially decrease the accuracy. 
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iii Overall Landsat Results 

Figures 2a-15 through 2a-20 summarize the effects of the ECHO 
parameter settings on the six dependent variables for nine Landsat 
data sets. Table 2a-5 summarizes the statistical significance of the 
effects of the ECHO parameter settings on these six dependent 
variables . 

The homogeneity parameter has statistically significant effects 
on CPU time, training field performance, RMS proportion error, and 
classification variability. As the cell homogeneity threshold is 
increased, the amount of CPU time required to perform a classification 
decreases (Figure 2a-15) . As the homogeneity parameter increases: 

* the training field performance increases (Figure 2a-18) , 

* RMS proportion estimate error (Figure 2a-19) increases, and 

* classification variability (Figure 2a-20) decreases. 

Parameter settings of the annexation thresholds have a significant 
effect on five dependent variables, CPU time, full field performance, 
training field performance, RMS proportion error, and classification 
variability. As the annexation threshold goes from 10“1 to 10 - ^ 

(the tendency for cell-to-cell annexation increases), the amount of 
CPU time, the full field performance, the training field performance, 
and the classification variability all decrease. As the annexation 
thresholds decrease, the RMS proportion estimate error increases. 

The interaction effects of the homogeneity and annexation 
parameters were significant for four of the dependent variables, CPU 
time, training field performance, RMS proportion error, and classifi- 
cation variability. For these four variables, as the cell homogeneity 
parameter was increased, the effects of the annexation parameters were 
increased. This result is expected since annexation may take place 
only when adjoining cells are homogeneous. 

Table 2a-6 presents the significant effects of the homogeneity and 
annexation parameters on the six dependent variables for each of the 
nine individual Landsat data sets. Although the effects of the homo- 
geneity threshold on the field center pixel and full field performances 
are not significant when the data sets are considered together, the 
effect of the homogeneity parameter is significant in seven and eight 
of the nine Landsat data sets, respectively, when the data sets are 
analyzed individually. This result indicates that the effects of the 
homogeneity parameter on the field center pixel and full field 
performances are opposite for differing data sets. 

On the other hand, the effects of the homogeneity and the 
annexation thresholds on CPU time are not statistically significant 
for any Landsat data set, when the data sets are examined individually. 
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Figure 2a-15 

Nonsupervised ECHO CPU Requirements 
for all Landsat Data Sets 
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Figure 2a-16 

Nonsupervised ECHO Field Center Pixel Performance 
for all Landsat Data Sets 
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Figure 2a-17 

Nonsupervised ECHO Full Field Performance 
for all Landsat Data Sets 


GRAPH OF TRAINING FOR 1 ARCAISI 

FA-V HU— III LE-IV SH-Vl WM-V 


GRAHAM GRANT HASKELL KEARNY 


CELL SUE IS 2X2 


TRAINING 0.050 


pppppppppppppppppppppppppppppppppppppp Pppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppp 


Figure 2a-18 

Nonsupervised ECHO Training Field Performance 
for all Landsat Data Sets 
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Figure 2a- 19 

Nonsupervised ECHO RMS Proportion Estimate Error 
for all Landsat Data Sets 
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Figure 2a-20 

Nonsupervised ECHO Classification Variability 
for all Landsat Data Sets 
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TABLE 2a-5 

Overall Landsat Comparisons 
Effects of Parameters 



Homogeneity 

Annexation 

Homogeneity X 

Variable 

Thresholds 

Thresholds 

Annexation 

Degrees of Freedom 

2,16 

2,16 

4,32 

CPU Time 

7.61** 

13.89*** 

19.08*** 

Field Center Pixel , 

Performance 

.45 

2.92 

.90 

Full Field 

Performance 

1.11 

3.89* 

2.06 

Training Field 

Performance 

4.56* 

9.46** 

7.54*** 

RMS Proportion 

Error 

6.48** 

8.26** 

3.68** 

Classification 

Variability 

43.08*** 

8.87** 

8.96*** 

Significance 

Levels 

Table 

entries are F-values 


+ 10 % 

* 5% 

** 1 % 

*** . 1 % 
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Table 2a-6 

Results of Individual Landsat Data Set Parameter Anovas 
(F Tests with 2,4 degrees of Freedom) 



CPU 

FCP 

FULL 

TRAIN 

RMS 

VARI 

GRAHAM 

Sell 

ns 

ns 

10% 

ns 

ns 

0.1% 

Annl 

10% 

ns 

ns 

ns 

25% 

ns 

GRANT 

Sell 

ns 

5% 

5% 

ns 

0.1% 

0.1% 

Annl 

25% 

10% 

10% 

ns 

ns 

5% 

HASKELL 

Sell 

ns 

25% 

25% 

ns 

25% 

0.1% 

Annl 

10% 

ns 

ns 

ns 

ns 

10% 

KEARNY 

Sell 

ns 

0.1% 

0.1% 

0.1% 

1% 

0.1% 

Annl 

25% 

ns 

ns 

25% 

ns 

25% 

FA-V 

Sell 

ns 

1% 

1% 

ns 

0.1% 

0.1% 

Annl 

ns 

5% 

10% 

25% 

25% 

25% 

HU-III 

Sell 

25% 

0.1% 

1% 

10% 

0.1% 

0.1% 

Annl 

25% 

ns 

ns 

ns 

25% 

ns 

LE-IV 

Sell 

ns 

0.1% 

0.1% 

25% 

0.1% 

0.1% 

Annl 

ns 

ns 

ns 

ns 

10% 

25% 

SH-VI 

Sell 

ns 

0.1% 

0.1% 

1% 

0.1% 

0.1% 

Annl 

25% 

25% 

25% 

25% 

25% 

ns 

WH-V 

Sell 

ns 

1% 

0.1% 

0.1% 

0.1% 

0.1% 

Annl 

ns 

25% 

25% 

1% 

25% 

25% 


Annl Annexation Thresholds 

Sell Homogeneity Thresholds 
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However, the effects are consistent and the cumulative effects of 
both the annexation and the homogeneity parameters are statistically 
significant when the data sets are considered together. 
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iv Classifier Comparisons 

Differences in the results achieved by the Supervised ECHO, 
the Nonsupervised ECHO, and the perpoint classifiers over the Landsat 
data sets will be examined in two ways. First, the classifiers 
will be compared by examining the results of all 768 observations of 
the Supervised ECHO classifier, versus all 81 observations of the 
Nonsupervised ECHO classifier, versus 10 observations of the perpoint 
classifier over the Landsat data. Then, the optimal results for each 
classifier, on each data set will be compared. 

Comparison of all Landsat Observations 

Table 2a-7 presents the effects of classifier on each of the six 
dependent variables. As can be seen from this table, the choice of 
classifier has a significant statistical effect on each of the six 
dependent variables. Table 2a-8 summarizes the results of pair wise 
comparisons of the classifiers. By examining this table, we may 
identify statistically significant differences between specific pairs 
of classifiers. By examining the average response over all ECHO 
parameters settings measured for the three classifiers, the following 
conclusions may be drawn: 

* The Supervised ECHO classifier requires significantly less 
computer time than either the perpoint or Nonsupervised 
ECHO classifiers. There is not, however, a significant 
difference between the computer time required to classify 
an area by a perpoint classifier and the computer time 
required to classify the same area by the Nonsupervised 
ECHO classifier. 

* The field center pixel performance of the Nonsupervised 
ECHO classifier is significantly lower than the field 
center pixel performance of the perpoint classifier and 
the field center pixel performance of the Supervised 
ECHO classifier, when all parameters settings are 
considered. The field center pixel performance achieved 
by the Supervised ECHO classifier is not significantly 
differerit from the field center pixel performance of the 
perpoint classifier when the results of the 24 Supervised 
ECHO classifications produced with the differing parameters 
settings for each data set are averaged. 

* The full field performance of the Nonsupervised classifier 
is significantly lower than the full field performance 

of the Supervised ECHO and perpoint classifiers. There is 
no statistical difference between the performance of the 
Supervised ECHO classifier and the perpoint classifier 
when the results from all 24 ECHO parameter settings of the 
Supervised ECHO processor are averaged and compared to the 
perpoint result. 
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Table 2a-7 


Overall Landsat Comparisons: 
Comparison of Perpolnt, Supervised ECHO, 
and Nonsupervlsed ECHO Classifiers 


Variable 

Degrees 

of 

Freedom 

F-Value 

CPU Time 

2,16 

21.64*** 

FCP Performance 

2,16 

14.41*** 

Full Field Performance 

2,16 

12.85*** 

Training Performance 

2,16 

6.63** 

RMS Error 

2,16 

4.16* 

Classification Variability 

2,16 

10/54** 


Significance Level 

* 5% 
** 1 % 
*** . 1 % 
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Table 2a-8 


Overall Landsat Comparisons: 

Comparison of Perpoint, Supervised ECHO, and 
Nonsupervised ECHO Classifiers 
(859 Observations) 

1 = Supervised ECHO (768 Observations) 

2 = Nonsupervised ECHO ( 81 Observations) 

3 = Perpoint ( 10 Observations) 


Range Tests 

(T-tests with 17 degrees of Freedom) 


CPU Time 1 3 2 

1,2 6.61*** 

1.3 2.11* 

2.3 .30 

Field Center Pixel Performance 2 3 1^ 

1.2 6.59*** 

1.3 .04 

2.3 2.26* 

Full Field Performance 2 l 3 

1,2 6 . 12 *** 

1.3 .06 

2.3 2.19* 

Training Field Performance 2 3 1 

1.2 4.25*** 

1.3 .01 

2.3 1.47 

RMS Proportion Error 3 1 2 

1.2 3.35** 

1.3 .13 

2.3 1.29 

Classification Variability 1 2 3 

1.2 1.67 

1.3 4.99*** 

2.3 4.17*** 
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* The training field performance of the Supervised ECHO 
classifier is significantly higher than the training field 
performance of the Nonsupervised ECHO classifier. The 
training field performance of the perpoint classifier is 
not significantly different from the training field 
performance of the Nonsupervised ECHO classifier; nor is 
the training field performance perpoint classifier 
significantly different from the training field perfor- 
mance of the Supervised ECHO classifier. 

* The RMS proportion -estimate error of the Nonsupervised 
ECHO classifier is significantly larger than the RMS 
proportion estimate error of either the perpoint or the 
Supervised ECHO classifiers. Again, the RMS proportion 
errors for the perpoint and the Supervised ECHO 
classification results are not significantly different. 

* Both the Nonsupervised and the Supervised ECHO classifiers 
have significantly lower variability of classification 
results than the perpoint classifier. The variability of 
the classification results for the two ECHO algorithms 

are not significantly different, however. 

These results indicate that; with random selection of homogeneity 
and annexation parameters, both ECHO classifiers perform better than 
the perpoint classifier only with respect to classification varia- 
bility, and the Supervised ECHO classifier will require less time 
than the perpoint or the Nonsupervised ECHO classifier. On the other 
hand, with random selection of parameter settings, the Nonsupervised 
ECHO classifier has inferior performance with respect to field center 
pixel accuracy, full field accuracy and RMS proportion error to those 
of the Supervised ECHO and the perpoint classifiers. 

Comparison of the Optimal Landsat Results 

Table 2a-9 summarizes the optimal results for each dependent 
variable of each classifier for each Landsat data set. In addition, 
this table lists the difference between each pair of classifiers 
to be examined. The difference between the conclusions reached in 
this section and the conclusions reached in the previous section is 
that, in the previous section, the results for all the homogeneity and 
annexation parameter settings were considered for the Nonsupervised 
and Supervised ECHO classifiers and these classifiers were then 
compared to each other and the perpoint classifier; in this section 
only the results yielded by the optimal parameter settings of the 
Supervised and Nonsupervised ECHO classifiers are considered. 

Table 2a-10 summarizes the results of paired T-tests between the 
optimal performances of each pair of classifiers over the Landsat 
data sets. 



3a-40 


TABLE 2a-9 



OPTIMAL 

RESULTS 

OF LANDSAT 

DATA SETS 



Data 

Set 

Variable Perpoint 

Super. 

Nonsup . 

Perpoint 

-Sup. 

Perpoint 

-Nonsup 

Nonsup 
-Sup . 


GRAHAM 

CPU 

121.00 

74.00 

154.00 

47.00 

-33.00 

80.00 


FCP 

89.90 

93.90 

95.10 

-4.00 

-5.20 

1.20 


FF 

87.40 

90.00 

92.00 

-2.60 

-4.60 

2.00 


TRAIN 

90.40 

98.50 

94.50 

-8.10 

-4.10 

-4.00 


RMS 

1.20 

.10 

.10 

1.10 

1.10 

0.00 


VAR 

20.85 

10.80 

10.49 

10.05 

10.36 

-.31 

GRANT 

CPU 

122.00 

98.00 

149.00 

24.00 

-27.00 

51.00 


FCP 

64.50 

68.10 

60.80 

-3.60 

3.70 

-7.30 


FF 

66.40 

70.80 

64.90 

-4.40 

1.50 

-5.90 


TRAIN 

56.20 

58.40 

54.00 

-2.20 

2.20 

-4.40 


RMS 

6.50 

.10 

2.28 

-1.60 

4.22 

-5.82 


VAR 

24.62 

8.60 

12.81 

16.02 

11.82 

4.20 

KEARNY 

CPU 

389.00 

100.51 

205.00 

288.49 

184.00 

104.49 


FCP 

65.00 

65.40 

51.00 

-0.40 

14.00 

-14.40 


FF 

63.90 

64.40 

52.30 

-0.50 

11.60 

-12.10 


TRAIN 

67.60 

68.20 

45.20 

-0.60 

22.40 

-23.00 


RMS 

12.30 

12.28 

8.00 

.02 

4.30 

-4.28 


VAR 

44.12 

16.16 

11.36 

27.96 

32.76 

-4.80 

HASKELL 

CPU 

396.00 

166.82 

223.00 

229.18 

173.00 

56.18 


FCP 

67.30 

77.20 

65.40 

-9.90 

1.90 

-11.80 


FF 

67.00 

76.40 

64.50 

-9.40 

2.50 

-11.90 


TRAIN 

77.70 

88.50 

63.80 

-10.80 

13.90 

-24.70 


RMS 

3.50 

.98 

.10 

2.52 

3.40 

-0.88 


VAR 

45.73 

14.07 

13.55 

31.66 

32.18 

-0.52 

LE-IV 

CPU 

284.00 

180.22 

306.00 

103.78 

-22.00 

125.78 


FCP 

50.70 

57.20 

55.10 

-6.50 

-4.40 

-2.10 


FF 

50.60 

55.70 

54.20 

-5.10 

-3.60 

-1.50 


TRAIN 

69.00 

73.00 

60.60 

-4.00 

8.40 

-12.40 


RMS 

16.20 

19.37 

23.' 78 

-3.17 

-7.58 

4.41 


VAR 

46.38 

8.54 

3.29 

37.84 

43.09 

-5.25 

HU-III 

CPU 

852.00 

480.00 

696.00 

372.00 

156.00 

216.00 


FCP 

59.10 

61.70 

61.90 

-2.60 

-2.80 

.20 


FF 

59.50 

59.70 

52.00 

-0.20 

7.50 

-7.70 


TRAIN 

85.20 

90.80 

40.00 

-5.60 

45.20 

-50.80 


RMS 

14.70 

11.59 

6.83 

3.11 

7.87 

-4.76 


VAR 

48.90 

18.90 

17.80 

30.00 

31.10 

-1.10 
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TABLE 2a-9 (cont'd) 
OPTIMAL RESULTS OF LANDSAT DATA SETS 


Data 

Set 

Variable Perpoint 

Super. 

Non sup . 

Perpoint 
-Sup . 

Perpoint 
-Nonsup . 

Nonsup 

-Sup. 


SH-VI 

CPU 

526.00 

296.75 

473.00 

229.25 

53.00 

176.25 


FCP 

66.90 

63.02 

50.60 

3.80 

16.30 

-12.42 


FF 

62.40 

60.20 

50.10 

2.20 

12,30 

-10.10 


TRAIN 

59.20 

65.70 

53.10 

-6.50 

6.10 

-12.60 


RMS 

21.60 

20.50 

19.24 

1.10 

2.36 

-1.26 


VAR 

48.45 

13.09 

9.53 

35.36 

38.92 

-3.56 

WH-V 

CPU 

458.00 

239.13 

425.00 

218.87 

33.00 

185.87 


FCP 

74.60 

75.50 

57.30 

-0.90 

17.30 

-18.20 


FF 

. 70.60 

70.51 

57.60 

.09 

13.00 

-12.91 


TRAIN 

87.80 

91.00 

65.90 

-3.20 

21.90 

-25.10 


RMS 

6.90 

7.20 

23.13 

-0.30 

-16.23 

15.93 


VAR 

47.87 

19.30 

10.75 

28.57 

37.12 

-8.55 

LIVSTON 

CPU 

513.00 

253.00 


260.00 




FCP 

71.30 

75.20 


-3.90 




FF 

68.90 

69.40 


-0.50 




TRAIN 

85.10 

88.20 


-3.10 




RMS 

12.80 

7.00 


5.80 




VAR 

53.07 

18.00 


35.07 



FA-V 

CPU 

348.00 

198.88 

325.50 

149.12 

22.50 

126.62 


FCP 

86.60 

89.20 

89.50 

-2.60 

-2.90 

.30 


FF 

78.40 

80.10 

74.30 

-1.70 

4.10 

-5.80 


TRAIN 

72.20 

78.00 

88.50 

-5.80 

-16.30 

10.50 


RMS 

4.70 

5.03 

16.05 

-0.33 

-11.53 

-11.02 


VAR 

44.30 

18.28 

18.80 

26.02 

25.50 

.52 


CPU times are in seconds. 

Field center pixel (FCP) , full field (FF) , and training field (TRAIN) 
performances, RMS proportion estimate error, and classification 
variabilities (VAR) are in percentages. 
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TABLE 2a-9 (cont'd) 
OPTIMAL RESULTS OF LANDSAT DATA SETS 


Data 




Perpoint 

Perpoint 

Nonsup . 

Set 

Variable Perpoint 

Super. 

Nonsup . 

-Sup . 

-Nonsup . 

-Sup . 


SH-VI 

CPU 

526.00 

296.75 

473.00 

229.25 

53.00 

176.25 


FCP 

66.90 

63.02 

50.60 

3.80 

16.30 

-12.42 


FF 

62.40 

60.20 

50.10 

2.20 

12.30 

-10.10 


TRAIN 

59.20 

65.70 

53.10 

-6.50 

6.10 

-12.60 


RMS 

21.60 

20.50 

19.24 

1.10 

2.36 

-1.26 


VAR 

48.45 

13.09 

9.53 

35.36 

38.92 

-3.56 

WH-V 

CPU 

458.00 

239.13 

425.00 

218.87 

33.00 

185.87 


FCP 

74.60 

75.50 

57.30 

-0.90 

17.30 

-18.20 


FF 

70.60 

70.51 

57.60 

.09 

13.00 

-12.91 


TRAIN 

87.80 

91.00 

65.90 

-3.20 

21.90 

-25.10 


RMS 

6.90 

7.20 

23.13 

-0.30 

-16.23 

15.93 


VAR 

47.87 

19.30 

10.75 

28.57 

37.12 

-8.55 

LIVSTON 

CPU 

513.00 

253.00 


260.00 




FCP 

71.30 

75.20 


-3.90 




FF 

68.90 

69.40 


-0.50 




TRAIN 

85.10 

88.20 


-3.10 




RMS 

12.80 

7.00 


5.80 




VAR 

53.07 

18.00 


35.07 



FA-V 

CPU 

348.00 

198.88 

325.50 

149.12 

22.50 

126.62 


FCP 

• 86.60 

89.20 

89.50 

-2.60 

-2.90 

.30 


FF 

78.40 

80.10 

74.30 

-1.70 

4.10 

-5.80 


TRAIN 

72.20 

78.00 

88.50 

-5.80 

-16.30 

10.50 


RMS 

4.70 

5.03 

16.05 

-0.33 

-11.53 

-11.02 


VAR 

44.30 

18.28 

18.80 

26.02 

25.50 

.52 


CPU times are in seconds. 

Field center pixel (FCP), full field (FF), and training field (TRAIN) 
performances, RMS proportion estimate error, and classification 
variabilities (VAR) are in percentages. 
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Table 2a-10 

Comparison of the Optimal Landsat Results for 
the Supervised ECHO, Nonsupervised ECHO and Perpoint Classifier 



Perpoint 

Perpoint 

Nonsupervised 


versus 

versus 

versus 

Variable 

Supervised 

Nonsupervised 

Supervised 

Observations 

10 

9 

9 

CPU 

T 

5.53 

2.03 

6.41 

Significance Level 

.1% 

5% 

.1% 

FCP 

T 

2.65 

1.37 

2.94 

Significance Level 

1% 

NS 

1% 

FF 

T 

2.1 

2.22 

4.31 

Significance Level 

5% 

5% 

1% 

Training 

T 

5.24 

1.88 

2.81 

Significance Level 

.1% 

5% 

5% 

RMS Error 

T 

1.03 

.48 

.33 

Significance Level 

NS 

NS 

NS 

Classification Variability 

T 

10.08 

7.67 

1.71 

Significance Level 

.1% 

.1% 

NS 
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At optimal parameter settings, both the Nonsupervised and 
Supervised ECHO classifiers require significantly less CPU time 
than the perpoint classifier. This is true even though the 
Nonsupervised classifier was run with the "MAP" option specified 
for the Landsat data. In addition, the Supervised ECHO classifier 
requires significantly less time than the Nonsupervised ECHO 
classifier for these data sets. The specification of the "MAP" 
option for the Nonsupervised classifier contributed to this result. 

The Supervised classifier demonstrates a significantly superior 
field center pixel performance to either the perpoint or the Nonsuper- 
vised ECHO classifier. However, the Nonsupervised classifier, which 
identifies objects without the benefit of class statistics, demon- 
strates no statistically significant difference from the perpoint 
classifier with respect to field center pixel accuracy, at optimal 
annexation and homogeneity parameter settings. 

The full field performance of the Supervised ECHO classifier is 
significantly superior to that of the perpoint and Nonsupervised 
ECHO classifiers. On the other hand, the full field performance of 
the Nonsupervised ECHO classifier at its optimal parameter settings 
is still significantly inferior to that of the perpoint classifier. 

For the Landsat data sets, the Supervised ECHO classifier has 
a statistically significant advantage in training performance to both 
the Nonsupervised ECHO and the perpoint classifiers. Again, even 
at optimal parameter settings the Nonsupervised ECHO classifications 
are at a statistically significant disadvantage when compared to the 
perpoint results. 

There are no statistically significant differences between the 
perpoint, the Supervised ECHO and the Nonsupervised ECHO classifiers 
with respectto RMS proportion estimate errors. 

While not being significantly different from each other, the 
results of both the Supervised and the Nonsupervised ECHO classifiers 
at their optimal parameter settings demonstrate significantly less 
classification variability than the perpoint results. 
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b. Simulated Thematic Mapper Results 

The analyses of variance for the effects of the Nonsupervised 
ECHO parameters on the six simulated Thematic Mapper data sets (two 
resolutions from Williams County, ND, four resolutions from Finney 
County, Kansas) are presented in Table 2a-ll. Results for the six 
data sets considered together are presented in Table 2a-12. 

i Cell Width Parameter 

For the Landsat and aircraft data sets the Nonsupervised ECHO 
cell width parameter was set to two in all cases. For those data 
sets the cell width setting of two was selected to economize on the 
CPU time required to perform the verification tests based on the ob- 
servation that, for Supervised ECHO results, larger cell width parame- 
ters did not improve ECHO performances. The six simulated Thematic 
Mapper data sets were classified by the Nonsupervised ECHO classifier 
at cell widths of two and three. There were several statistically 
significant effects of the cell width parameter on the Nonsupervised 
ECHO classifications of simulated Thematic Mapper data: 

* When the cell width parameter is set to three, less CPU time 
is required than when it is set to two (significant at a 10% 
confidence level). 

* The training performance at cell width three is significantly 
poorer than at cell width two (5% confidence level) . 

* The RMS proportion estimate error is less at cell width three 
than at cell width two (10% confidence level). 

* The classification variability is less for cell width two than 
cell width three (10% confidence level). 

The fact that the proportion estimates are better for cell width 
three than for cell width two for the simulated Thematic Mapper data 
indicates that the partitioning using the larger cell size may, indeed, 
be of value with the Nonsupervised ECHO processor. 

Because the likelihood that a cell will be homogeneous when it 
contains nine scene elements is smaller than the likelihood that a 
cell containing only four scene elements will be homogeneous, it 
makes sense that, for a given homogeneity setting, more cells at cell 
width three, than at cell width two, would be identified as singular 
and split with constituent points classified individually. This 
situation may explain why the classification variability at cell width 
two is smaller than the classification variability at cell width three. 
On the other hand, when large homogeneous areas are present in a scene, 
fewer calculations of the classification equation will be necessary 
at larger cell widths. Thirty-six pixels in 2 by 2 pixel cells require 
nine calculations of class probability for each spectral class while, 
in 3 by 3 pixel cells, only four calculations of class probability 
for each spectral class are required. 
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Table 2a- 11 

Results of Individual Thematic Mapper Data Set 
Parameter Anovas 
(108 Observations) 



CPU 

FCP 

FULL 

TRAIN 

RMS 

VARI 

DF 

1730 








Celw 

1% 

ns 

5% 

25% 

1% 

1% 

1,4 

Sell 

0.1% 

5% 

0.1% 

1% 

1% 

0.1% 

2,4 

Annl 

1% 

25% 

25% 

ns 

10% 

0.1% 

2,4 

Celw X Sell 

1% 

ns 

10% 

10% 

ns 

0.1% 

2,4 

Celw X Annl 

1% 

ns 

ns 

ns 

25% 

5% 

2,4 

Sell X Annl 

25% 

ns 

ns 

ns 

10% 

5% 

4,4 

1740 








Celw 

0.1% 

5% 

10% 

ns 

1% 

1% 

1,4 

Sell 

0.1% 

5% 

1% 

1% 

0.1% 

0.1% 

2,4 

Annl 

1% 

25% 

25% 

5% 

1% 

0.1% 

4,4 

Celw X Sell 

5% 

ns 

ns 

ns 

5% 

0.1% 

2,4 

Celw X Annl 

5% 

ns 

ns 

25% 

ns 

5% 

2,4 

Sell X Annl 

5% 

25% 

25% 

25% 

1% 

1% 

4,4 

3730 








Celw 

0.1% 

ns 

ns 

ns 

ns 

ns 

1,4 

Sell 

0.1% 

0.1% 

0.1% 

0.1% 

0.1% 

0.1% 

2,4 

Annl 

0.1% 

5% 

5% 

25% 

ns 

5% 

4,4 

Celw X Sell 

1% 

ns 

ns 

ns 

10% 

5% 

2,4 

Celw X Annl 

1% 

25% 

10% 

25% 

25% 

25% 

2,4 

Sell X Annl 

5% 

5% 

5% 

25% 

ns 

25% 

4,4 

3740 








Celw 

1% 

ns 

10% 

ns 

25% 

5% 

1,4 

Sell 

0.1% 

ns 

0.1% 

0.1% 

0.1% 

0.1% 

2,4 

Annl 

0.1% 

ns 

5% 

ns 

25% 

1% 

2,4 

Celw X Sell 

1% 

ns 

1% 

ns 

0.1% 

1% 

2,4 

Celw X Annl 

1% 

ns 

ns 

ns 

ns 

5% 

2,4 

Sell X Annl 

5% 

ns 

25% 

ns 

ns 

10% 

4,4 

3750 








Celw 

5% 

0.1% 

0.1% 

0.1% 

0 . 1% 

0.1% 

1,4 

Sell 

0.1% 

0.1% 

0.1% 

0.1% 

0.1% 

0.1% 

2,4 

Annl 

1% 

ns 

25% 

5% 

10% 

10% 

2,4 

Celw X Sell 

5% 

1% 

5% 

5% 

5% 

1% 

2,4 

Celw X Annl 

1% 

ns 

ns 

ns 

ns 

5% 

2,4 

Sell X Annl 

10% 

25% 

25% 

10% 

25% 

25% 

4,4 
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Table 2a- 11 (Continued) 


CPU 


3760 

Celw 25% 

Sell 0.1% 

Annl 1% 

Celw X Sell 5% 

Celw X Annl 5% 

Sell X Annl 5% 


FCP 

FULL 

TRAIN 

ns 

ns 

ns 

1% 

1% 

0.1% 

ns 

ns 

ns 

ns 

ns 

ns 

25% 

25% 

25% 

ns 

ns 

ns 


RMS 

VARI 

DF 

10% 

ns 

1,4 

1% 

5% 

2,4 

ns 

ns 

2,4 

ns 

ns 

2,4 

ns 

ns 

2,4 

ns 

ns 

. 4,4 


o 
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Table 2a- 12 


Parameters 


OVERALL THEMATIC MAPPER RESULTS 
Significant Effects of Parameters on Dependent Variables 

(108 Observations) 


Degrees 

of 

Freedom 


Variables 


Resolution 

3,1 

Cell Width 

1,1 

Homogeneity 

Threshold 

2,2 

Annexation 

Threshold 

2,2 

Resolution X 
Cell Width 

3,1 

Resolution X 
Homogeneity 

6,2 

Resolution X 
Annexation 

6,2 

Cell Width X 
Homogeneity 

2,2 

Cell Width X 
Annexation 

2,2 

Homogeneity X 
Annexation 

4,4 

Significance 

: Level 

+ 

10% 

* 

5% 

** 

1% 

*** 

.1% 


CPU FCP 


Time 

Perform. 

8.35 

57.904- 

91.73+ 

3.34 

25.51* 

32.50* 

178.21** 

4.05 

92.01+ 

14.51 

3.68 

50.21* 

.68 

1.72 

45.59* 

479.36** 

30.27* 

19.66* 

227.36*** 

2.89 


Full Field 
Perform. 

Training 

Perform. 

19.44 

509.24* 

0.76 

345.90* 

16.26+ 

164.70** 

3.45 

13.46+ 

4440.00* 

64.71+ 

22.59* 

12.70+ 

3.05 

.31 

7.68 

44.74* 

5.74 

11.05+ 

3.2 

29.19** 


RMS 

Error 

Classi. 

Vari. 

5.41 

2.85 

119.46+ 

72.44+ 

68.63* 

73.41* 

9.78+ 

65.95* 

19.90 

38.07 

43.59* 

48.58* 

8.07 

1.87 

5.83 

12.87+ 

286.00** 

63.04** 


4.21+ 96.25*** 


Table entries are for F-values 
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There are several interaction effects of cell width and resolution. 

As resolution increases, the effect of the cell width parameter on CPU 
time decreases (10% confidence level) and the effects of the cell width 
parameter on the full field and the training performances decrease (10% 
and 5% confidence levels, respectively). As resolution elements become 
larger, fields will be sampled by fewer pixels. Therefore, the number 
of truly homogeneous cells at a given cell size will decrease. Since 
there will be fewer homogeneous cells at 60 meter resolution than at 30 
meter resolution, more cell splitting will take place and, hence, the 
effects of aggregating pixels into cells will be smaller at 60 meter 
resolution than at 30 meter resolution. This effect expands as the cell 
size becomes larger since the number of homogeneous 3 by 3 cells will 
shrink faster than the number of homogeneous 2 by 2 cells as the resolu- 
tion elements become larger. 

The interaction of the cell width parameter and the homogeneity 
parameter is statistically significant with respect to CPU time, field center 
pixel performance, training field performance, and variability of classi- 
fication results. In all cases, as the cell width parameter increases, 
the effect of the homogeneity parameter increases. This results is logical 
since, as the cell become larger, the number of pixels affected by the out- 
come of each homogeneity test is larger. 

There are also significant cell width-annexation threshold interactions. 
As the cell width becomes larger, the effect of the annexation parameter 
becomes smaller for CPU time, training field performance, field center pixel 
performance, RMS proportion estimate error, and classification variability. 
This results reflects the facts that as the partition becomes larger, there 
will be fewer and fewer homogeneous cells, and that, as the partition becomes 
larger, it becomes less likely that a neighboring homogeneous cell will be 
contained in a single object. 
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ii Resolution 

As resolution increases, both field center pixel and traiing field 
performances show a statistically significant decrease. The 50 meter 
Thematic Mapper data is probably responsible for this circumstance. Both 
the field center pixel and the training field performances on the 50 meter 
data sets are significantly below the performance levels of the 40 and 
the 60 meter data sets. Reasons for this degradation in performance at 
50 meters are unclear. One would expect a tailing off of classification 
accuracy as the size of the cell approaches the size of the objects on the 
ground. However, in this data set, even at 60 meter resolution, cell sizes 
do not approach the size of the agricultural fields. 

The interaction of the resolution with the homogeneity parameter is 
statistically significant for the field center pixel performance, the full 
field performance, the training performance, the RMS proportion estimate 
error, and the classification variability. For the 30 meter resolution the 
field center pixel and the full field performances improve as the cell homo- 
geneity parameter increases from .05 to .10, and then fall somewhat as the 
homogeneity parameter goes from .10 to .25. At the larger resolutions, 
however, these variables fall steadily as the homogeneity parameter is in- 
creased. The effects of the homogeneity parameter increase for: the RMS 

proportion estimate error, the training field performance, and the classi- 
fication variability as the resolution element size increases. 

Interaction of the resolution and the annexation thresholds is not 
significant for any of the six dependent variables monitored. 
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iii Homogeneity Threshold 

The homogeneity threshold has a significant effect on all six depen- 
dent variables. As the homogeneity threshold increases: 

* the CPU time required to perform the Nonsupervised ECHO classifi- 
cation decreases, 

* the field center pixel performance decreases, 

* the full field performance decreases, 

* the training field performance decreases, 

* the RMS proportion estimate error increases, and 

* the classification variability decreases. 

Figures 2a-21 through 2a-26 graph the effects of the homogeneity and 
annexation parameters on the six dependent variables for 2 by 2 pixel cells. 
Unlike the Landsat results (See figure 2a-15), for the simulated Thematic 
Mapper data sets, the effect of homogeneity parameter on CPU time is statis- 
tically significant. The "MAP" option was not used on the simulated Thematic 
Mapper data. Therefore, there was no data replacement when homogeneous 
objects were identified in the field extraction phase, and consequently, 
the CPU time required goes down as the number of homogeneous cells increases 
due to the reduction in the number of passes through the classification 
equation made possible by classifying a field, rather than a point at a 
t ime . 


The results of the field center pixel performance, full field perfor- 
mance, training field performance, and RMS proportion estimate error measure- 
ments indicate that for simulated Thematic Mapper data, at resolutions above 
30 meters, homogeneity parameter specifications of .05 or less are appropriate. 

Interaction between the homogeneity parameter and the annexation 
thresholds is significant for CPU time, training field performance, RMS 
proportion estimate error, and classification variability. As the homo- 
geneity parameter increases, the effects of the annexation also increase for 
each of the dependent variables listed above. Since annexation can take 
place only when adjoining cells are homogeneous and since the number of 
homogeneous fields increases as the homogeneity parameter increases, it makes 
sense for the annexation parameter to have a greater effect when the homo- 
geneity parameter is larger. 
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Figure 2a-21 

Nonsupervised ECHO CPU Requirements 

for the Simulated Thematic Mapper Data Sets 
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Figure 2a-22 


Nonsupervised ECHO Field Center Pixel Performance 
for the Simulated Thematic Mapper Data Sets 
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Figure 2a-23 

Nonsupervised ECHO Full Field Performance 
for the Simulated Thematic Mapper Data Sets 
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Figure 2a-24 

Nonsupervised ECHO Training Field Performance 
for the Simulated Thematic Mapper Data Sets 



GRAPH OF RMS ERR FOR 6 AREA I S I 

17 50 1740 3730 3740 

CELL SUE IS 2X2 


RMS ERR 0.0*0 
0. 78 

a. i? 


ppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppprpppppppppppppppppppppppppp 


Figure 2a-25 

Nonsupervised ECHO RMS Proportion Estimate Error 
for the Simulated Thematic Mapper Data Sets 


GRAPH OF VAHHlfV FOR 6 ART.AISl 

17 10 1 74.' 1710 

CELL SI/1 IS / X / 


V AHiit TV 0.0*0 


PPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPP 


Figure 2a-26 

Nonsupervised ECHO Classification Variability 
for the Simulated Thematic Mapper Data Sets 
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iv Annexation Threshold 

The annexation threshold has a significant effect on the CPU time 
required to perform classification, on the training performance, on the 
RMS proportion estimate error, and on the classification variability. As 
the annexation threshold goes from 1 0 — 1 to 10“3 } the time required to per- 
form a classification decreases, the training performance decreases, the 
RMS proportion estimate error increases, and the classification variability 
decreases. This result indicates that for the simulated Thematic Mapper 
data, annexation parameters in the neighborhood of 10” are most appropriate. 


v Classifier Comparisons 

Differences in the results achieved by the Supervised ECHO, the Non- 
supervised ECHO, and the perpoint classifiers over the simulated Thematic 
Mapper data sets will be examined in two ways. First, the classifiers will 
be compared by examining the results of 8 observations of the perpoint 
classifier versus the 108 observations of the Nonsupervised ECHO processor 
versus the 576 observations of the Supervised ECHO processor over the 
simulated Thematic Mapper data sets. Then, the optimal results for each 
classifier, on each data set will be compared. 

Comparison of all Simulated Thematic Mapper Observations 

Table 2a-13 presents the effects of classifier and resolution on the 
simulated Thematic Mapper data sets for the six dependent variables when 
all observations are considered for the ECHO classifiers. The effect of 
the classifier is significant for all six of these variables. Resolution 
has a statistically significant effect on field center pixel performance, 
full field performance, and classification variability. The interaction 
of the classifier and the resolution is statistically significant for 
training field performance, field center pixel performance, full field 
performance, RMS proportion estimate error, and classification variability. 
Table 2a-14 illustrates some of the effects of the three classifiers. The 
CPU time required by the Nonsupervised classifier is significantly less 
than the CPU time required by either the Supervised ECHO or the perpoint 
classifiers. On the other hand, the Nonsupervised ECHO classifier demon- 
strates inferior field center pixel, full field, and training field perfor- 
mances when all ECHO observations are considered. 

The poor performance of the Nonsupervised ECHO processor on the 50 
meter data may be responsible for many of the significant effects. By 
examining the graphs of the average field center pixel and full field per- 
formance in Table 2a-14, it can be seen that the 50 meter resolution is 
causing the significant interaction between resolution and classifier and 
the statistically significant differences between classifiers, for these 
two variables. The 50 meter results are also causing the significant effect 
of resolution on field center pixel performance, and the interaction effect 
between resolution and classifier for the training field performance and 
classification variability variables. 



Table 2a-13 


Overall Thematic Mapper Results 
692 Observations 


EFFECT 

DEGREES 

OF 

FREEDOM 

CPU 

TIME 

FCP 

FULL 

FIELDS 

TRAINING 

FIELDS 

RMS 

PROPORTION 

CLASSIFICATION 

VARIABILITY 

CLASSIFIER 

2,2 

18.72+ 

229.11** 

122.61** 

514.46** 

80.47** 

69.45* 

RESOLUTION 

3,3 

5.02 

9.51* 

10.35* 

2.06 

2.90 

7.409+ 

CLASSIFIER X 
RESOLUTION 

6,4 

.65 

9.94* 

6.82* 

10.77* 

5.37+ 

7.989* 

Confidence Level: 




Table 

entries are 

F-values 


+ 10 % 

* 5% 

** 1 % 

*** . 1 % 
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Table 2a- 14 

Overall Thematic Mapper 
Effects of Classifier on Variables 


1 = Supervised ECHO Results 

2 = Nonsupervised ECHO Results 

3 = Perpoint Results 


CPU Time 

1,2 

1.3 

2.3 


RMS Proportion Error 1 3 2 


1.2 4.49** 

1.3 .03 

2.3 1.26 


Significance Level Table entries are T-values 

* 5% 

** 1 % 


2 1 3 

3.52* 

2.48 

3.42* 



Table 2a-14 


Field Center Pixel Performance 



Full Field Performance 



Training Field Performance 





Table 2a- 14 (Continued) 


Classification Variability 



Looks like all three 
could be sig. diff. 
Res = 50 removed. 
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As the resolution elements increase in size, the classification varia- 
bility increases. This reflects the fact that it becomes less likely that 
adjacent pixels will be of the same class as the size of resolution elements 
increase. 

There is a statistically significant difference between classifiers 
with respect to classification variability. The perpoint classification 
results are more variable than Supervised ECHO classification results 
which are, in turn, more variable than the Nonsupervised ECHO classifica- 
tion results. The difference between the two ECHO processors indicates 
that the homogeneity parameters used for the Nonsupervised ECHO tests 
result in more cells being identified as homogeneous than the homogeneity 
parameters in the Supervised ECHO experiments. It is possible that for the 
simulated Thematic Mapper data sets, the optimal homogeneity parameter may 
be less than the smallest Nonsupervised ECHO homogeneity parameter tested 
(which was 0.05). 

The RMS proportion estimate error for the Supervised ECHO processor 
is significantly lower than that of the Nonsupervised processor. The pro- 
portion estimate error of the perpoint classifier is significantly dif- 
ferent from neither the Supervised nor the Nonsupervised ECHO processor. 

Comparison of Classifier Optimums 

Table 2a-15 presents the optimal dependent variable measurement (lowest 
CPU, RMS proportion estimate error, classification variability, and highest 
field center pixel, full field and training field performance) for each of 
the three classifiers, together with the differences between pairs of 
optimal responses, for each of the simulated Thematic Mapper data sets. 

Table 2a-16 presents the results of a paired T-test to identify the signifi- 
cant differences between each pair of classifiers with respect to each de- 
pendent variable. 

Comparing optimal results, the ECHO classifiers require significantly 
less CPU time to classify an area than the perpoint classifier requires. 

The CPU time required to produce a classification by the Nonsupervised 
and the Supervised ECHO processors are not statistically different, however. 

At optimal parameter settings, the Supervised classifier has signifi- 
cantly higher field center pixel and full field classification performances 
than either the Nonsupervised ECHO or the perpoint classifier. The Non- 
supervised ECHO classifier is not significantly different from the perpoint 
classifier at optimal parameter settings of the Nonsupervised processor with 
respect to these two variables. 

The Supervised ECHO results have training field performance statistically 
superior to those of the perpoint classifier, which in turn, have the 
training field performance statistically higher than the results of the Non- 
supervised ECHO processor. 

The RMS proportion estimate error of the perpoint classifier was not 
statistically different from either of the ECHO classifiers. The Nonsuper- 
vised ECHO classifier had significantly higher proportion estimate error than 
the Supervised processor, however. 
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Table 2a- 15 




OPTIMAL 

RESULTS OF 

THEMATIC MAPPER 

DATA SETS 



Data 





Perpoint 

Perpoint 

Nonsup . 

Set 

Variable 

Perpoint 

Super. 

Nonsup . 

-Sup. 

-Nonsup . 

-Sup. 

1730 

CPU 

1170.00 

335.20 

370.90 

834.80 

799.10 

35.70 


FCP 

92.10 

94.70 

91.20 

-2.60 

.90 

-3.50 

* 

FF 

88.10 

89.60 

87.40 

-1.50 

.70 

-2.20 


TRAIN 

97.60 

98.50 

96.20 

-0.90 

1.40 

-2.30 


RMS 

3.40 

3.61 

3.77 

-0.21 

-0.37 

.16 


VAR 

35.93 

14.27 

12.66 

21.66 

23.27 

-1.61 

1740 

CPU 

732.00 

228.94 

256.80 

503.06 

475.20 

27.86 


FCP 

88.60 

91.20 

91.70 

-2.60 

-3.10 

.50 


FF 

82.60 

84.20 

85.00 

-1.60 

-2.40 

.80 


TRAIN 

97.20 

97.70 

97.30 

-0.50 

-0.10 

-0.40 


RMS 

8.10 

7.80 

7.98 

.30 

.12 

.18 


VAR 

37.68 

18.31 

13.28 

19.37 

24.40 

-5.03 

1750 

CPU 

436.40 

144.97 


291.43 




FCP 

88.20 

88.60 


-0.40 




FF 

74.70 

74.90 


-0.20 




TRAIN 

98.70 

98.90 


-0.20 




RMS 

7.60 

7.66 


-0.06 




VAR 

35.22 

15.71 


19.51 



1760 

CPU 

358.80 

111.66 


247.14 




FCP 

94.50 

95.50 


-1.00 




FF 

81.00 

81.90 


-0.90 




TRAIN 

98.70 

98.70 


0.00 




RMS 

8.50 

8.60 


-0.10 




VAR 

42.31 

22.30 


20.01 



3730 

CPU 

1200.00 

353.16 

402.50 

846.84 

797.50 

49.34 


FCP 

87.50 

91.40 

89.30 

-3.90 

-1.80 

-2.10 


FF 

84.10 

87.70 

85.80 

-3.60 

-1.70 

-1.90 


TRAIN 

97.90 

98.60 

97.60 

-7.00 

.30 

-1.00 


RMS 

7.30 

2.40 

2.34 

4.90 

4.96 

-0.06 


VAR 

34.75 

18.47 

17.51 

16.28 

17.24 

-0.96 

3740 

CPU 

1222.00 

418.43 

312.00 

803.57 

910.00 

-106.43 


FCP 

91.60 

92.80 

90.50 

-1.20 

1.10 

-2.30 


FF 

86.20 

87.60 

85.70 

-1.40 

.50 

-1.90 


TRAIN 

96.60 

98.00 

95.90 

-1.40 

.70 

-2.10 


RMS 

2.70 

2.71 

3.17 

-0.01 

-0.47 

.46 


VAR 

45.93 

27.02 

24.63 

18.91 

21.30 

-2.39 
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Table 2a-15 (cont'd) 
OPTIMAL RESULTS OF THEMATIC MAPPER DATA SETS 


Data 





Perpoint 

Perpoint 

Nonsup 

Set 

Variable 

Perpoint 

Super. 

Nonsup . 

-Sup . 

-Nonsup . 

-Sup. 

3750 

CPU 

855.20 

233.47 

218.60 

621.73 

636.60 

-14.87 


FCP 

91.50 

91.46 

78.10 

.04 

13.40 

-13.36 


FF 

85.00 

84.68 

75.80 

.32 

9.20 

-8.88 


TRAIN 

92.90 

92.90 

84.90 

0.00 

8.00 

-8.00 


RMS 

3.80 

3.79 

7.69 

.01 

-3.89 

3.90 


VAR 

46.08 

21.20 

14.70 

24.88 

31.38 

-6.50 

3760 

CPU 

779.60 

322.56 

170.40 

457.04 

609.20 

-152.16 


FCP 

93.20 

93.80 

92.10 

-0.60 

1.10 

-1.70 


FF 

84.10 

84.70 

83.30 

-0.60 

.80 

-1.40 


TRAIN 

97.70 

98.20 

96.20 

-0.50 

1.50 

-2.00 


RMS 

3.00 

2.61 

3.53 

.39 

-0.53 

.92 


VAR 

50.35 

33.97 

2.73 

16.38 

47.62 

-31.24 
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Table 2a- 16 


Comparison of the Optimal 

Simulated Thematic Mapper 

Results 

for the Supervised ECHO, Nonsupervised ECHO, 

and Perpoint 

Classifiers 


Perpoint 

Perpoint 

Nonsupervised 


versus 

versus 

versus 

Variable 

Supervised 

Nonsupervised Supervised 

Observations 

CPU 

T 

6.79 

10.87 

1.87 

Significance Level 

.1% 

.1% 

NS 

Field Center Pixel Performance 

T 

3.18 

.80 

1.87 

Significance Level 

1% 

NS 

10% 

Full Field Performance 

T 

2.85 

.70 

1.93 

Significance Level 

1% 

NS 

10% 

Training Field Performance 

T 

1.58 

1.60 

2.36 

Significance Level 

10% 

10% 

5% 

RMS Proportion Estimate Error 

T 

1.08 

.03 

1.52 

Significance Level 

NS 

NS 

10% 

Classification Variability 

T 

19.97 

6.21 

1.68 

Significance Level 

.1% 

1% 

10% 
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The results of the Nonsupervised ECHO processor were significantly 
less variable than those of either the perpoint or the Supervised ECHO 
processor. In addition, the Supervised ECHO classification results were 
less variable than the perpoint classifier’s results at optimal parameter 
settings. 

For the optimal parameter settings of the ECHO classifiers, the per- 
formance of the Nonsupervised ECHO processor is superior to the perfor- 
mance of the perpoint classifier only with respect to the CPU time required 
to perform the classification and the variability of the classification 
results. The Nonsupervised ECHO processor had training field performance 
which was worse than the perpoint classifier, for all the parameter 
settings which were tested (performances would have been the same if no 
homogeneous cells had been identified). These results indicate that 
there may be little or not advantage in classifying with the Nonsupervised 
ECHO processor with respect to classification performances. However, the 
classification variability results indicate that for the five channel 
simulated Thematic Mapper data, the homogeneity parameters selected for 
testing of the Nonsupervised ECHO caused more cells to be identified as 
homogeneous than were identifed by the Supervised processor as homogeneous. 
All homogeneity parameters for the Nonsupervised ECHO tests over this data 
set may have been too high. Simply lowering the homogeneity parameter may 
improve the performance of the Nonsupervised processor. In any case, it 
seems fair to conclude that the Supervised ECHO processor performs better 
than both the Nonsupervised ECHO and the perpoint classifiers. Its super- 
iority may be attributable to the use of class statistics in object iden- 
tification. 


9 
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c. Aircraft Results 

Figures la-21 through 2a-31 present the CPU time, the field 
center pixel performance, the full field performance, the training 
field performance, and the classification variability achieved by 
the Nonsupervised ECHO processor for the single aircraft run tested. 

RMS proportion estimate error values were not calculated. 

There are insufficient samples to perform an analysis of variance 
comparing these results to the single perpoint classification of the 
area. However, results of T-tests of the probability of each dependent 
variable’s perpoint result falling in the distribution of the nine 
measurements of that dependent variable for the Nonsupervised ECHO 
results is presented in Table 2a-17. It appears from these results that, 
for the aircraft data set, the Nonsupervised ECHO classifier: 

* requires less CPU time than the perpoint classifier, 

* has inferior field center pixel performance to the perpoint 
classifier, and 

* has superior full field performance when compared to the per- 
point classifier. 

These results are weak since so few samples are available. On the 
basis of Landsat and simulated Thematic Mapper results, it seems un- 
likely that the field center pixel performance is significantly superior 
for the same data set in which full field performance is significantly 
inferior. 

Table 2a-18 presents the analysis of variance results for the 
effects of the Nonsupervised ECHO homogeneity and annexation parameters 
on each of the dependent variables. There are insufficient data points 
for any effects to be statistically significant at a 10% confidence level. 
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TABLE 2a-17 


T-tests Comparing the Perpoint Results to the Distribution of 
Nonsupervised ECHO. Results for Six Variables 


Mean of the 
Nonsupervised 

CPU 

Time 

FCP 

Full 

Field 

Training 

Field 

Classification 

Variability 

Results 

Standard Devia- 
tion of Nonsuper- 

279.22 

70.53 

70.52 

85.22 

14.94 

vised Results 

39.43 

.77 

.72 

1.01 

5.65 

Perpoint Results 

434.00 

71.80 

69.50 

86.40 

22.29 

Perpoint result 

yes 

yes 

yes 

no 

no 


is different than 
the ECHO result at 
a 10% confidence 
level 


TABLE 2a-18 

Effects of Nonsupervised ECHO Homogeneity and 
Annexation Thresholds on Six Dependent Variables for 
One Aircraft Data Set 



CPU 

FCP 

Full 

Field 

Training 

Field 

. Classification 
Variability 

Homogeneity 

NS 

NS 

NS 

NS 

NS 

Annexation 

NS 

NS 

NS 

NS 

NS 
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Figure 2a-27 

Nonsupervised ECHO CPU Requirements 
for Aitcraft Data- Set 
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Figure 2a-28 

Nonsupervised ECHO. Field Center Pixel Performance 
for Aircraft Data Set 
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Figure 2a-2'9 

Nonsupervised ECHO Full Field Performance 
for Aircraft Data Set 
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Figure 2a- 30 

Nonsupervised ECHO Training Field Performance 
for Aircraft Data Set 
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Figure 2a-31 


Nonsupervised ECHO Classification Variability 
for aircraft data set 
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Nonsupervised ECHO Object Map Integrity Assessment 

The Object Map Integrity Assessment task called for the comparison 
of the objects found by the Nonsupervised ECHO with the known fields in 
the scene. Several problems arose which made it impossible to obtain 
quantitative results. The program written to construct object boundaries 
from the intermediate results tape took longer than anticipated to debug, 
and the delay strained the time and personnel resources allotted to this 
task. The display of the objects in a form suitable for comparison with 
the field maps of the LACIE segments also encountered difficulties which 
have not been completely resolved. As a result, no assessment of the 
object map integrity has been possible to date. 
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Figure 2a-32 


Nonsupervised ECHO Object Map 
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CONCLUSIONS 


Summary 

The Supervised ECHO processor (which utilizes class statistics 
for object identification) successfully exploits the redundancy of 
states characteristics of sampled imagery of ground scenes to achieve 
better classification accuracy, reduce the number of classifications 
required and reduce the variability of classification results. The 
Supervised ECHO processor requires cell size, cell-to-f ield annexation, 
and cell homogeneity parameters, input data, and a class-conditional 
marginal density statistics deck for both object identification and 
classification. The improvement in classification performance the 
Supervised ECHO classifier provides (over a perpoint classifier) ex- 
pands as the size of the objects expands in terms of numbers of pixels. 

The Nonsupervised ECHO processor (which identifies objects with- 
out the benefit of class statistics) successfully reduces the number 
of classifications required and the variability of the classification 
results. It is unsuccessful in improving classification performance, 
however. The Nonsupervised ECHO processor runs in two phases, a field 
extraction phase, then a classification phase. Cell size, cell-to-f ield 
annexation, and cell homogeneity parameters, along with input data are 
required by the field extraction phase of the Nonsupervised processor. 
The classification phase requires the intermediate results produced by 
the field extraction phase of the Nonsupervised ECHO processor and a 
class-conditioned marginal density statistics deck to classify the 
objects identified in the field extraction phase. 

Both ECHO processors provide information which may be of value 
in the training process. The Supervised processor produces a singular 
cell map which may be used to assess the adequacy of training for the 
area classified. Cells are identified as singular by the Supervised 
processor when the likelihood of the cell belonging to the most likely 
of the available classes falls below a threshold. Therefore, when a 
cell is categorized as "singular", it either contains pixels from more 
than one class or the spectral class of the pixels of the cell is not 
represented in the available class statistics. Groups of contiguous 
singular cells will occur when one or more spectral classes have been 
omitted. The singular cell map may indicate where additional training 
statistics should be collected. 

Using the "MAP" option causes the intermediate results tape pro- 
duced by the field extraction phase of the Nonsupervised ECHO processor 
to replace the response value in each channel of a data vector by the 
average response value for that channel of the field to which the data 
vector belongs. A false color photo of these data channels from such an 
intermediate results tape may help the analyst identify areas which are 
spectrally homogeneous and may reduce such noise effects as the six 
line banding which is encountered in some Landsat data sets. Figure 
2a-32 is an example object map produced by the field extraction phase 
of the Nonsupervised ECHO processor. Blanks on this map indicate 
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singular cells. Each pixel falling in an object identified by the 
Nonsupervised field extraction algorithm is represented by a symbol 
which has arbitrarily been assigned for that object. 

The greatest potential aid to the training process is the inter- 
mediate tape produced by the Nonsupervised ECHO field extraction routine. 
This tape contains the covariance and mean matrices of each field 
identified (without the benefit of class statistics) by the Nonsupervised 
ECHO processor. This spectral/spatial cluster may contain all the infor- 
mation an analyst needs to produce a more representative class statistics 
deck for a given area in less time than would be required to produce 
a statistics deck by conventional methods. 

Recommendations 


The ECHO processors are documented in this report, the Final 
Technical Report for May, 1977 [ 5 ], the ECHO User's Guide [6], and 
an ECHO case study [7] . All the ECHO algorithms are available to 
JSC (and other remote terminal sites) via the LARS remote terminal. 
It is recommended that : 

* JSC personnel should be encouraged (or assigned) to use the 
ECHO algorithms via the remote terminal. Only through this 
type of experience will NASA personnel develop experience 
and confidence in using this approach to classification. 

Only in this way will they gain sufficient insight into the 
characteristics of the ECHO classifiers to appreciate the 
potential impact of the ECHO approach in the context of large 
area surveys. 

* In order to further evaluate the utility of the Nonsupervised 
field extraction algorithm as a method for spectral/spatial 
clustering, tests should be initiated involving it and ERIM's 
"Glob and Blob" algorithm. These tests could be carried 

out either by LARS or by personnel at JSC via the LARS remote 
terminal . 

* The general utility of the ECHO processors and training pro- 
cedures and classification analyses should be investigated. 

For example, the singular cell map produced by the Supervised 
ECHO processor provides an indicator of the adequacy of the 
training statistics for a given area; groups of contiguous 
singular cells, where fields or other objects are known to 
exist, indicate the omission of one or more spectral classes 
from the statistics. Such investigation could be included 

in the training area selection/pixel labeling investigations. 
Also needed are training procedures geared specifically for 
the training of sample classifiers. Although these procedures 
have been evolving gradually, present methods are still more 
appropriate for training perpoint classifiers. 
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MAIN SUBROUTINES TREE OF THE FIELD EXTRACTION PHASE 
OF THE NONSUPERVISF-D ECHO PROCESSOR 


NS1SUP 


NS1INT 





NSECHO 

RDWRTE 


GATHER* 


VARCAL* 


BEGFLD 


*Delivered with Supervised ECHO Documentation in the Final Technical Report on 
NASA Contract NAS9-14970, May 1977 







MODULE IDENTIFICATION 
Module Name: NS IS UP 


Purpose: Supervisor for NS1ECH0 


System/Language : CMS /FORTRAN 
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Latest Revisor: 


LARS Program Abstract 


Function Name: NSlECHO 


Date: 8/20/77 
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MODULE ABSTRACT 

Supervisor for the field extraction phase (phase 1) of Nonsupervised 
ECHO (Extraction and Classification of Homogeneous Objects). 
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1 . Module Usage 
NS1SUP 

CALL NS 1 SUP 

NS1SUP is called by LARSMN to execute the field extraction phase of 
nonsupervised ECHO. There are no parameters involved. 


2. Internal Description 

NS1SUP receives control from LARSMN to perform Nonsupervised field 
extraction. NS1SUP calls NS1RDR to read and interpret the function 
control cards, then calls NS1INT to complete the initialization and 
compute array bases used to separate storage in the arrays for the sub- 
routine that performs the field extraction process (NSECHO). 


3. Input Description 

Not applicable 

4. Output Description 

Two messages are produced and written to unit TYPEWR (the console) . 
UNSUPERVISED ECHO FUNCTION (PHASE I) REQUESTED 
signifies beginning of the function. 

UNSUPERVISED ECHO FUNCTION (PHASE 1) COMPLETED 
signifies end of the function. 

5 . Supplemental Information 

See LARSYS System Manual for supervisor requirements. 
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FORTRAN IV G LEVEL 
FILE NSISUP 


0001 

0002 


0001 


=r 


0004 

0005 

0006 
000 / 
0008 

0009 

0010 


oou 


0012 

001) 

0014 

0015 

0016 


0017 

0018 


0019 


0020 


20.7 


NS 1 SUP 


NS1SUP 


OATS - 77228 


12011041 


PACE 0001 


LARS XXXX 
NR I T TEN 


By C. A. POMALAZA 


C SUPERVISOR ROUTINE FOR FIELD EXTRACTION PHASE OF NONSUPERV l SED ECHO 


subroutine nsisup 

IMPLICIT INTEGER • 4 CA-Z) 


COMMON /OLCCC'*/ BLANK, CAR0I20) • CHKOUT, COPFIL, CLASSR, CLASSX, 

1 CLUSTX, CCNPUT. CPVOUT, CRCROR, CROSEO, DATAPE, 

2 CU PL TP , OUPR'JN, CRR.MSO, FBP.NT, 

J FILESV, FLOliNO, HOATA, MCAOISB), 10(2001 , IMAGSX, 

4 (MARK, KLYfiO, MAP TAP HAXCHA, MAXCLS, 

5 P AGS I Z , PNCH, POINT, PRESUX. PRNTR, RE AO IN , 

6 RESTKT, RUNFIL, RUNTA8I 10,11 , 

■7 SUATA, SCPARX, SEPTPX, SPARE! 1C), TENPAS ( 30 ) • 

8 TPSTAK6I, UFLOX .TYPEWR, 

9 TOP, ARRAY! 125001 


REAL • 8 ARRAY 
REAL * 4 F RUCAL (5,101 

integer ♦ 4 co-entoai, oateisi, kcohiai, heo 2 ( 161 , timeisi 

INTEGER • 2 ULANK2 
LOGICAL * 4 CHKOUT 
LOGICAL • 1 BLANK l 

EUUI VALENCE IOATSAV, 10! II I . (CURRUN, ICOll , (FROCAL 1 1) , ID! 51 1 1 , 

1 (HFOl l).MEAO(el), (DATE! 11.HEA0I26I I . ( KED21 ll , HEAD! 191 

I !Mi?!Av|TPSTAT(lli: ‘FOMENT! 11.HEAIM7211, 

4 JSuPSCR.TPSTAT (2 1 ) , ( OUP I N, T<>SrAT ( 3) 1 , ( OASTAT, TPSTAT ( 4 1 

5 (CCPSER.TPSTAT (5) ) , (TRAUUT,TPSIAT(6II, 

6 (BLANK ,0L ANK2 .BLANK l | 


22 


C 

c 

c 

c 

c 

l 

i 

100 

l 

c 

4000 


COMMON /NS ICQ */ 

1 ANFLAG, ANN1. ANN2 # BUFPTS, BUFRCZ, 

2 Cl. C<-, CtLSlZ, CFLWTH. 

3 CSEU 3, 30) , CS6T3l3 f 30j, I NF 0(17). JP.'S, 

4 K9 f Ll» MAPFLC • MAXSIZ. HI.NSIZ* MTXSIZ* NDC200)* 

5 NEWRUN. NCClSt NCF6AT, NSR • NUMCLS. NVR f 

6 , ftQFlU, RQT APt * SELECU30), 

7 TT I TLE ( 30 ) • TOTPTSt VARSIZ, VECSIZ. 

8 CSEU30K o CSEl3l30)» FETVC3I30). FETVECC30), ZDUM, 

REAL*H array* 

REAL *4 CSET fC$ET3 v SELECl »ANN1,ANN2 
INTEGER • 2 CSEL.CSEL3, FETVEC* FETVC3 

logical • 4 mapflg 

EQUIVALENCE I VECS l Z .N0FET3) »lHTXSlZ f VARSZ3)t (CELUTH»GRSIZE I 


ECM00010 
ECM0032C 
ECMOOUIO 
•*6CH000*0 
ECHUC350 
ECH00060 
ECH00070 
• * E CHOCUd 0 
ECHC0090 
ECHOOIOO 
ECHOOUO 
ECH00120 
ECHOOI 30 
••ECHOOIOO 
ECH00I50 
ECH0ni60 
ECMOOl 70 
ECHOOI 8 0 
ECMOCl'JO 
ECH00200 
ECH00210 
ECH00220 
ECH00230 
ECHC0240 
ECH00250 
ECH00260 
♦ECH00270 
ECH00280 
ECH00290 
ECH00300 
ECH00310 
ECH00320 
ECH00330 
ECM00340 
1 1 E CHOO 35 0 
ECH00360 
ECM00370 
) .ECH00330 
ECH00390 
ECH00400 
ECH00410 


CALL READER AFTER INITIAL MESSAGE PRODUCED 
WRITE! TYPE WR , 1001 

FORMAT! • 1 0000 UNSUPERVISED ECHO FUNCTION {PHASED *t 

• •RECUESTEC (NSISUP)* ) 

CALL NSIWDR 

CALL PROCESSING ROUTINE • THEN PRINT TERMINATION MESSAGE 

CALL NSIINT 

MR I TE< TYPE WR. 9000) 

J CQ 0° UNSUPERVISEO ECHO FUNCTION (PHASED COMPLETED •, 

• M NSISUP)*) 


ECH005SC 
ECH00590 
ECH00600 
ECHC0610 
ECH0OS20 
ECH006 30 
ECi-00o*Q 
•ECHP06S0 
ECH00660 
ECH03b70 
ECH00680 
ECM00690 
CCM00700 
ECH00710 
•ECM00720 
•ECH00730 
•ECM00740 
CCH00750 
6CH00T60 
ECH00770 
ECH00780 
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Purpose : Reads function control cards 

System/Language : CMS /FORTRAN 

Author : C. A. Pomalaza Date : 8/21/77 

Latest Revlsor: Date: 


MODULE ABSTRACT 


NS1RDR reads and interprets all function control for the field extraction 
phase of Nonsupervised ECHO (Extraction and Classification of Homogeneous 
Objects). Checks are made for data validity. Also an intermediate results 
tape is readied for passing results to the classification phase of the 
Nonsupervised ECHO algorithm (phase 2). 


PURDUE UNIVERSITY 
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1 . Module Usage 
NS1RDR 

CALL NS1RDR 


This section lists the actions taken when the following control cards 
are read. 

INTERMEDIATE-TAPE The variable RQTAPE is set to the given 

tape number. 

INTERMEDIATE-FILE The variable RQFILE is set to the given 

file number. 


INTERMEDIATE- INITIALIZE 

INTERMEDIATE-NEWRUN 

CHANNELS 

ANNEXATION-MEANS 

ANNEXATION-VARIANCE 

CELL-HOMOGENEITY 

CELL-WIDTH 

PRINT-MAP 


The local flag INITFG is set to TRUE. 

This flag is used to trigger a call to 
TAPMNT to initialize the tape. 

A new run number is submitted to be used 
in the ID record of the intermediate 
tape. 

Subroutine channel is called to return 
VECSIZ, CSEL3, CSET3 and FETVC3 based 
upon interpretation of the channels card. 

Significance level 1 of unsupervised mode. 

Significance level 2 of unsupervised mode. 

Cell selection (Level 1) thresholds. 

Value of the width of a (square) cell. 

The flag MAPFLG in NS1C0M is set to 
TRUE. 


2 . Internal Description 

NSlRDR uses standard card reader logic in using LARSYS system sub- 
routines CTLWRD, CTLPRM, IVAL and FVAL in reading and interpreting the 
control cards. Subroutine CHANNEL is used to interpret the CHANNELS card. 

NSlRDR begins by calling TSTREQ to clear the stop/suspend flag. Then 
flags and arrays which will convey control card information are initialized. 
From this point the program, functions in a loop of reading and interpreting 
control cards until a DATA or END card is read indicating the card of func- 
tion control cards. After the control cards have been read, several checks 
are made on the data. A CHANNELS card must be supplied and the value of the 
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cell width is checked to be equal to or greater than 2. The intermediate 
tape is then mounted and positioned. Finally a list of selected options is 
printed. 


3. Input Description 

Function control cards for NS1ECH0 are read by LARSYS system routine 
CTLWRD . 


4. Output Description 

Control card error messages are written to both the printer (PRNTR) 
and the console (TYPEWR). A brief list of these follows: 

ERROR ON RESULTS CARD (TAPE OR FILE PARAMETER) - TYPE IN CORRECT CARD 

Syntax error in the TAPE or FILE specification. Standard corrective action 
is taken (requesting the user retype the card from the interactive terminal). 

ERROR ON CHANNELS CARD - TYPE CORRECT CARD 

Standard corrective action is requested. 

NO AUTO OR CHANNELS CARD - TYPE IN CHANNELS CARD 

Channels must be given via the CHANNELS card. After this message the key- 
board will unlock to accept a CHANNELS card. A response of carriage return 
will cause the keyboard to unlock again. 

CELL SIZE MUST BE GREATER THAN OR EQUAL TO 2 - DEFAULT OF 2 ASSUMED 

BOTH FILE AND INITIALIZATION OPTION REQUESTED FILE REQUESTED IGNORED. 
FUNCTION CONTINUES. 

The intermediate tape is initialized. Only file 1 can be initialized. 

In addition to this message the list of options selected is printed. 


5. Supplemental Information 

See LARSYS System Manual for card reader requirements. 
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6. Flowchart 


START 


INITIALIZE 

VARIABLES 



KEYBD 


READ 

CONTROL 

CARDS 

(CTLWRD) 


^ DATA ^ 
OR END 
vCARD> 



EXTRACT REST 
OF THE INFO 
FROM CARD 
CHECK VALIDITY 



CHECK VARIOUS 
)PTION COMBINATIONS 
FOR VALIDITY 


' ANY v 
ERROR 


ANY \ YES 
ERRORS? 


WRITE ERROR 

MESSAGE 

GET CORRECTION 




PRNTR 


MOUNT INTERMEDIATE 
RESULTS TAPE 
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•VAR I •/, 
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* S A ■ , * “PL C ' , * CLA«,«SSIF*,»ICAT* 
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r OR TRAN IV 
ILE NSlROa 


0*11 = // 22 *» 
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CECOCE INTERMEDIATE TAPE ESPECI F I CAT IONS 


IE I COL .CO. 72) 00 TO 101 

CALL C TLPRMtCAJtO. COL. INICOD, 6, COD 

IF ( CODE .EC. 41 CO TO 31U 

VECSZ > 1 

CALL IVALICAKO. COL. I. VECSZ, C 950) 

IFIVtCSZ .CO. 01 00 TO 350 

IF (CODE .EC. I) NEWRUN ■ I 

IF I CODE .EC. 2IRUTAPC « I 

IFICUOE .EC. 9)MQFILE * I 

CO TO 300 
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CO 10 900 

ERROR CM TAPE CARD 
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DATA CARO LNCOUN TIRED. . CHECK CCLWTH AND CHANNELS 


0102 

0101 

0104 

0105 

0106 

010 / 
0 1 Ob 

0109 

0110 
0111 
0112 

0113 

0114 

0115 


01 16 

0117 

0118 

0119 

0120 
0121 
0122 
0123 


IFIVfcCSIZ .GT. 01 Cl) TO 620 
ASS I ON 600 TO COVLC 
ERRNU" * 943 
ERRCCiH * 3 
CO TO 116 

1FIVECSI2 .IE. ISEL) CO TO 6»0 

IF! ISEL .GT. 01 CC TO 640 

ASSIGN 620 HI COVLC 

GO TO 645 

I * ISEL » 1 

DC 642 J * 1, VECSI 2 

SELECKJ) - SEIEC1 1 ISCl ) 

CONTINUE 
GO TC 630 

error cn Channel cr parameter cakc 

ERRNUM * 124 
E8RCUR * 3 
CO IT. 115 

IFICCIWTH . oE . 2)00 TO 650 

CELWTH * 2 

WKI TE I TYPE MR, 9620) 

MRITtlPANTR, 9620) 

FORMAT l* EXXX CELL SIZE MUST BE GREATER IMA) OR tOl'AL TO 2- 
•CEFAULT Of 2 ASSUMtO I.NECRDR) •! 


MOUNT INTERMEDIATE TAPE 


0124 
Cl 25 
0126 

0127 

0128 
0129 

0190 
0131 
01 32 
01 33 


MODE * l 

1F( .NCT. INI TFC) CO TO 680 
IFI R JF ILE .CO. 0) CO TO 660 
WH I TE I TYPE MR, 96501 
WRITEIPRNTr, 9650) 

FORMAT I' 1C050 BOTH FILE AND INI I IAL IZ AT ITN OPTION RECUES1EG. 
• * F I L C REQUESTED IGNORED. FUNCTION CCN1 INUCS INSIRCB)*) 

ROFILE * 0 

CALL T APMNT (ROT APC , ROFILE. MODE) 


CO 


700 


CALL TAPMNf IROT APE , ROFILE, MCDE I 


PRINT CUT PROCESSING PARAMETERS 


01 34 
01 15 


■R I TL I PRNTR .9/00) VECSIZ, CELuTh, NLwbijn.KUI API. HOT ILL, 

• ANN2.4NNI . ISELEC II J) . J - 1. VLCSIZ) 

FORMAT l//*OREAOER INFORM AT ION 0*// 

• ' NUMLCB OF CHANNELS **. I 9// 

• * CELL WIDTH *' .12// 

• * NE. RUN NUMBEfi * • . 19, • TAPE **,I5.‘ Flit 

• • MEAN AND VARIANCE THRESHOLDS A.«L *,FI0.5.* AN!)*,F 10.6/7 

• * CELL HOMOGENEITY TIIRESHOlCS ARE • , I T40. 5F 1 0. 5 1 I 
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MODULE IDENTIFICATION 


LARS Program Abstract 

Non-Supervised ECHO 

Module Name: TAPMNT Function Name: Functional Support 

Purpose : Mounts and positions Intermediate tapes 

System/Language ; CMS /FORTRAN 

Author : C. A. Pomalaza Date : 8/20/77 

Latest Revisor: Date: 


MODULE ABSTRACT 

TAPMNT mounts and positions the intermediate tape used in Phase 1 and 
Phase 2 of Nonsupervised ECHO. 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 


TAPMNT-2 


1 . Module Usage 
TAPMNT 

CALL TAPMNT (RQTAPE , RQFILE , MODE) 

Input Arguments 

1*4 Number of requested tape. A tape number of 0 Is 
a request for a scratch tape. 

1*4 Number of requested file. If RQFILE = 0 then 
the tape will be initialized by writing a record type 
1 (ID record) on the intermediate tape with filetype 
= 0 . 

1*4 Indicating usage of TAPMNT. MODE = 0 indicates 
that an intermediate results tape is being mounted 
with the ring out (only for reading) . 

Output Arguments 

RQFILE Current file position of tape. 


RQTAPE 

RQFILE 

MODE 


2. Internal Description 

See Output description. 


3. Input Description 

The record type of the intermediate results tape is read for each file 
up to and including the file needed. 


4. Output Description 

The following information messages are issued under the circumstances 
listed : 

10042 - is typed when a tape has been mounted and before TAPMNT posi- 

tions it. This message is not typed when the tape is being 
initialized or when the correct tape number was already mounted. 

10043 - is typed when MODE = 1 and the file tape has results in it 

(Check ^ 0) . 

10045 - is typed when the tape is correctly positioned. This is not 
typed when initializing a tape. 

After 10043 the user is asked whether he wished to overwrite the file, 
respecify a new intermediate results card or terminate the function. 
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10100 - is typed to allow entry of the new intermediate results card. 
This occurs when the user requests to respecify the inter- 
mediate results card. 

The following error messages are typed under the conditions listed: 

E361 - is written when the tape is being filed forward and a file is 

encountered with filetype other than zero (Check £ 0) before 
the requested file is reached and MODE = 0. 

E362 - is written when the circumstances for E361 occurs and MODE = 1. 

For message texts refer to the User’s Manual. 


5. Supplemental Information 
Input 

If a tape is mounted on the device and it is the incorrect tape number 
(as noted from the appropriate status words in CLODOM) TOPRV is called to 
unload the tape before the correct tape is mounted. If the correct tape 
is mounted TAPMNT assumes that the file number (as recorded in GLOCOM) is 
correct and moves the tape backwards or forwards to find the requested file. 

Output 

The tape is mounted with ring in for MODE = 1 and with ring out for 
MODE - 0. The tape is left positioned at the beginning of the requested 
file. When the tape is initialized LARSYS System routine TOPRW is used to 
do this. 
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FORTRAN IV 
FILE I4PMNT 


0004 

oocs 

0006 

ooo; 

0(100 

0009 

0010 


Oul I 
On l 2 

oo n 


Doll 

0018 


0019 

0020 
0021 
0022 
0021 


SUBROUTINE TAP**NI ( RQI APE , RQF IL E , PCD*. » 


C I .'.Prut EARS 0 *»x K" 1 0001 0 

c mm r 0002c 

C* ••.PMTOOO 

C ('•*100040 

C lAPHKT MOUNTS AND POSITIONS IN! t RMF. f) I A If RESULTS lAPfcb “MIOOOSO 

L WRITTEN BY C. A. PUMA LATA MMT 00060 

C MM T 000/0 

£••••••••••• ••••••••••••••••••••••■••••••••«••••«••••••••••••••* ••••••• *PMTOOOHO 

MM T 00090 

MMiooino 

MMfom i o 

■M 100120 
MM I 1)01 10 
MM I 00 14 0 
mm i no i so 
MMlODIf.n 
MM J OO I 10 
MM ion lit II 
MMIOOINO 
MM I 00200 
MM 1(1021(1 
«M I (10220 
MM I 002 (O 
MM T 0(1240 
MM t 0(12 Ml 
MMPO(l24n 
MM T DO,' 1 1) 
MM I 002 M (I 
MM I 002 Ml 
MM I 01) till* 
i MM I OO 110 
MMI00120 
MM I Oil I 10 
,MMt 00 (40 
MM1 00 I’.O 
MM1 00 (Ml 
mm i not/ 0 
MM I 00 (HO 
MM 1 00 100 


- CHKT1UT, ClIPFll, Cl AS SR i CLASS*. 

CL US T* • CONPUT. CPYOUT. CRORUR. CRDSFQ. T)ATAPf, 

L-UPLTP, DUPRU.N, rn.'MSC. T HPNt . 

F1LLSV, flOBNO, MDATA, HI ATM flfil , 101201)1, IMAUl*. 

IMAM*. AlYHD, MAPIAP. PA*CHA, MAXCIS, 

racsiz, pnch. point, prtsu*. prntr, rtaoin. 

(F. ST4T , HUNT I L , RUNTABI 1,1, II , 

SOATA. SL P AK * » si PIP* , SPARCIIOI. I i MPAS I (PI, 

IPMAllfal. llrLQi ,TYP!W4, 


A .(44V 

I HOCAL | S. (III 
• *. LIIMl Nl( )(, ) , OAII IS). 


DIMM, 111 021 16) 


1C M • 4 (.IIK'tU I 
I C A I • | III AH* | 

IIVAI l NC< I DA I SA V, (Dili I , I C UK 4 UN i I IM HI, 1 1 KUL Al II I , 1 1 
III 1>1 I I I .HI ilHHI ) . IIIAIEI I I. HI ADI ,’M I , {>1)02111.1 
I IMLI II .HE AIM SKI I , (CnMfMC I 1,10 ADI /2 I I . 


ML 1 II .HE AOI Ml I , (C 
PSAV. IPSIATII |1, 

I SlPSlR .IPSTAII2II , 10 

ic<ipsm. ipstaiis 1 1 , ii- 

1 Rl. ANA .01 AN*2.ni AN* I I 


!Un 


DA 1 A 1 UO/2 Cll*l)/ 

LA I A INTCIII)/* INTI */. ML SCUD/ *< 

■ u acih:/* • apt • .•( in •.•nl 


MIOl 


(Oil 


MM I 004 l 0 
MM I 004 2 0 

mm 1 004 in 

“M 1 00440 


? 004SO 


....APE 
ifc 
I NO 
UNI T 


' REPLY F/W - USER AS 10 WHETHER W CONTINUI WITH NUN 
IRKUR CUOF FROM CPFUNC I NOT USLO I 
U St D AS DUMMY RE AO I NO VLRSIOK NUMBER FROM TAPE 
INDICATES FUNCTION OF TAPMNI , SUPPLIED FROM CALLER 
* l MEANS WISH ro GEI TAPE K{ADY FOR WRITING RESULTS 
= O MFAlS WISH rn GFT TAPF READY FOR READING RESULTS 

■ RfcUUE STEP FILE NUMBER IN TAPE TU Of READIED. If RQFIlc 
IS INPUT AS A o. I MEN TAPMNT IS TO WRITE AN INITIAL 
RECORD ONTO THE TAPE. 

1 REQUESTED TAPE NUMBER ID HE MOUNTED. A ZERO MEANS SCRATCH 

■ FUC NUMBER AS He AO FHO" INTERME Cl ATE RESULTS TAPE 

■ TAPE NUMBER AS READ FROM TNT ERMF C I AT E RESULTS TAPE 

1 0S4-N GF TAPE BEING MOUNTED. THIS WILL BC HAPfAP OR CPYOUT 


FUND 

IF CORRECT TAPE IS MOUNTED GO TO POSITIONING 

IF IMAPSAV .tU. RCTAPEI GO Til 110 
CALL TCPAUIMAPTAPI 
GO TO 110 
UNIT « CPYOUT 
FJLNC * COPFIL 


MM I 004 Ml 
MM T 004 20 
MMt004<10 
MM T 004 90 
MM TOO SOO 
MM 1 004 | 0 
MMT00S20 
MMTOOS (0 
MMT00S40 
PMT00550 
T AP t MM T 0OS60 
MWTOOS/O 
MMTOOS60 
MMTOOS90 
MM T 00600 
••MMI006IO 
MMT00620 
MMTOOMO 
MM T 0064 0 
MMT006SO 
MMT00660 
MMTOOoTO 
MMT 00600 
WM| 00690 
MMT 00 TOO 
MMTOOMO 
MMT 00 T20 
MMT00/10 
MMTOOMO 
mm TOO ISO 
MMIOOT60 


PA Sf COO I 


FQR IRAN |V 
TILE IAPMNI 

0024 

U0?5 

0026 

002 / 


0029 
00 10 
0(1 U 

8 0 12 
0)1 
00 14 
00 IS 
00 16 
01)17 
00 18 

0(1 19 


0040 

0041 

0042 
Oo 4 i 
004 4 
004 S 
0046 


004 7 


00* a 
0049 

ooso 

00 S| 


00S2 

005) 


0GS4 

0055 

005 


0U5 8 
0049 
0060 
0061 
0062 
006) 


006 

0066 


OuM 


G LEvfl 20.7 


IAPMNIC DATF * //?! I 220)1004 


ii i cop Si h .10. 40i Am i gh rn in 


I MUM .10. O) C A | | MOUNT T ROT API , 


I It Nil ■ | 

RUT III l» 


WLANS I Ml f. Al | IS III mil MI/I 
01 Ml III 120 


n Mur ice ,ni . 

JCIIUNT • 8 10 
IDIItll • 4 0 1 A P I 
101)12) • I 

CALL TT.PUR I UN I I . JCOUNT.I RHIJH. IUHI 

RUT ILL » 1 

CALL UPRMTUMiri 

i.n in (no 

0 WRIIE llrp) W«, 91201 Vl/TAPC , -loritl 

O FORMAT <• 10042 POSITIONING I NTLHMLT) I Al I Hi SUI TS 1 

I * TO Ml^* j 14, ■ IIAPMNTIM 


API*. I*,, 

14. ■ I IAPMNI I * > 

CALL TCPRWTUNl 

t AU FIRST 4LC0RD AND SEl IF CORRECT T APT IS POITI II U 


JCOUNT « 000 

CALL T CPRO I UN I T . JCOUNT, ER4UR, ICO./, bll 

CALL TCPRWCUNIM 

IF (INC .EC. RQTAPEj GO TO I 10 

CALL CPFUNl 114.*M CP WRONG T API •, ERROR 1 

CALL TIJP.1UCUNI T I 

gc to no 

POSITION RFSULTS TAPC HI CORRECI T ILF 
IF IROFllE .01. 1) GO 10 140 


I )5 I 
I 

140 


ill TCPRWIUNITI 
0 TU 200 

« .(OFILE - FILNO 

J IS THE NUMBER CF FILES THE TAPC "USI 8E MOVED 

IF I J I 150. 200. 170 
I J * -J 

THE 00 loO 10 CP SACKS THE TAPE UP TC THE CORRECT FILE 

CC 160 I - l,J 
FILNO * FUND -I 

i call rcPHFiuiin 

GC IC 200 


> DO MO I » I, / 

JCOUNT . fl;o 

CALL TCPROIU-II T. JCUUNT . ERROR. ICO. 2. WIFI 
IF {CHECK .CV. 01 GO TO 210 
FILNO * FILNO » 1 
i CALL T CPFF lU'il T I 

IF FILNO AND ROT 111 ARI UNfUUAI , ll(T PROGRAM HAN GCClIltO. 
I IT IFIINU .Nl. 4(J*lll) CALI I H PR’, I || Oil, • SIHP 1 I 

NIX. 4EAU llir IfR*, | 41 CUHli Til llll Mil A*ll, (.HI I.n V Al I U| | v 


«M I<)||/-. 
MMMiONI 
MMIOllMI 

mm I non/ 


iMIlllIrt* 
mi (mu* 
M I Oil <| I 


MM I ()f)*(l 
MMIOO'I 
MM 1 llO’lj 
MM I 1)11*1 
PMI009< 
MM I 009 1 
«M I 0091 
MM 1009 
PM I 0091 
MMT 009* 
MMtOIOI 

mmidio; 

MMIOlOi 
MM TO 10 
MMTOIO- 


MMfO 
MMTO. 
MM T 0 1 


IT I MOTH .10. ni r.D Til ion 


PAGE 0002 
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F 0-<t«A j | V > 
I ILF IAP. 6 NI 


li*/:l 20./ 


OML * //?t| 


PAGE a out 


OU'.fl 
0069 
00 /O 
Oil/l 


oo m 
oo /u 

00 HO 
Oo ‘j | 
0002 
OUFl } 
O'JbA 
03<,5 
Oo R6 
O'JR 7 
O.itt 
Oo 1.1 

oiyio 
cum 
o:j ■./ 


OlCU 

0101 

0»i:2 

OlOt 

on* 

0105 


0106 

0107 

0106 

01C4 

Olio 

om 
01 12 
:n i 


m check . n t' . ni r.c ic ?io 
(.0 fn too 

/lu CALL E RPKNI I 161 *“i:cc. *00* I 
v«n ( i 
»i< I I. I lYPt *•*■ >/|0) | 

4/10 1 IW*A 111 OX « •/ J 1 1: • • I I) 

Cali r imai . 

7 10 NHlfi IIYP IRK, u 40 1 
kM I 1' (PM.'ITR, 9/40) 

4/10 F l.itMA I | ’ 1006 1 rill MAS HINDUS IN II - «l YOU NISH If. •, 

l 1 i;v..<n.<i m rut flit, rkspi.cifv t hi .Nina iusuiis card,*. 

/ • i:a riR“iN»it inf ji!it iiivirwpiif, ursPFciry, t r h>* imai i ) 

t , • i iapmni ) • | 

/IS CiJL * <1 

i .4, It iH > 1 

CA| L CIlNRI ICAHU. COl . HtSC’il). », C'lci, KIYRO, fRHC.lM 
IF 1 1 HP. COR .HI. 6 | CAtt R THAI FI 
•’ IIOO.MI.^O , COOF 


IN io 
■ 56 0 
. NSO 
1560 
5/ 0 
5B0 
S'lO 
6110 
610 
6/n 


4/50 » OR HA I I 

Cl:i * r 

i-RRC <R 

/*/ CALL CFL-Rr ICARO 
IF ILHRcOR rn *■ 

261 IF U.ni .E .. 

Call c il hr*' icaro, col 


v. ILI'R*' (vMU, COL, > 

. 1/50 ,/60, / /O) , COCE 

r. Fc.o'I jai; cun 

c 

*66 CALL Kt“Al . 

r. AifTHi .< slits Caro aftfk error 


r-|lFR l.'JURMCDIATC CARI- 

. COL, IN1C0D. 

I CALL RTMAIN 
•' *C 100 

CLACwO. 


, cocc. 


, Ci ICE , 


i r apnni i • i 
MVRO. F.RRCCR I 
1/651 


C 

C •'RC1CCSS PARAREIER ERRORS 

926BI CARD 
•GO* I 


/nr k*| IL CTYPlNR, 4/ 

4/6 H FORMAT | OK ,20 “6 I 

CALL EkPRNI I ERR Nil**, 


; liECUOfc TAPE SP..CIF ICATtON 
/Nil NECS.’ « t 

CALL I VAL 1 C ARI), COL, ROTAPC. VECSE, 
IF IVLCS2 .£0. 0) (,() TO 25? 


I II . SPEC I F I C A I ION 


. 26 4 

2o5 URWJM . 

CCL -• il 
IRRCUR * I 
UC I /6/ 

l INIFlALl/l SPECIF ICAF ION 

.*F4» iMiro * i 


VECS2. 125?) 


*00 kRI TL 1 1 YP 1 UR, 41001 

■Fiji’ F.VIkaI I • IOU65 RESULTS TAPI RUuMfO AND POSITIONER’, 
I • (TAPPNTI'I 


MMTO 
MM TO 
MM T 0 
MM I 0 
MM TO 

mm I n 
mm m 
MM T O 
MM TO 
MM 10 
r MM IP 
•MM TP 
MM I O 
MM 1 1* 

MM TO 
MM Ml' 

MM I 0 
MM | 0 
MM TO 
MM 1 0 
MM 10 
MMfO, 
MMTol 
MMT01 
MMT01 
MMT01 
MM 1 0 1 
mmTo | 
«MT01 
MM T 0 ] 

MM | 0 I 
MM 1 0 | 
MMTO 
MMtlT 
MMTO 
MMTO 
MMTO 
MMTO 
MM10 
MMTO 
MMTO 
MMTO 
MMTO 
MMTO 
MMTO 
MMfO, 
MMfOi 
mmto; 
mm r o i 
mmto; 
mmto; 
mmto; 
mmto; 
mmio; 
mmto; 
mmto; 
"M102110 
“MI02I20 
mmto? 1 10 
MMT0/I60 
MMI02I50 
MMIQ/160 
MMIO/I 10 
MMIO/UO 
MM10/IN0 
MMIO/200 
MMTO//10 
MMI07220 
MMT0Z/10 
“MI0//60 
MMTO 2/50 
MMT02/60 
mm T 07/ 1 0 
MM TO 2 2d 0 


b'iO 

TOO 

no 

//o 

no 

/60 


HOO 

VO 
HIO 
R6 0 
UNO 
H60 
1*10 
RM 0 
H'FO 


4/0 

<»to 

96 0 
450 
960 
410 
4U0 
14 0 
000 
010 
020 
010 

sss 

060 

010 


F:|.4TXA’4 IV 6 LlVI l 20./ 


T APMNIt 


HA If = / 12 1 1 2/011006 


PACE 0046 


FILE TAP- FT 


01 IB 
0| 14 
Ol/O 
0121 
0122 
.1123 
0|/6 
0 I 25 
012fc 
H 77 
Jl«H 


CALL ICP.4F ION I f » 

wo iriMcoF .Eh. oi cc t:i aoo 

MAP SA V * RJTA»F 
FILESV • F (l.‘4(l 
U (U4‘l 

6 CO CCPSLR * RCTAPE 

CCPF II = F 1LNC 
JtUHN 

sop call lrptn t "1 165, •stdp* i 

5 10 CALL EPPR’IT I 26 1 , • STOP* I 
F‘lC 


MM102240 
MM 107 100 
MMT07 110 
MMT0217O 
MM 102 HO 
MHI02 J60 
MM102350 
MMT02 160 
MMIO? 370 
MM102380 
MMT02190 
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MODULE IDENTIFICATION 


LARS Program Abstract 


Module Name : NS 1 INT Function Name: NS 1 ECHO 

Purpose : Initialization Including array base computation 

System/Language : CMS /FORTRAN 

Author : C. A. Pomalaza Date : 8/23/77 

Latest Revisor: Date: 


MODULE ABSTRACT 

NS1INT reads the area to be processed and computes the array bases for 
storage allocation of the variables used by the subroutine that performs 
the field extraction phase (NSECHO) . 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 
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1 . Module Usage 
NS1INT 

CALL NSlINT 

There are no arguments passed to NSlINT. It Is called by NS1SUP to 
allocate storage prior to the call of NSECHO. 

2. Internal Description 

NSlINT performs the following functions: 

1 - The data card of the area to be processed is read and interpreted 

using LARSYS system routine LAREAD. 

2 - The new ID record to be used in the intermediate results tape 

is built using part of the ID record from the multispectral image 
storage tape, (MIST) and information read by NS1RDR. 

3 - Check the validity of the channels requested comparing them with 

the ones present on the MIST. 

4 - Compute the space needed for arrays to be used. 

5 - Prints data summary. 

6 - Call NSECHO to perform the field extraction phase. 

7 - Call RDWRTE which transfer results written on disk by NSECHO to 

the intermediate results tape. 


3 . Input Description 

The field description card of the area to be processed is read by a 
call to LARSYS system routine LAREAD. 


4. Output Description 

Information concerning the parameters and the field selected are 
written on the printer (PRNTR) . 

Information and error messages are also printed. 

LAST SAMPLE NUMBER XXXXX OF FIELD XXXXXXXX EXCEEDS LAST SAMPLE ON TAPE. 
LAST SAMPLE RESET TO XXXXX. 

The last sample number is reset to the one present on the MIST tape. 
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YOU HAVE REQUESTED A CHANNEL NOT AVAILABLE IN THIS RUN. REQUEST 
CANCELLED. 

A channel number was requested which does not exist on the data run. 

FIELDS EXCEEDS LIMITS OF DATA. FIELD IGNORED. FIELD DESIGNATION 
FOLLOWS : 

The first column of the field is greater than the last column of data on 
the tape. 

A DATA DECK IS MISSING - FUNCTION TERMINATED 
The user has forgotten to include the field description card. 


Supplemental Information 


Not applicable. 
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. LLVI.l 70 ./ 


t>K I I T F N HYC. 


i ARRAY ALLOCATIONS FOH NCNSUPL HV I SI 0 I III b f * IHACI ION 


ECHOOOIO 

ECH00020 


001-4 
00 C 5 
OUvIb 
00C7 
00 . *8 
ooco 
0010 


> 4 I A- /I 

CO“MQN /GlCCOM/ PLANK. CARD! 20 I • CHKOUT. CORFU. CL AS SR ■ CLASS! ■ 
l CLUSTX . CONV»UT , CPYOUT * C«DRDR, CRUSEO, DATAPE, 

? (JIJPLTP, UUPRUN. CRRMSG. FBPNT. 

I IlLCSV. HDBND, HOAIA. HEADIB8I. 101700). |HAGEX, 

4 | NARK ■ KEYBD. “APTAP. -AXCHA, MAxCLS, 

5 PAGSI2. PNCH, POINT. PRESUX. PRNTR, R E AO I N , 

6 RESTRf, RUNFtL, RUNTABI 10 .)! , 

7 SO A r A , coiaci 


I , RUNFIL, RUNI AB( IU. II , 
iy,m a , St'PAHX, SEPIPX, SPARE) 101, Tt'MPASIIOI. 
IPSTAII6I. FtrtOX .IYPCWR, 

TOP. 4RRAYII2S00) 


8 1 AL • d ARRAY 
RLAL • 4 F .LOCAL IS. SOI 

INTEGER • 4 CO“ENT < 16) , OATEIS), HED 1(16). HE02I16I, TINE(5> 
INTEGER • 2 BLANK 2 
LOGICAL • 4 CMKOUT 
LOGICAL • 1 BLANK! 

COtllVALENC: I G A r S A V . I D I I ) ) r I CURH UN, I C( 3 1 » . I F ROC AL II ) , I Dl 5 III , 

1 (MED1I1) .MEA0I8) ) , IOATEI ll.HEAOl 26) I. (HED2I I) .HEADI )9) ) . 

2 I T (Mel 1 1 , ME ADI 58 I I . < COMLNTI | I.HEADI 72> I. 

J 1NAPSAV, TPSTAI 11) I . 

* ISCPSCR.TPSTAT 12) I • I DUPIN.TPSTAT I II), I DAS T A T, TPSTAT I 4 I ) , 

*» (CUPSER, TPSTAT (4»l , C TRAOUT, T PSTAT Ifcl » , 

b (BLANK, 8LANK2.PLANK1) 


CUP NUN /NSICJM/ 


OJ 1 2 
Og 1 3 
0014 
0-1 IS 
0 J 1 6 


I A-IFLAG. ANNI. ANN2. BUFPTS, BUFR02. 
i C) , C4, CLLSll . CELWTH, 

1 CSttlt.JO), LSET»tl,30t. INFOttT). JPTS. 
•. K9, LI, MlPFLG, M A x S 1 2 > M1NSIZ, MfXSIZ. 

6 .'IE-RUN, NCCLS. NOFEAT, NSR, NUMCLS, AVI), 


. ... 12. NDI 200 1 • 

•» NE-RUN. NCCLS 

t PIS. RESULT, ruaiui iwr.u, «yi»rc, mcimui. 

/ ITIILEIIO), TOTPTS. VARS 1 2 • VCCS 12 , 

bCSELlSOI, CSEUflOi. FETVC3I30I, FETVCCOO). ZOUM, 

9 A R R A Y 4 1 40C0 ) 

REAL*B A8QAY4 

Kf Al«4 CSFI.CSfclJ.SFLECl. ANN I , AN 62 
INTEGER * £ CSEL , CSEL3, FETVEC, FETVCI 
LOGICAL • 4 HAPFLO 

EOUIVAlENCi. (VECSIZ.NOFET 31,1 "TXSIZ.VARSZS I, (CELWTH, GRS I ZE I 


C Cl F INC LOCAL VARIABLES 


0021 

0022 

0023 


RLM»4 F L0C2IS, 30) 

t SSWSaifS E .|{?»f)U!: mIuI! 

I iFsitciIi'vFomi : lis'coulwo'?!! 


COL I NT , INFO (9 1 1 , 


FORTRAN I’ 
FILE NS1INT 


G LEVEL 20.T 


0077 

002B 

0029 

00)0 

0031 


00 36 
OG * 7 
OJ IS 
Od 39 
0040 
004 1 


ENO/'END •/ 


DATE » 77229 


(UNEB. 10(2211 . (COLA, 10(24)1 , I COl B. I C< 2S) I 


VARIABLES CE SC R IP T I ON 


ANNCP » 2 — GNSUPERVlSEO MUV G.L.R. ANNEXATION 
RECSKP ■ HCW MANY RECORDS TO DESIREO ARIA 
GCCK • NUMBER UF COLUMNS TO NEXT CELL DATA 
NOCCLL * NC OF CELLS 


ARTAY SIZES 


; stack 

. CLUSEO 
; uPt n 

unsiz 

EDA T A 
HOAIA 


I HALF aORO FOR EACH CELL HORIZONTALLY 
l VYTE FOR EACH CELL HORIZONTALLY 
- PY IE FOR EACH CELL HORIZONTALLY 


1 0 A T A BUFFER) 


CALL I SIRE Cl Nl)T L ) 


C READ IN ARIA TO PC CLASSIFIED 


INITIALIZE NEW 10 


00 too I * 1.200 

ADI I I * JCIll 
CCN 1 1 NUS 
SAV10 * NO (61 


SET UP CALIBRATION VECTOR 


m 


C SElOP * 2 — UNSU p ERV I SCO MUV CEIL SELECTION VARIANCE 


iLL LAREACI INFO. CODE. CRDkOH, 6. 0) 

ICCCE .N t . 0 1 CALL ERPRNT ( 104, • S TOP • , 

'IATAPE, ID. CRROR, R UNf AQ , (MARK) 

RUNERRI ERROR . INFOI | ) ) 


•ECHO 

:fSSg 

•ECHO 

ECHO 


CU (20 J*l, IQ 
OU 130 I«1 ,3 

IFICSET3M ,JJ .EU. -50000.1 CSEI3II.J) > F ROC AL ( I »2 , J I 
IFICSETJII.J) .NE. -50000.1 FRQC2 I ! '7.J1 ■ C SE T 3 I I * J I 
CONTINUE 
CONTINUE 


700 
210 
720 
7»0 
25 0 
25 0 
260 
270 
2S0 
290 
100 
310 
32 0 
330 
340 
ISO 
360 
170 
380 
190 
400 


CHECK AKEA COORDINATES 


004 7 
004) 
0044 
004 5 


LSAMNO * IUI6I- 6 

(FlIMFCm .GT. LS4MN0F GO TO 160 

IF I INF C 1 8 1 .IE. LSAN.NOI GO TC 140 

HR I TC I TYPE HR, 91401 I.NF0I8) . INT0I2), INFO 131. LSAMNO 
FORMAT I • IC04B LAST SAMPLE NUM3cK I • .IS, ' I OF FIELD', 

(.AST '/ax,' SAMPLE ON TAPE. LAST SAMPLE RESET 


EXCEEDS, LAST 
• 'TO '.IS,' INI SINT I 

1NF0I81 * LSAMNO 


' I 


ECHO 

echo: 

ECHOi 


1450 

:460 

470 

400 

490 

500 

510 

1520 


A-23 



r OR FRA 4 IV 

mi nsmnf 

00*8 


OUNtf 

Oil'll 

Oil'.? 


00 S* 
Oil'S 5 
00 •*6 


□010 
Oi.bl 
006? 
006 1 


006* 
0 j 6 S 
0066 
Oi'6 7 
UH 6 H 
0069 
00/0 


OHM 
0(1 / 1 
0 U>3 


0 )/* 
Ou/S 
00 76 
00 7/ 
00 7« 
00 79 
OCHO 
OOH! 
00 d 2 
OOu l 
OOB* 


OH MB 
0 UB 6 
008 7 
OOfld 
0089 
0D90 
00 >1 
0092 


ilvti ;*n./ nsmnfc Dot * //2/9 i?nsic>iM 

i*-i irii.iFbi*! .Gi. ii)<?on oo r* tbo * cum sic 

IMI.JFUtsl .CF. 10(2011 INFIX 51 * I0(2i>l LCH01S*0 

•*•••••«••«••*••••••••••••••••••• 

c •ECH0I560 

£ CHFC*. CMAN.LlS •ICM01S70 

c • ECMCXSBO 

(-.••••••••••••••••••••••••••••••••••••••••••••••••••••••••a* 

r, > JCFtFl f CHOI 600 

ID |’»t1 1 * | ,:i (CHOIMO 

ifi iitvcim .Li. iois)) u c ro iso icxntszc 

MJFlI I * .NCUFI -I I CHOI 6 10 

*9 1 r* i prnfn. -iisoi icnoibso 

uin nri’e.a, *»|S0l lOOlbso 

1 1 >11 FllRMAF I * 100*6 YOU HAVE REOUESILU A CliANKI l 901 AVA II AHl L •, ltnfl16f><l 

• 'in tiiis u'ip muurst csncrnm inisinii*i ichom/o 

SF:ll» 1C LCM0I6H0 

I -.1 LCNII-NIE ECHO 1 690 

co rn zou i Cun i /oo 

C«» •••••••••»•*•♦*•♦••*•••••••••••• •♦•••••••*•••••»• 710 

C *(0101720 

C ILLEGAL AREA REUUCSICU f 08 CL »SS 1 * I C A 7 lIJN *10101710 

c •I'Liini 7*0 

c«*«« •••••••••■•>••*•••••••••••••••••••••••••••«*•••»••••••••••• ••••••••ECHOI 7S0 

MO VHt II l OR'IM, 9 1 60) i NF X ECliOl 760 

kHllt 1 1YPE »H , »16')1 INFO l CI'Ol 7 7 0- 

•Uni rCRM.lt I'o* , |X. Ib.2x. 2A*,2X,6I S,9A*,A2 1 [CHOt/RO 

CALL ERPHNII 2*2, 'SlOP* I ECHO! 790 

€•••••••••••••••••*•••••••••••••••••••••••••••••••••••••••••••••• •••••••(□! ifM boo 

f. *1010 I Bio 

«. St F Ul» THE NUCbfeR OF PI) INFS AND UCCFINI INrr VLCKW *(0x01820 

C •ECMOIUJO 

(;»•••«>•••••••••••••••«••••••••••••••••••••*••••••••••••••••••••••••••• *ECH<) 1 8*0 

/co jpfs * (iNroid) - iNroi/i • nroiv i >/< i.nhii imchwthi icnoiaso 

prs * JPFS*CtLWlH ECMOH60 

IIIFCMI > INroi/i • FPTS-1 I • I TIP 01 91 (CHOI870 

UCA = I NF Cl 61 *CtLWFH (CHOIRHO 

INFUISI • | Ml 0 ( * I ♦ 'SICA* I I INrm-»l-INFMt* 1 *IM I1IMI /CLC« >- I NFIX* ) LCMO l !190 
ASH • PFS * 6 ECM0I900 

AYR * N$.«* /ECS1 2 ECliomo 



C *ECM0I9J0 

; MJDirv THE Ntr. 10 itCURO (01019*0 

C *t CHO I960 

£•••••••••••••••••••••••*••••••••••••*••••••••••••••••••••••••••♦•••••• *ECHO 1 960 

Nv.ll) * HO r Apt LCMO 1970 

NCI 21 * ftQrilC ECHO 1980 

NCI II * NEnMUN (0101/90 

c •••••••••••«••••••••••••••••••••••••••«••••••••••»> ••••••••••••••••••••ECH02000 

C *(0102010 

t COMPUM AH A AY SIZES IFOR ARRAY*) (CH02020 

C *ECM020 10 

r», ••••••••••*• ••••••••tCH020*0 

CUSIZ *w€ Lb IH*CCL 61 H ECM02060 

CLA r 1 » 1 tCH02i)60 

MHA7A I <10 AT A * I <0161 *VcCSI / ♦ i I /* ECH02070 

1AP0UF «ROAFA ♦ ;iVR*CELWFM CCHtlZObO 

CEL SUM » FAPbUF ♦ 2 ECH02090 

CClvCR » C 1 L SUM >VECSIZ*JPIS ECM02I00 

MiXSIZ • VI CSI Z* I VECS I 2 • 1 1/2 CCH021 1 0 

PlLVU • CiLCOR ♦ “FKSM*JPFS ECM02120 

PIKCJ* • P I AVAL * VECSM«CELSI2 CCH02110 

CATliSS * P I XCOR ♦ *FxSM*CELSM £01021*0 

IFIUAFI'AS-I .Gr. 8C00) CALI ERPRMTI 176, 'SFOP'I ECM02150 

£•••••**••*•*•••*••*••••••••••••••*•••*••••*•*••*••••••*•••♦••*••••*•*•*£0)102160 

C *tCH02l/0 

C ALLOC A FE cent run ARRAY USCM IN ANNCXAIION *£0102180 

C »tCMO? 190 

C 

C* • VCCSI 2 *CELS J 2 CCM 02210 

X 0 WSI 2 • J»IS*CELSI 2 ECH 02220 

$121 « I JP FS* 1 ) /* ECM 02210 

SI 2 / * IJPISWI /8 ECMO? 2*0 

SMI » IJPTSMI /2 ECH 022 SO 

$ 12 * * IVECSM ♦ 11/2 rcii 0?260 

SMS - I VECS M • J«M$ ♦ 1 W? It)! 02?70 

SI 26 * INF ISM ♦ I M 2 LCH 02280 


I Cl IRAN IV 0 LEVFL 20./ 


00-11 
009* 
0095 
009b 
009 7 

0098 

0099 

0100 
oi nj 
0102 
01 Vi 
010* 
0105 
0196 

0107 

0108 
01C9 
01 10 
01 II 
01 12 
01 I 1 
ON* 
1)1 IS 
0116 


0111 

0120 

0121 


Slid * INU^CL S* JP F S ♦ 11/2 
STACK • 1 

Closed « siack ♦ smi 

OMLN • CLOSCI. * SM2 
FLDSl l • OPEN ♦ S M2 
BOA F A I • FLOSM « SMI 
HDAfAl • OI.'AFAI » I PFS* 11/2 
PIXCXH « H . A T A I 

»»j<VAL * PIXCOH » (MTXSI/*CFLSM 
SIZE • C* 

FI.DVM. • P 1X9AL ♦ ISME*ll/2 
LI'llO * ILI.VAL 
FlOSUM • SMI • I • tlNMO 

nocu » «riLSuM • sir . 


M7 


CbLVAK . FLCCCR » 51 
FlOVAR « C *L VAR • Si.. 

CLLA/t . FLCVA1 ♦ SMS 
•UXSU- > CflAVF ♦ SM* 
auxva.i » .iuxsum ♦ srz* 

■OR A * AUXVAk * SI/* 

HUFF IN » Wullb 

8U«H» ./ * I ll)//<) - IlillFTfR -1M/JPT-. 

IflboFtU/ .L). 21 CALL LHPiCII I 1 76 , 'SINP'I 


MMFiLISSINC PAR AMI I 


CCM02*20 
L CM02* 10 


MX I IIIPRNFH, ll 70) r.UFP02,«AXiM 
r UR m* | | // / / • PROCESSING PARA-EFERS'// 

• * NUMBER OF HC*S Of CELLS IN f.CFFFR > 

• • MAXIMUM FlttC SME INO. OF PC INF S I 


PN 1 IF RUN HEADER 


CALL CE F (Mr I 1 I “E 1 
UKI Ft I PR IF 1.HC10) 
w.(| re I PKNF.t, FI ILL I 

CALL MEAOE <IVECSM.FFFVC1,CSf.l l.CSEF » . I0.FRJC4L I 


CLASSIFY THE 4ESICN 


012 * 

0125 

0126 
JL27 
0126 
01 29 
01 10 
•Jill 
.11 12 


call n<e:m:i 

• ARHLY ISFALKI , 4R R A V I Cl OS 1 0 1 . 4RR AY I C PEN I . 

• AR IAY I FLDSl/ I . A HR AY I 80 A fill, AM H 1Y ( RC 1 F A I ) . A*H * V * I P I /C'..1 I , 

• ARRAYS <p| XVAL ) . ARRAY8CCELSU“I. ARRAYrICiL COP I , 

• AHHAY IFLOSU") « ARRAY I • LOCOP I • AROAYI CFLVAPI • A9> L r I F ( I,/A-< I , 

• ARRAY ICCLAVL I . APR A Y ( AU* SU**! > AH P AY I ANY V A/ I . 

• ARKlY(aOA<t, A AH AY i HUFF t A I , 

• AQHlYILIN'.CI. AHRAYRI8DAIAI , A Rrf AYn I F A I , <0 i A Y H I I APlUJF ) I 

CALL F',F <C.I .OF* I 
:ua1 > 1 

HLAIi > I /LCSM* l-JIbl ♦ II/* • 1 
• C F p F S * I vECSI 2* 13I6F » 11/2 • 2 
•TRObPF * JhF5*CEL*FM ♦ 6 
MFPIS * JPIS** 

1-UFPrS » JYF'(** 

•PFS = JPFS • I 
CALL FShRI: I ARRAY I eO*I» J , 


ARRAY*! Y.P IS) .AP3AY*l*IPf$ I | 


AR? AYI-OAT A ) , NNU.PT. ARRAY*! II. 


FOR IRA'. |V 
FIVE NSMNt 


. LC'rftt 20./ 


O'F € • / 7229 


•ECHOIOSO 
•ECM7MO60 
ECHO 10 70 
ECMOIOSO 
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MODULE IDENTIFICATION 


LARS Program Abstract 


Module Name; NS ECHO Function Name; NS 1 ECHO 

Purpose ; Perform field extraction and disk filing for NS1ECHO 

System/Language ; CMS /FORTRAN 

Author : Date : 

Latest Revisor: C. A. Pomalaza Date: 8/21/77 


MODULE ABSTRACT 

NSECHO computes cell mean and variances row by row, performs cell splitting, 
annexation and field closures. As fields are annexed this information is 
sequentially stored on disk. As fields are closed, field statistics infor- 
mation is stored in a random access file on disk. 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 
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Module Usage 
NSECHO 


CALL NSECHO (STACK, CLOSED, OPEN, FLDSIZ, PRTBUF , RDATAT , PIXCOR, PIXVAL, 
CELSUM, CELCOR, FLDSUM, FLDCOR, CELVAR, FLDVAR, CELAVE, AUXSUM, AUXVAR, 

WORK, BUFFER, LINBUF, BDATA, RDATA, TAPBUF) 


All arguments are array bases computed by NSlINT. 


STACK 


CLOSED 

OPEN 

FLDSIZ 


1*4 Array dimensioned JPTS, where JPTS is the number 
of cell widths/line being processed. It is a revolving 
queue to store flags about the field status (which 
field storage areas are in use). 

L*1 Arrays each dimensioned JPTS. Logical flags for 
each field, used to ascertain closure status. 

1*4 Array dimensioned JPTS. Keeps record of number 
of pixels in each field. 


PRTBUF 1*2 Array used to store output symbols to print 

object map (dimensioned by core constraints). 

PIXVAL R*4 Array dimensioned VECSIZ * CELSIZ where VECSIZ 

is the number of channels used and CELSIZ is the 
number of pixels/cell used to store cell values 
during processing. 

PIXCOR R*4 Dimensioned VECSIZ * MTXSIZ where MTXSIZ - VECSIZ 

x (VECSIZ +1) . Used to store intermediate corre- 
2 

lation matrices. 


CELSUM R*4 Array dimensioned VECSIZ * JPTS and used to 

store cell means one row at a time. 

CELCOR R*4 Array dimensioned MTXSIZ * JPTS and used to store 

cell correlation matrices one row at a time. 

FLDSUM R*4 Array dimensioned VECSIZ * JPTS used to accumu- 

late field means. 


FLDCOR R*4 Array dimensioned MTXSIZ * JPTS used to accumu- 

late field correlation. 


CELVAR R*4 Arrays dimensioned JPTS * VECSIZ and used to 

FLDVAR store cell and field variances. 


CELAVE 


R*4 Array dimensioned VECSIZ used as intermediate 
cell mean buffer. 
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AUXCOR R*4 Arrays dimensioned JPTS * VECSIZ used as auxi- 

liary storage during variance processing for candi- 
date annexations. 

AUXSUM R*4 Array dimensioned JPTS * VECSIZ used to store 

means for candidate annexations. 

WORK R*4 Working area. What remains of ARRAY in GLOCOM 

after ZCOV. 

BUFFER 1*2 Hold field assignments from current and previous 

line of data. Dimensioned 2*JPTS*2. Buffer (l,x) 
contains the field number of cell x, BUFFER (2,x) 
contains the relative field pointer into FLDCOV 
(covariance of open fields). 

BDATA L*1 Array dimensioned VECSIZ * ID (6). Receives raw 

data values. 

RDATA R*4 Dimensioned VECSIZ * JPTS * CELSIZ holds calibrated 

data values from the tape. 

TAPBUF 1*4 Symbol storage area. 

LINBUF 1*4 Output storage area. 


2. Internal Description 


NSECHO operates in a loop over the lines requested. It retrieves a 
row of data, annexes cells into fields on a column basis. It then writes 
the buffer line to disk. Reads the next row (cell width lines) annexes 
it by columns and by rows to the first row. Fields that are closed are 
written with the field number to the intermediate tape and the buffer is 
written to disk. The slots in FLDCOV are then released and a new row is 
read. This processing continues until the lines are completed. NSECHO 
uses GATHER to compute cell variances and means, ANNEX to perform annexa- 
tion and BEGFLD to open new fields. 


3 . Input Description 

Data is read from tape via GADLIN. 


4. Output Description 

Information is written on disk using unformatted write. 


5. Supplemental Information 

See BUFFER and SCRATCH file description. 
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Fll < IRAN IV i 
111 NSIU3U 


son? 

Ou'l 3 


0J04 

0JC5 

OUOb 

000 / 

OoOH 

ooc« 

O-.MO 


*»:»f (,«*>' rc 'LilFGRN FIELD E«nuCNCN PIIA5L in NONSUPI Kv igeu fchg 


SU0RUUH IE MSCCnC I 
• r »C * . CICSED. CPt- 

<. v I * V 4 L ■ l L SUM . CtLLUK, 'Lllllil. HlllUK. rlT/VOK 

L CL L A Ve , 4 U * SUM > «U*V4R, H>IKK, BUFFER, l INPUT, HD A f A . 

L HDAIA, lAUiiUF) 

IMPLICIT l\ICCC-l*4 i A-/ » 

CC“MUN /GL'-COM/ BLANK, CARDl/OI. CMKlUl , CCPFIL, CLAS5R. CLASS* , 
I CLCSIX, CCNPUt, CPVQUt, CHDKfM, CHOSE C, DAI API, 
l UIPLlP, CUPRUM, CRRNSG. FP.RrU, 

( I I L L SV , FLDU.NU, MDA 14, HC ADCItrt I » I0I20UI, I MACE * , 

4 IMARK, fktVB'J, 9APIAP, MAX CMA , MAXCLS, 

j I'lTSU, PNCH, POINT, PRCStlX, PKNlH, 9 E AC I U . 

6 RLSTRI. iUNFU, ".U A I A R I 10 , 3 ) , 

7 SUH», SlPA-U, Sr. P IPX, SPAKE I 11)1 , If “PAS I to ) ■ 

H IPSlAIIb). HFlDX tlVPFUR, 

•J TUP. A-UAYI 1/400 I 

Real • h a may 

HI Al • 4 F H.JCAI 1 1, <01 

l.i<M,i« • C0»lNtt IM . OAIUSI. Him MM. MllWUM, «»*•!<•>» 

IN I IGF R « . PLANK/ 

LOGICAL • CHkUUI 
LOGICAL • I BLANK I 

LUll (VAIL ?1C LUAI SAV.ICtll I . I C1MH UN. t CL L I 
1 IHI III I I I .Hi AIM HI I . ( '• Al l | | 1 , ill Al: f 

/ I I |Mf II I .HI ADISNI ) , I COM* Nil I I, I'l 


ill, (I KOCAt L L I, l()(M 1 1 , 
l:( /hit, lit: 0/1 I l .HI AIM HI I I 
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I Coooil/o 

I COU0010 
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•EChOOOSO 
CCMIlUObO 
LCHOOO/O 
( CHOOOHO 
•LCHOOONO 
ECHOOLOO 
ICHOOUO 
ECHOOI/O 
LCHOOI 10 
FCH00140 
rcoooiso 

LCOOOINO 

Lcoom/o 

tCHOOMO 

CCHOOHO 

f cooo/oo 

CCHOO/IO 
ECM0022P 
fChoo/ in 
ICM0O/40 

rcooo/so 

LCHOO/bO 
I CHOU / tO 
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1 CHOI) 350 
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,i conn i/o 

I CHfll) mo 
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I CnOOMH) 
I CH0O4I0 
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LCHOONin 


0.SI2 
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0/14 
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1,14 


CtlMMUN /NSICCH/ 
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■I AKRAY4I4000I 
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I .IE SC J • / CSEl.CSELl.fEIVEC.FEIVCI 
LOGICAL • 4 “APfiO 

EGUIVALE.NC: IVECSIZ.NOFEHI.IMIXSI/.VAKSZJI. ICELWIm.CRSIZF J 


ouzo 

OUZl 

Oo22 


iSFINE LJCAl VARIABLES 

r.lLGE-4 • 4 11N3UFI1). IAP4UFII), SIATJZ). SCKICHIZI 
INIE-CR • / BLCCk ( 4 I > BOA I A I 1 I > 4L?I?) 

lcuivalenc: ialziii. alp » 
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LOGICAL • \ UOSEOL J Q IS* . OPCNtJPISI, SlOP 
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0061 
0062 

0063 

0064 

0065 

0066 
006 7 
0068 
Ot)69 

0070 
OJ 71 


LtVll /0.7 NScCHOC OAIE . 7 7/2 F* 1201100/ 

INIeGER • 2 GyrFtRLII.PaiBUFC JI.SIACKIJPISI. CLASS/. JPIS2. LChOOFJO 

6 SINFUL. V ER r I . mCRI l • PLUS. F X . ClOP-UF(5l)iM FCHOOMOO 

INIEGtR • 2 S VM I 36 J . CCMOOaiC 

INICGEKP4 FLOSIZ(JPtS) FCM00820 

N t A L • 4 tCHOOH 30 

3 F L USU“ I VE C S I / , JP IS ) » FtDC3R(M7xSIZ.JPTSI. Ep‘00B40 

5 CElvARIVECSlx.JPISI . FLOVAR IVECSli, JPlil. C(>00a50 

£ CEL4VL1VE ‘.SI L ) . taflRKUI. FCH00860 

e AUXSUH vEESIZI . AU*VAR | vECSI Z ) . KOAIA|r;vR, CLLOlH), (ChOORZO 

A CHIS/R. MAXllA. Al, r.LJFRM, fl&AI. R ) | 10 ) ICH008»0 

IwUlVAUNC’ IFSiLlN.lNF 0(4)1. ILSILIN, INrO I 5 1 I , 1 1 I N I M , INH) ( 6 I I LCMOOSTO 

CAT a SING /* '/ IChOONOO 

SAIA SIAI /•SIAI*.* ’ / , SCI I CM /•SCIH*,' •/ ECMOU910 

CAfA SINvUL . vFAl I .HORIZ.PLUS.E* / *0* .• I X* / CChOOV/0 

C a I A S tm/ • a •,*. *,'6 •,*, ',*C ',•/ '.'O fCmOO'MO 

I •• • , * F ' . *G ■,•( •,*« »,*♦ • , • | •,••• ICuno )4 n 

.* * I •.*( ' , ' K • , ■ I * ■ *1 •,•? ','t *,*0 *, F ChOONSO 

l *1 '.'P '.’I '.'C •.*'* ’.’R ',*7 ‘.‘S •/ LCMOOlfcO 

C f.CM0O'#/0 

£ EChOOnro 



6 FCnOIOJO 

: st i up vfcior i or pRCPf.R read fchoiolo 

C LCmOIO/O 

••ECHO 10 3 0 

C ICM01040 

G LCHOI050 

LINLS * ILSILIN-FSILINI/LININr ♦ | EChOIOM; 

C. ECH0IO7O 



C CCH010NC 

C PH I I HEADLL. FChOIIOO 

C ECmOIIIO 

u 2n 

C iCFtOUIO 

(.2 * 4CJ LCHOII40 

FChOJISO 


CSR-r? . 
cefine FILL 
ME W I ,0 30 
K9.| 

Si 


221 10000. 200. L.K9) 


J K 7-1 1 400 
6 WR I I 1 I 22* ALP) A 

*H1 IE I PRNrR,)l?5l I I NF01 II . 1-4,91 

91/5 FORMAT MOANEA PROCESS'D LINES ' 

l I 7 1 • ' CGIU“ IS 1 ,I4,*-* ,IS,* IHY «, 
C - 1.NFUI9I*CCI«IH 

•I * I Nf 0 1 9 I 

l • INF 01/1 

JUS* .EL-IH 


IFIN .Cl. I'W) I 1*|0 
ca’in K-i.jptsi 

•PftlttUF LKl ■> 1/(100*111 

PKtbUF (K'JPISI) • MOCtl. 100*1 ll/I 10*1 II 
PR I DDF IK*2*JPIS1 I - MODI 1,10*111/11 
30 I « I »N 

II - / *JPT S I • 1 
LI - 0 

ca 40 i-i. ii.jpisi 

J - 1 ♦ JP7S1- 1 _ „ 

40 .Al It! PRNI^, /040I I PR IBUF L 4 I , K - l , J I 
9040 FCH-AII ' lllll) 

IFIN .LE, 9J9» u3 ID 50 
I * l 

CO 45 JM . JPISI 
PR T 6UF IJ) > MODI I .101 
45 I « l *N 

-Rl IF I PRNt (.90401 IPRTBUF II), I- I, JPISI I 
50 UR I T£ I PRHl 1 .9 1 2l)l 
9120 FJRM A I I *0 * I 


LCHO 

ECHO! 

ECHO 

ECHO 

tCMO 

c ChO 

EChO 

cChO 

:ChO 

Ur,3 


190 
200 
210 
220 
230 
24 0 
25C 

.:gc 

.CCtZ/r 
l CHOI 260 
Z £601290 
CChCI »0C 
:c«oi mo 

C CHOI 320 
ECHO I 330 
tCHOI 340 
ECHO! 15G 
ECHO I 360 
EChOI 3 TO 
E CHO 1330 

echoi no 

ECHOI400 
ECHO 14 10 
ECHOI4ZO 
ECHOI 430 
ECH0I440 
FCH0I45C 
ECMOI460 
ECHOI 5 70 
ECM014B0 
CCHOI4NO 
FCH01500 
FCH0I5IC 
CCHOI520 
CChO I 5 3C 
CChO I 54 0 
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FORTRAN IV (, LEVEL 
FILE NSFCHO 


00 72 
0073 
OOM 

0075 

0076 

0077 
0018 

0079 

0080 
OOd 1 
0082 
0083 
008*. 
OOdS 
0086 

0087 

0088 

0089 

0090 

0091 

0092 
0091 

0094 

0095 

0096 

0097 

0098 

0099 

moo 

0101 

0102 

0103 

0104 

0103 

0106 

01C7 


C INI I t A L I2E 


JPTS2 « JPTS- 
NEX T* I 

bottom.jpts 

00 100 1=2,4 
100 8L0CR III * INFO! I *51 
DO 110 I ■ l.JPTS 
STACK! I ) 

PRTBUF 

CLOSE 0 

110 OPENII I * .FALSE. 

CLCR * L I N I NT*CELWTH 
Cl ■ CELSI l - 1 
CETCON 


H » FLOATICl**VEtSm 

(VECS l Z • ! 2*V£CSI2*3|-| 1/ 16.* IVECSI2* III 

ivecs iz-n*^v|csu«?>/6. 


kvii) • SELEcim • sELEcnn 
120 CONTINUE 

NOLINt » 0 

PLINI * FSTIIN 

NEXT * 1 

I5VN • 0 

PR f KM » 1 

JPt * JPtS - I 

OLDRW » 0 

NEWRW > I 

CUR L 1 N • F STL I N 

C THIS IS THE BEGINNING OF THE OLO SIM LOOP 
ISO J * NOLINE •CELWTN 
SOLtNE = HCLISB ♦ I 
I - CURLIN 
CO 1S5 J * l.CELWTH 
BLOCK 111 * 1 

CALL GAOLlNl BLOCK 11 I , CSEL 3 . CSET 3, 1 0, OAT APE , VECS I Z ,NSR ,80* T A 
L RDAIAlljJI, ROLL, ERROR) 

IF (ERROR. GT.O) CALL L INE RR I ERROR , I , £ 1SS, C 155, t 1 SS ) 

155 I - 1 ♦ INF0I6I 
160 H - I 

CO 201 CELNUH - l .JPTS 

CALL GMHERUOATAIH, 1) .P1XVAL >P|XCOR«CEtSliM I.CEINUMI , 

. .CELWTH.VECblZ.MTxSIZ.NSR.NV.U 


201 


ELWTH 


C ENTER "AIN LOOP — R F AD CELL SUNS ANC CORRELATIONS • 


0111 

oil) 

Oil* 


OUT 

0L18 

0119 

811 ? 

fill 


151 J ■ NOL I NE *CC l w th 

1FIM0D1N0L1NE.2S) .EO. 0) MR I TE I T YPEWR . 92351 J 
9235 FORMAT (16.* LINES CLASSIFIED. IECH083I 1 ) 
CROkBS * I • (ClORk-1 )*I**JPTS ) 

KROwPS • I ♦ INEkRk-l )*IJPTS**I 
FRO kd S * ( NE wR W* JP T S •* I 

C GET CNE ROW OF DATA. 

C 

£ COMPUTE CELL STATISTICS ANO SCREEN OUT BAO CELLS. 

C J * NRCkOS 


12 * l 

I : i 

is * 1 


FORTRAN IV G LEVEL 20.7 


FILE NSECMO 


0128 
0129 
01 30 


0132 

0133 
01 3* 
01 35 

0136 

0137 
01 38 
0139 
01*0 


01*3 
01** 
01*5 
01*6 
01*7 
01*8 
0| *9 
01 SO 

0151 

0152 

0153 
015* 
01 55 

0156 

0157 

0158 


0159 

0160 
0161 
0162 
0163 
016* 

0165 

0166 

0167 

0168 
0169 
01 70 


DO 200 1*1 ,JPTS ECH023I0 

SUFFER IJt * 0 ECM0232O 

ECH02130 

RETURN VARIANCES FROM SUM AND CROSS PRODUCT MTX ECh023*0 

6CM021S0 

170 CALL V ARCAL I CEL SUM ( t , I2I,CELAVE,CELC0R( l« 121 « ECH02360 

C CELVARIl, ISI.CEISIZ.VECSIZI ECH02 370 

176 CO 177 K» l • VECS 1 1 ECM0238O 

IFICEL VARIK, IS) .ST. R9 1 K I »C EL SUMI K , 12 1 *CEL AVC IK II CC TO 190 ECH02390 

177 CONTINUE ECH02*00 

GC TO 198 ECM02*10 

ECh02*20 

Ml CH02* 30 

E CM02**0 

INMOMOGENtCUS CELL. CLASSIFY IT SEPARATELY. ECH02*S0 

ECH02*60 
••••••••lCh02* 70 

190 CONTINUE ECHO2A90 

K * J ♦ 1 ECM02500 

00 196 l*J,K ECH0251O 

196 BUFFER IL I = -1 ECH02520 

PRTBUFIll = SING ECH02S*0 

198 12 * 12 ♦ l ECM025A0 

15 * 15*1 ECH02550 

200 J * 4 * * CCH02560 

IFIULORU .EO. 01 GO TO 265 ECH025 7C 

ECM02580 

)• ••••••• •••••*eCH02590 

ECH02600 

SCAN OL DRW FOR VCRTICAL BOUNDARIES ANO SINGULARITIES. ECH02610 

ECH02620 

i# •••••••••••** E CH026 30 

„ ECH02640 

K ■ CROWDS ECM02650 

GO TO 228 ECH02660 

ECH02670 

ICH02680 

ECMf)2&90 

COMPARE CELLS WITH FIELDS. ECH02700 

ECHO2T10 



ECH027J0 

COMPARE CURRENT ROW WITH PREVIOUS ROW. ECH027*0 

. „ ECHO 2 750 

228 J » ORCwBS ECH02760 

K ■ NR0N8S ECH02T7C 

CO 2*0 1*1 .JPTS ECH027dC 

H * BUFFER l J> ECH02790 

■ * BUFFER I J*l I ECH32R03 

IF1M . LE . JPTS .ANO. H .GE. 01 CL CSEO (HI * .FAL SE . ECH02d|C 

I F 1 BUF FER I K I .NE. 0) GO TO 239 ECHC282C 

IF I H .LI. 0) GO TO 230 ECM02BIO 

ASSIGN 229 TO RETURN ECHC28*C 

GO TO 1000 6CH02360 

229 | F I ANF LAG I 230,230,227 ECH0286C 

227 CPEN1HI « .TRUE. ECH02870 

BUFFER IK) = H ECH02800 

BUFFER 1K*1 »“M EC M 02890 

PRTBUFIll = QLD8UF 111 CCH0Z900 

GC TO 239 ECH029I0 

ECH02920 

COMPARE WITH FIELD TO LEFT, If ANY. ECH02730 

230 IF! I .60. || GO TO 239 CCH02950 

H * BUFFER IK-*) ECH02960 

P * SUFFER IK-31 ECH02T70 

IF 1H. LE.O) GO TO 239 ECH0298D 

ASSIGN 232 TO RETURN ECH02990 

GO TU 1000 ECH03000 

2*2 I F I ANF LAG I 2385,2385,238 ECHOIOIO 

235 CPENCHI « .TRUE. ECH01020 

238 BUFFER (Kl * H ECH030J0 

BUFFER !K*I I * w ECHO 30*0 

PRTBUF (II = PRIBUF ! l-l I ECHfUOSO 

GO TO 239 ECHO 1060 
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tcriran iv . 

me nsl c»u 


pace coos 


c cnisi rules 1 1 * 4 T arc in niDRw ikii niii in m wr*. 

C 

IASI ' JRIUh;. * ! Jp I N - l I • A 

CU 2i.« rt'l.jfis 

ll IClOSEOIn) .Ox. CM'NIHII Cl Ill 2'.ll 
BMt.l" * I • fliOlHIlI t()N. JPISI 
r>iu.mii:rr-;Mi * •« 


Mr v 


CiGsE : I II i*is. 


* t a 

A l S | 


a i s ukmuwhs.i ast,* 
ii i-»uff» k«wi .1 *j.mi r.oiu a i si 

C*»*n milt 

AlP'.lUlf IK ll»* 1 1 
»lf«RUfl E4(u*ll 

w<i it (?/■•! mihuf i tm.:i .Bum <io*i i.fl esi/ ihi. ifldsuhi ii , 
t n.vrcsm .rfioccRi ri ,hi , ri = i ,* rxsr/i 

HHI Ii I FK.il 4, | I KUF | Ffc I U) , HUM 1 9 I O* I I . I l US 171 H I 

I CLOSLDlM) - .IRUl. 

I CPE. Kill S .FALSE. 

KAPP*C«OWB4* I 
KUPP=N»3WDS 


*SCA.N ACROSS kCu. CJIjIN NEW f ItUJS. 


0146 

cm 

onr 

01 IV 

02 CC 
1)2)1 
0207 
02>:i 
02OA 

02!is 

0? n b 


1*1 3 
u*n - 
w7l? 


s fc * vlOW'lS ♦ IJATs-ll • A 
< * 3**TS 
Cl AM.J-'IS 

; * i 
I * I - I 

I * • - A 

II * li'jFfERlK) 

CL = nlFFEMt. *1) 

ir im .n. .ij on re ?do 
ir i.» ,nc. >i co re ?r« 

II s L I • I 

h - STACK I .[ill 

NLXI ■ 1 * -UCI *1 , JPIS) 

toff FR IK) * II 

BUFFER IK* l I * LI 

CO * l I 

I sv** * i*N.-:ei isy**, ibi 
wure-jMii = syp(isyp) 

CALL 9 tor L CICELSUHI 1 ,NI , F LOSU**l 1 . H , C El COR t I . N I , F L DCOR I 1,1 
l CFLVAMI,.). FLOVAR I I .Ml .FlDSIZ I III 
c IIII.iU .04. DUF F c 9 1 1 > .NE . D) CD 10 ?80 

assign m re -RtruRN 
Cl i j iooo 

? IFUNFLAyl 2dO.280.27S 
5 BUFF t Alt! * II 
PCFFS-R!l*l> * CJ 
PR 1 CUE 11) = PWUUflN) 


Oil 1 
0210 
0219 
01? 0 


tl It OLII PRINTER LINE 


hrtcel * r jptsi*i i/cf twrw 

CU 120 III « . CEL W TH 

L = CURLIN ♦ II 1 1 - 1 ) • I NFD 1 6 I 
170 *R|rt«PHifR.V2’»OI l,C IPRTOUFl II . J > I , C El WT h ) , l*),PRTCEL> 


FILE NSFCHO 

0221 

0222 

0221 

022A 

0225 

0226 
0227 


8151 

02 30 


L LA EL 20./ 


52' 


OATE * 77228 


fcr mat ii/> ix, nun 
CU 3 JO JK = l.JPTs 
MO CLOrtUFIJKI = PRIHUF(JK) 

I AO CCRlIN * CURL I N ♦ GLC* 

C THIS IS THE END CONDI licit 

1A2 hRlTM 10H BUFFER IM»I*NROwBS, FR0W8S) 

IF I CUR LIN .uT. LSILINl GO TO AOO 

C POl.t I TO NEXT R3w IN BUFFER 
CL DRW * NEw.U 

360 NLWRw * H03neWRW,BUF;i02) ♦ I 
C IMIS IS THE AC TURN OR END OF ROW PROCESSING CONDITION 
IF (NEMAN . NE . MTrw) GO TO ISO 


0231 
0212 
0233 
02 3A 
02 IS 
02 16 
02 1/ 
02 18 

02 IN 


02A6 
02 A 7 
0?A8 


k 1 i ♦ iprtrw-iimowsu 

l - K ♦ ROWSU-CELSU 
OJ 165 I *K ,L ,CEL S I 2 
H * t.tiFFCMM 

JPTSI GO TO 165 


if.S CUN II NUE 
- ASSIGN ISO FC RCTURN 
GC Tii aooo 

AOO CUNT I Mil 


C CLASSIFY All 0**lN MELOS. 


A|«> L * NRIIWHS • IJPIS - |l • A 
1)11 *20 I (NrOw'IS t L , A 
M - BUFFER! I I 
IF IH . L? . 01 GO Ii) *20 
AL2I2I -bUF F LR I 1*1) 

WRI III 22 'ALP I HUFF ER| I I .BUFFER 1 1*1 l,FL ESI 2 IHI , IFLOSU"! II , 
t *1 .VFCSI2 I , IFLOCCR! 1 1 .H > , I L = 1 . HI *S I Z I 
A I A RQW * KERR" 

FIELC * H 
■ A2P CONTINUE 


C FINISH PRINTING tJUf FE I. 


0250 
07 SI 
02 S2 


C ANNLXAF 1CN ROUTINE 


I F ( FLOS I Z I HI . CE . WAXSIZI GO FO 1006 

CALL A6NCXtCElSy l '|l,| I , FIOSUH 1 1 »H I , AUXSUH , CELCOR I 1,1), 
t FLOCCR ( I .HI .CEL VAR 1 1 , I ) .FLOVAR I I ,M| , AUXV AR, FL US 1 Z I H I I 
lOHb GU 1G TElWi, I27<»,7S2.2T2» 


TINT Oil .NOW CF CCllG. 


02 55 

0256 

0257 

0258 
02 SV 


A 01 10 11 = 1* « .*RI»W-| I.ROWSIZ 
12 = II » CEISIZ - CELWFH 
CC AU02 I 1-11. 12. CELWFH 
IS = 0 

I* * II* -0WS1Z-CELS1Z 
CC AOO 1 16-1 )» I A .CELSIZ 


FORTRAN |v 

FILE :;sccmo 

07 40 
0261 
02 62 


0?6S 

0266 

0267 


DATE * 7722H 


I 7 * 1 1 ♦ CELWTH - 1 
CC A001 18 =16,1 7 
IS * |S ♦ I 
GC T*» *002 

AOOl PRTBUf I Is) = bUFFERCiei - JPFS7 

A0C2 PUN.- = °L I ME ♦ L INI NT 

PRTR« * •ROCIPRIRW.puFROZ) ♦ l 

CC TN ISO 

E.\0 


ECHOA5NO 
ECH0A600 
ECM0A6I0 
ECH0A620 
ECM0A6J0 
ECM0A6A0 
ECH0A650 
ECHOS 66 0 
ECM0A670 
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LARS Program Abstract 

MODULE IDENTIFICATION 

Module Name : BEGFLD Function Name: NS 1 ECHO 

Purpose : Initialization of field statistics 

System/Language : CMS /FORTRAN 

Author : Date : 

Latest Revlsor : C. A. Pomalaza Date: 8/21/77 


MODULE ABSTRACT 

BEGFLD initializes field statistics from cell statistics for the field 
extraction phase of unsupervised ECHO. 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 
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BEGFLD-2 


1 . Module Usage 
BEGFLD 

CALL BEGFLD (CELSUM, FLDSUM, CELCOR, FLDCOR, CELVAR, FLDVAR, FLDSIZ) 


Input Arguments 

CELSUM R*4 Array dimensioned VECSIZ and used to store cell 

means . 

CELCOR R*4 Array dimensioned MTXSIZ and used to store cell 

correlation matrix. 

CELVAR R*4 Array dimensioned VECSIZ and used to store cell 

variances. 


Output Arguments 

FLDSUM R*4 Array dimensioned VECSIZ and used to store field 

means. 


FLDCOR R*4 Array dimensioned MTXSIZ and used to store field 

correlation matrix. 

FLDVAR R*4 Array dimensioned VECSIZ and used to store field 

variances. 


FLDSIZ 1*4 Variable to keep record of the number of pixels 

in each field. 


2 . Internal Description 

BEGFLD initializes FLDSUM, FLDCOR, FLDVAR and FLDSIZ to be equal to 
CELSUM, CELCOR, CELVAR, and CELSIZ. 


3. Input Description 

Not applicable. 

4. Output Description 

Not applicable. 

5. Supplemental Information 

Not applicable. 
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FORTRAN IV G LEVEL 
FILE BEGFLD 


20 . / 


BEGFLO 


DATE = 77228 


12012047 


PAGE 0001 


0001 

0002 

0003 


0004 

0005 

0006 

0007 

0008 
0009 


0010 
oou 
0012 
00 13 

0014 

0015 

0016 
0017 


C- 

C 


EC40001 0 
EC400020 
EC400030 
EC400040 
EC400050 
-EC400060 
EC400070 


BEGFLD MOCIFIEO BY C.A. PO M AL A2A 

INITIALIZES STATISTICS FOR A NFW FIELC 

SUBROUTINE BEGFLD (CEL SUM, F LDSUM , C ELCOR , FL DCOR , CEL VAR, FLOVAR, 

& FLDSIZ> 

IMPLICIT INTEGER * 4 (A-Z) 

COMMON /NS1C0M/ 

1 ANFLAC, AN.N1, ANN2 . BUFPTS, BUFRQZ , 

2 Cl . C4, CELSIZ, CELWTH, 

3 CSE T ( 3 f 30 ) • CSET 3 < 3 » 30 ) , INF0U71, JPTS. 

, 4 K9, LI, MAPFLG , MAXS1Z, MINSIZ, MTXSU, NC(200U 

5 5 NEWRUN, NOCLS, NOFEAT, NSR , NUMCLS, NVR, 

6 PTS, RESULT, ROWSIZ, RQFILE, RQTAPE, SELEC1 ( 30 ) , 

7 TT1TLEI30), TOTPTS, VARS1Z, VECSIZ, 

8 CSELI30), CSEL3I30), FETVC3I 30 1 , FETVECI 30) , ZDUM, 

9 ARRAY 4 1 4000 ) 

C 

REALMS ARRAY4 
C 

REAL * 4 CSET.CSET3.SELEC1, ANN1, ANN2 

C EC400240 

„ INTEGER * 2 CSEL , CSEL 3 , FET VEC , FE T VC3 EC400250 

C EC400260 

LOGICAL * 4 MAPFLG EC400270 

C EC400280 

, EQUIVALENCE ( VECS I Z , NOFET 3 ) , I MTXS I Z , V ARSZ 3 ) , ( CELWTH, GRS1 ZE > EC400290 

C EC400300 

REAL <■ 4 CELSUM (VECSIZ) ,FLDSUM< VE CS I Z ) , CELCOR ( MTXSI Z ) , EC400310 

C FLOCORIMTXSIZ) , CELVAR ( VECS I Z ) , FL DVAR I VECS I Z ) EC400320 

EC400330 

FLDSIZ = CELSIZ EC400340 

CO 1 1=1, VECSIZ EC400350 

FLOVAR ( I ) = CELVAR ( I ) EC400360 

1 FLOSUM(I) = CELSUMt I ) EC400370 

GO 2 1=1 • MTXS I Z EC400380 

2 FLOCORII) = CELCOR ( I ) EC400390 

RETURN EC400400 

END EC40041 0 



LARS Program Abstract 

MODULE IDENTIFICATION 

Module Name: ANNEX Function Name: NS 1 ECHO 

Purpose : Carry out annexation tests and processing 

System/Language : CMS /FORTRAN 

Author : Date : 

Latest Revlsor : C. A. Pomalaza Date: 8/21/77 


MODULE ABSTRACT 

ANNEX performs the multiple-univariate tests for annexation of cells to 
fields and updates field statistics when the annexation criteria are 
achieved . 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 




ANNEX- 2 


1 . Module Usage 
ANNEX 

CALL ANNEX (CELSUM, FLDSUM, AUXSUM, CELCOR, FLDCOR , CELVAR, FLDVAR, 
AUXVAR, FLDSIZ) 


Input Argument 



CELSUM 

R*4 Arrays dimensioned VECSIZ and used to 

store 

FLDSUM 

cell and field means respectively. 


CELCOR 

R*4 Arrays dimensioned MTXSIZ and used to 

store 

FLDCOR 

cell and field correlation matrices. 


CELVAR 

R*4 Arrays dimensioned VECSIZ and used to 

store 

FLDVAR 

cell and field variances. 


FLDSIZ 

Output Arguments 

1*4 Variable that records the number of pixels 
each field. 

FLDSUM 

R*4 See Above. 


AUXSUM 

R*4 Array dimensioned VECSIZ and used to 
means for candidate field annexations. 

store 

FLDCOR 

R*4 Array. See above. 


FLDVAR 

R*4 Array. See above. # 


AUXVAR 

R*4 Array dimensioned VECSIZ and used to 
variances for candidate annexations. 

store 


2 . Internal Descriptions 

ANNEX performs the multiple-univariate tests for equivalent mean vectors 
and for equivalent covariance matrices between a cell and a field. If both 
tests are successful the cell is annexed to the field and the field statis- 
tics are updated. ANNEX calls the function FDIST to find the decision 
threshold values of the F-distribution necessary to perform the tests. 


3 . Input Description 

Not applicable 


4. Output Description 
Not applicable. 
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Supplemental Information 


Not applicable. 


ANNEX -4 


6 . Flowchart 
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FORTRAN IV G LEVEL 20.7 ANNEX OATE • 77228 12012030 

FILE ANNEX 

C JECOOOl 0 

C JEC00020 

C ANNEX M0DIF1E0 BY C.A. POMALAZA JEC0Q030 

C JECOOOAO 

0001 SUBROUTINE ANNEX (CEL SUN f F LOS UH . AUXSUH.CEL COR tFLDCOR.CEL VAR* JECOOOSO 

t FLOVAR.AUXVAR.FLOSIZ) JEC00060 

0002 IMPLICIT 1NTEGCR * A (A-Z) JEC00070 

C JECOOOSO 

C JECQ0090 

0003 CONNON /NS ICON/ 

1 ANFLAG, ANNl, ANN2, BUFPTS, BUFROZ, 

2 Cl, CA, CCLSIZ. CELWTH. 

3 CSET13.30), CSET313.30), INF0I17I, JPTS, 

A *9, LI. MAPFLO, MAXS1Z, MINSIZ. MTXS1Z. N0I200I. 

5 5 NEMRUN, nccls. nofeat, nsr, numcls. nvr, 

6 PIS, RESULT, ROASIZ, ROFILE. ROTAPE, SELECII30). 

7 TTITLEIJ01. TOTPTS, VARSIZ, VECSIZ, 

8 CSELI30I, CSEL3I 301 , FETVCJ<30i, FETVECI 30) • 20UH, 

9 ARRAYA(AOOO) 

c 

OOOA REAL'S ARRAYA 

0005 REAL'A CSE7 ,C$E7 3.SELECI , ANN 1 . ANN2 

C JEC00260 

0006 INTEGER • 2 CSEL .CSEL 3 ,FET VEC ,FET VC3 JEC00270 

0007 LOGICAL • A PAPFLG JEC00290 

C JEC00300 

0008 EQUIVALENCE (VECSIZ ,NOFE T 31 * ( NTXS IZ.VARSZ3), 1C EL WTH ,GRS I ZE I JEC00310 

c _ JECOOtZO 

0009 REAL • A CELSUyTvECSIZ),FLOSUN(VECSIZI.AUXSUN(vIcSIZiT - ~ JEC003A0 

t CELCORINTXSIZI ,FL0C0R(NTXS1Z1. JEC003S0 

t CELVARIVECSlZl.FLCVARtVECSlZl* AUXVAR (VECSIZ). JEC0OI6O 

C A1.A2.AJ.AA, 40, A6, FOIST, 0£ TERM, ALOG , SORT , FLOA T JECD0370 

C-------- ------- "2 2 — — — — JECP03 8 0 

C UNSUPfcRVIStO MUV G~L~K~'rEST 'jECOOAOO 

c multiple-univariate test for equivalent mean VECTORS JECOOAIO 

0010 C5 » F LOS I L - 1 JECCIK20 

OUll CS - CELSIZ ♦ FLOSIZ JECO0A30 

0012 120 J » 0 JECOOAAO 

0013 A | = 1. ♦ FO 1ST 1 1 ,C1*C5» ANN2) /ICl'CSl JECO0A50 

00 l A cc 130 1=1, VECSIZ JEC00A60 

0015 J « J»l JECO0A70 

0016 AUXSUHIII = CELSUNIil 'FLOSUMUI JEC0DA80 

0017 AUXVARIII = CELCORIJI ♦ FLOCURIJj -AUXSUN ( I) • AUXSUM I 1) /C8 JEC00A90 

0018 IF ( AUXV4RI I ) .GT. I CELVAR ( I I * FL0VARI1 ) ) *A1 ) RETURN JEC00500 

0019 130 CONTINUE JfcCOOSlO 

C MULTIPLE-UNIVARIATE TEST FOR EQUIVALENT COVARIANCE MATRICES. JECO0S2O 

0020 100 IFIANNI ,LE. l.E-251 GO TO 135 JEC00530 

0021 Al = F LQA 1 (Cl ♦ C5I JEC005A0 

0U22 A2 = (l./Cl ♦ l./CS - 1./A1I/3. JEC0055 0 

0023 Al = 3./<A2*A2l JEC00560 

002A A2 • II. - A2 * 2./A3I/A3 JEC00570 

0025 AA > FOISTII.INTIAJ'.OD , ANN1 )/A3 JEC005S0 

0026 00 110 1=1, VECSIZ JEC00590 

0027 Al « C EL VAR I I 1/Cl JECOOGOO 

0028 IFIA3 ,LC. 0.1 RETURN JEC00610 

0029 A5 » FLNVARI I )/C5 JEC00620 

0030 IFIAS .LE. 0.) RETURN JEC00630 

0031 A3® A2.UI«ALOGUCELVAR(I>*FLOVAR(UI/A11-CI*ALOG(AJ»-C5*ALOG(A5))JEC006AO 

0032 I F I A 3 .GE. l.l RETURN JEC00650 

0031 AS » 1. - Al JEC00660 

001A IFIA1 .GT. A5«AAI RETURN JEC00670 

00 35 HO CONTINUE JEC00680 

0036 IJS CU IAO 1=1 .VECSIZ JEC00690 

0037 FL0SUMI1) • AUXSUM! I ) JEC00700 

0018 LAO FLOVARIll m AUXVARIII JEC00710 

OOAO 150 FLOC3Rll* l i H FLicORm ♦ CELCORtl) JEC00730 

°o°oi\ Stffo : { 8 jt&ggfcS 

0061 RETURN JEC00760 

0066 ENO JEC00770 
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MODULE IDENTIFICATION 


LARS Program Abstract 


Module Name: FDTST Function Name: NS 1 ECHO 

Purpose ; Computes decision threshold of a F distribution 

System/Language : CMS /FORTRAN 

Author : Date : 

Latest Revisor: C. A. Pomalaza Date: 8/21/77 


MODULE ABSTRACT 

FDIST returns the decision threshold of the F distribution for a given 
significance level and a given number of degrees of freedom. 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 






FDI.ST-2 


1 . Module Usage 
FDIST 

CALL FDIST (NUMBER, DENOM, SIGLEV) 

Arguments : 

NUMBER 1*4 degrees of freedom of the F distribution 

DENOM 

SIGLEV R*4 significance level 

FDIST is a R*4 function. 

2. Internal Description 

FDIST looks at a table to get the value corresponding to the signifi- 
cance level. The significane level may be one of the following values: 

.1, .05, .025, .01, .005, .001 

3. Input Description 

Not applicable 

4 . Output Description 

Not applicable. 

o 

5. Supplemental Information 

Not applicable. 
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6 . Flowchart 
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LARS Program Abstract 

MODULE IDENTIFICATION 

Module Name : RPWRTE Function Name: NS 1 ECHO 

Purpose : Writes intermediate tape 

System/Language : CMS /FORTRAN 

Author : Date : 

Latest Revisor: C.A. Pomalaza Date: 8/ 20/77 


MODULE ABSTRACT 

RDWRTE converts the disk files written by NSECHO into a intermediate tape 
file which is the input to NS1ECH0 (the classification phase). 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 
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1 . Module Usage 
RDWRTE 


CALL RDWRTE (BDATA, RDATA, NROWPT, BUFFER, MEAN, COVAN) 


All arguments but NROWPT are array bases computed by NSlINT. 

BDATA 1*2 Array dimensioned as VECSIZ * ID(6) and used 

as in GADLIN (for raw data). 

RDATA R*4 Array dimensioned as NROWPT * VECSIZ and used 

as in GADLIN (holding calibrated data values from 
MIST) . 


NROWPT 


1*4 Number of pixels per line + 6. 


BUFFER 1*2 Array dimensioned JPTS * 4 where JPTS is the 

number of cell widths/line. It is used to read 
, the flag buffer from disk. 


MEAN R*4 Array dimensioned as VECSIZ. Used to store 

cell means. 

COVAN R*4 Array dimensioned as VECSIZ * (VECSIZ + l)/2. 

Used to store cell covariance matrices while pro- 
cessing. 


2 . Internal Description 

RDWRTE performs the following functions: 

1. The new ID record, the processing parameters and the field descrip- 
tion card are written on the intermediate tape file using LARSYS 
System routine TOPWR. 

2. The statistics file (mean and covariance) fromeeach field are read 
from a disk file and written on tape using TOPWR. 

3. The program starts a loop getting a line of data via LARSYS system 
routine GADLIN. It then checks the MAP option and if it is on, 
the appropriate flag line is read from disk, and the data read by 
GADLIN is altered to reflect field means if the pixels were annexed 
to fields. The means are read from the disk stat file. If the 
MAP. option is off the data read by GADLIN is unaltered. 

4. The line when finished is output to tape using TOPWR. The flag 
buffer line is also written using TOPWR. 

5. The loop started in 5 continues until the area is finished. 

6. An additional file (only the ID record) is written in the inter- 
mediate tape. 
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3. Input Description 

Multispectral Data is read from the tape using GADLIN. 

4. Output Description 

Intermediate tape for input to NS2ECH0 (classification phase) is 
written. 


5 . Supplemental Information 

See intermediate tape file description. 
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BEGIN 
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6. Flowchart 


WRITE NEW 
ID TO IT / 


/ WRITE PARAMETERS 
I AND FIELD DESCRIPTION y 
ON INTERMEDIATE TAPE / 
(IT) / 


STATS 

SCRATCH 


' READ STATS FROMy 
DISK AND WRITE / 
THEM TO IT / 


READ LINE K OF 
DATA FROM MIST, 


MAP 

OPTION? 


BUFFER 

FILE 


/ READ FLAG 
'BUFFER FOR LINE^ 
K FROM DISK / 


IS\ 

PIXEL IN r 
SINGULAR 
FIELD? 


f READ FIELD 
STATS FROM/ 
DISK / 


W 


REPLACE PIXEL 
VALUE BY MEAN 
VALUE FOR THE 
FIELD 
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RETURN 
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ROwRTE 


DATE * tn?t 


12010012 
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FILE RDwRTE 


000 * 

000 * 

0006 

0001 

0008 

0009 

0010 


0012 

0013 

0014 
001* 
0016. 


001 7 
0018 
0010 
0020 
0021 
0022 
002 ) 

0024 


CDMN3N /GLCCfJM/ BLANK, CARDI/O), CMKQUT, Cl 
t CLUSTX, CCNPUt, CPYOUT, CODA OR. CROSEQ. 0 
L DUPLTP. OUPHUN, ERRMSG. FBPNT. 

t FILLS*, FLDHNO, KDATA, HE A D( 88 I . IOC 200), INA..EK. 

t i Kim i keybo, maptap, makcha, maxCls, 
t PAGS1Z, PNCM, POINT, PRCSUX, PRNfR. “EA01N, 

L HESTHI, RUNFIl. HUNT ABI 10 , 3 I , , „ , 

L S0AI4, SEPAHX. StPJPX, SPAREllOl, Il-PASIIO), 
t IPSTAH6I, T TF LDX , TyPEWR. 

L TOP. AHA AY 1 12500 I 


HFAL • 6 ARRAY 
HEAL • 4 F HOC AL (5,101 
INTEGER *4 COMFMTI 16) , 
INTCGLH ♦ 2 BLANK 2 


ClASSR, CLASSI. 


LOGICAL 


i chkout 


DATL'IM, IlFOl I 1 6 1 i HID2IIM. HMD*) 


LODI VALENC i IDA) LAV, 101)1 ) • ICURRLIN, IOC II ) .IIRUCAL (II, lOIM )l , 

L I HU) 1 (| ) ,HLAD(8))fiOATEI 1 I , HI ADI /«> I I ■ I ME C/I D.HCAOI IV)), 
t CTIMt ClI.HtADCBBI j.lCUHENtil I, HI ADI 7>)J, 
t (MAPSAV.TPSTATKII, 

t I SkPLCR . T PLIAT I 2 I | , (DUPIN, TPSIATI I ) I , I OAST AT , TP ST A I 1 4 I I , 

L ICOPSCR.TPSIAT CSI ) , I TAAOUI , TPST A T I 6) ) , 

C (BLANK, BL ANK2, BLANK 11 

COMMON /NS I COM/ 

1 ANFLAG, ANN1, ANN2 , BUFPTS, (IUFHCZ, 

2 Cl, C A , CiLSlZ, CELWTH, 

) CSEUI.IO). CSETm.301. iNFOlin. JPTS, 

4 *9, LI, HAPFLG, MAXSIZ, Ml NS 1/ , wTXSIZ. NOI200), 

* NEWRUN, NCCLS. NOTfcAT, NSR , NUMCLS. NVR, 

6 PTS, RESULT. ROWS 1 Z , AOFJLE, ROT APE , SELECIIIO), 

1 TTITLECIO), TOTPtS, VARSIZ, VFCSIZ, 

Q CSELOOI. CSEL IDO). FETVCH30I, rETVCCIJOI. ZOUM, 

1 ARRAY4I4000) 

R E A L *6 ARRAY* 

REAL *4 CSCT.CSETJ.SfctCCt . ANN1.ANK2 
INTEGER • 2 CSEL , CSEL l.FE T VEC ,FET VC) 

LOGICAL • 4 HAPFiG 

EQUIVALENCE I VECS I Z . NOFE T ) I , 1 nj <S IZ , V ARSZ 3 ) , C CELW TH.GHS I ZE ) 

INTEGER • 2 BUFFER (BUFPTS > . K1212I, 0. ft, NNICC400I, AR2C2S0) 
INTEGER • 2 BLOCK 141, BDAT Al TOTP T S I 
RE At *4 RD* IA( NRU.PT .VECS I Z ) ■ NEANI VECS 1 Z 1 , COV AN I Ml xS I Z ) 

REAL* 4 S, L.ARC 12S) 

INTEGER* 4 ARM 125) 

EQUIVALENCE 1 K 1 2 I I I , KI6), INNIDI II, NT) I 11) 

EQUIVALENCE CARCII. ARIIll.ARZIll) 

LOGICAL*! C ( 4 1 


ROW 000 | 0 
ROUO0D20 
ROWOCO 10 
HOW00040 
ROWOOOlO 
MOW 00060 

HOwoooro 
ROWOOOBO 
R0W00040 
ROWOOIOO 
RDW001 1 0 
HDW00I20 
KDWOOl 10 
ROW 0 0 14 0 
RD-OOISO 
HDWOO | 60 
KDWOOl 70 
H0WOOIHO 
RDWOUI *0 
KOwOOJUO 
RDW002 | 0 
HDW00220 
RDwOl)/ III 
HUwOO/41) 
HUW00240 
RUW00260 

ROuoa/tu 

Ruwno/eo 

kuwoo/90 
ROWOOIOO 
KUhOIMIO 
HUWOO I/O 
HOWOOUO 
HOW 00140 
ROwOO 1*0 
HUWOO 160 
HUWOO wo 
ROWOO IM 0 
ROWOO 190 
R0WOO40O 
KUW004I0 
HDW00420 
RDW004IO 


ROW00600 

ROW0O6IO 

RDW00620 

RDW00630 

R0W00640 

RDW00650 

ROW00660 

RDW00670 

ROW00680 

RDW00690 

rowoozoo 

RDWD0710 
KOW00720 
ROwOO? )0 
ROwOO 74 0 
ROWOO 750 
ROWOO 760 
ROwOO 770 
ROwOO 780 


FORTRAN IV 
FILE HOWRTE 
0025 


002? 

0020 

8 029 
010 
00)1 
00)2 
00 I) 
00)4 
00)5 
00)6 
00)7 
0018 
00)9 
0040 
004 1 
00*2 
004) 
0044 
0045 
00*6 
0047 
0048 
0049 
0050 
0051 
0052 
005) 


0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 
0062 
006) 

0064 

0065 

0066 
006 7 
0060 

0069 

0070 

0071 

0072 


00 77 
0076 
0079 
0000 
OUBI 
0082 
OOB) 

0084 

0085 

00 86 

0087 

0088 


G LEVEL 20.7 


CUUIVALENCLIC4.CI t)> 
H2*JPTS«CEIWTH*6 

SET U* INTE IMLOIATE TAPF 10 RECORD 


NO I * I ■ 1 1 


§i* 


NO 161*11 NFCC81-INF0I 71*11/1 NF 0(9) ♦ 
.NDI20)*IINFUI5I-1NF0I4I»|)/INF0(6I 
CO I 10 X>1 .VcCSIZ 
6N1U1 )9*2*K|*FETVC1IK) 
i 6NIDI40*2*K|«CSEL1IFETVC)CK)) 

COUNT. 800 
LRRUR-0 

CALL I0PWRIM6PTAP. COUNT .ERROR ,MD ) 

CALL TOPwR IMAPTAP, 60. ERROR. INFO) 

ARCII.CELWTH 

AR I 2 1 » ANN 1 

AH I 1 ) * ANN2 

L 2 * I NF 0 1 4 1 

L )• I NF C I 5 ) 

CO 120 J • 1,30 
ARI)*JI > SELECt I J I 
CONTINUE 


1C 


l* 


BLOCK t 2) * I NFO C 7 1 
BLOCK! ) ) » 1 NF 0 1 8 ) 

BLOCK I 4 I * 1 NFO ( 9 | 

RSZ *JP TS*4 

STATISTICS WRITTEN TO TAPE 
i HEAD AND WHITE OBJECT COV ANO MEAN MATRICES 
CO 1000 K I 9. 1 , L 1 

REAOIJZ’KDJD.B.F.IMEANI l ) , I * 1 , VE CSI 2 ) , ICOVANI J 1 , J. 1 , I 
CC0UNf*l2*VEC$U*MtxSm*4 
AR2I I ) *0 
AK2I2MB 
AK l 1 2 ) *F 
00 451 *0*1. 

ARI2*KC)*MEAi 

DO 452 KO* 1 »HT XS I Z 
452 ARC 2»VECSI Z»KOI«COVANIKO) 

1000 CALL TOPwR IMAPTAP, OCOUNT, ERROR, AR ) 

KLIN » 0 

L*. CELWTH* INF0C6) 

00 100 I A*L2, L ) . L 4 
REAOC 30) (SUFFER! 1 1,1*1 ,HSZI 
DC 1O0 JA ■ 1, CELWTH 


IA • 


I 


,ROU .ERROR ) 

K 3 ■ 0 
INN. 2 

02 CO 122 A2M.RSZ.4 
K3»K)* 1 

IF CBUFFERCA2I .EQ.-I.OR. .NOT.maPFLG) GC TO lift 
K16* BUFFER I A2 * 1 1 

READI?2'K16»0,B,F.(MEANC 1 I . I * 1 . VE CSI Z I , CCOVAN I J) , J> l, MT XS I Z ) 
K 1 6 - l 

DC 115 K2«l. CELWTH 
CU 115 K>1 , VECS 1 2 
MN.FETVCJIKI 
CALL RLCP ( "N, S ,L ) 

1 15 ROATAI K)*C EL wTH-CELwTH*K2 , K ) • CMEAMKl/F-ll/S 
GO TO 122 

110 DO 121 K2«l, CELWTH 
DO 121 K* 1 .VECS I Z 


ROwOO 790 
ROWOOHOO 
RDWOOBIO 
RDW00820 
RDWOOBIO 
RDW00840 
HDW00850 
R0W00R60 
RDWOOB 7 0 

Rowooaao 

ROwOO a 9 0 
ROWOOJOO 
ROW 009 1 0 
RDW00920 
RDW009)0 
ROW00940 
RDW00950 
RDW00960 
RDW00970 
KDWOO‘78 0 
RDW009VO 
RDWOIOOO 
RDWOIOIO 
ROwO 1020 
RDWOIOIO 
ROW01O40 
RDW01050 
ROW 01060 
ROWOIO/O 
ROwOlOHO 
RDWOI090 
ROwOllOO 
RDWOllIO 
HDW01 1 20 
RDWOllIO 
RDW01140 
RDWOI 150 
RDwOl 160 
RDWOI l 70 
RDWOI 180 
RDWOI 190 
HDW0I20O 
RD-01210 
RDW0I220 
RDW0I230 
RDW0I240 
ROW01250 
R0W0I260 
RDW01270 
RDWO j 280 
ROWO 1 290 
RDWOI 100 
ROWODIO 
ROW01320 
RDWOI 1 30 
RDWOI 140 
, OOAI A ,RDAT A I 1R0W0II50 
RDWOI 160 
RDWOI 170 
RDwOl 380 
ROWO l 390 
RDU01400 
RDWOI* 1 0 
ROwO 1*20 
RDWO 1430 
ROwO 1440 
RDW01450 
RDW01460 
RDwOl* 70 
HOW014BO 
RDW01490 
RDWO 1500 
ROWO 151 0 
RDW01520 
RDW01530 
ROW01540 


FORTRAN IV G LEVEL 20.7 


PAGE 0001 


FILE RDWRTE 

0089 

0090 

0091 

0092 

009) 

0094 

0095 

0096 

0097 

0098 

0099 

0100 
0101 
0102 

010 ) 

010 * 

0105 

OC Ob 

0107 

0108 

0109 

0110 
0111 


0112 
Oil) 
01 14 

0115 

0116 
0117 

8 1 18 
1 19 
0120 
0121 
0122 
0123 
0124 


Sill 


"N.FEIVCJCK1 
CALL RLCPCMN.S.L) 

121 RDA T A I K3*CELWTM-CEL WTH 

122 CONIINUE 
INN. 3 
C4-A 

KLIN « KLIN ♦ 1 
BOATA I 1 1 »KL IN 
BDAT A I 2 1 *RCLL 
DO MO K7* 1 , VECSI Z 
DO 1)0 K8M.RZ 

DO BOATA C 4*1 K 7-1 1 *RZ*K8 I *R0*TA|k8,K7 I 
INN.* 

C0UNT.VECSIZ*K8*4 
C0UNT2* VEC 5 1 Z *K8 
COUNT* * 4 

CALL MOVBTT I BOA I A, 0, 1 , ROA T A , 0, 1 , COUNT* I 

call MOvarr(8O4r*.<),j,K0Ar4,4, i.ckmtzi 

INM.5 

110 CALL TOPwHIMAPTAP, COUNT. ERROR, ROAIA! 


100 CONTINUE 

EMC INTERMEDIATE TAPE 

'CALL TOPEFINAPTAP) 

JCOUNT * 800 
DO 101 JJ * ) , 200 
NC IJJ) > 0 
101 CONI I NUE 

NDI 2 I > NDC2I * 1 
FlLESV * FlLESV ♦ 1 

Eitt ic?"Fl!!p'SEi JCOU " T,E “ lloa, '‘ D ’ 

CALL T CPRF ( NAPTAP ) 

CALL T OPE F IMAPTAP) 

CALL TCPBF IMAPTAP) 

CALL TCPBF IMAPTAP) 

RETURN 
END 


>K2,K>.IR0ATAIK)*CELWTH-CELWTH*K2,K|-L)/S 


550 
560 
570 
580 
590 
600 
610 
620 
6)0 
6*0 
650 
660 
6 70 
680 
690 
700 


Ml 


RDwo: 

rowo: 

ROWO 

rowo: 

rowo: 
ROWO 
ROWO 
RDWO 
ROwO 
RDWO 
ROWO 
RD.O 
ROWO. 
rowo: 
ROWO 
RDWO 
HOMO 
RCwO 
rowo: 

HOWO 
ROWOI 
RDwo 
ROWO 
MDW0 
RDwO 

RDwo 
ROWO 
ROWO 
ROWO 

rowo: 
rowo: 
rowo: 
rowo: 
rdwo: 

HOWO ! 
RDWOI 
RDWO! 

rdwo: 

aowoi . 

ROWO I 9*0 
ROWO 1 950 


750 

76C 

770 

780 

790 

800 

810 

820 

8)0 

8*0 

850 

860 

870 

880 

890 

900 

910 

920 

9)0 


















LARS Program Abstract 

MODULE IDENTIFICATION 

Module Name: NS 2 SUP Function Name: NS2ECHO 

Purpose : Supervisor for NS2ECH0 

System/Language : CMC /FORTRAN 

Author : C. A. Pomalaza Date : 8/21/77 

Latest Revlsor: Date: 


MODULE ABSTRACT 

N52SUP receives control from LARSMN. This supervisor performs no computation, 
but instead makes call to the subroutines which really make up the processor. 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 
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1 . Module Usage 
NS2SUP 


There are no parameters to NS2SUP. It is called by LARSMN when 
NS2ECH0 is requested. Control returns to LARSMN when the function is 
completed. 

2. Internal Description ; 

NS2SUP first calls the subroutine NS2RDR to read in all the control 
cards, and then calls NS2INT to read the statistics data and compute the 
array bases. NS2INT calls NS2ECH0 which performs the classification. 


3. Input Description 

Not applicable. 


4. Output Description 

Two messages are produced and written to unit TYPEWR (the console) . 
UNSUPERVISED ECHO FUNCTION (PHASE 2) REQUESTED 
Signifies beginning of the function. 

UNSUPERVISED ECHO FUNCTION (PHASE 2) COMPLETED 
Signifies end of function. 

5 . Supplemental Information 

Not applicable. 
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Flowchart 



* 

CALL 

HS2RDR 






•r 

i_n 


FORTRAN IV G LEVEL 
FILE NS2SUP 


20. 7 


NS2SUP 


DATE « 77229 


13051040 


0001 

0002 

0003 


0004 

0005 

0006 

0007 

0008 

0009 

0010 


0017 

0018 


0019 


0020 

0021 

0022 


0023 

0024 


C *NS I 0001 0 
C NS2SUP LARS XXXX *NS100020 
C WRITTEN BY C.A. POMALAZA NS100030 
C *NS 1 00040 
C *************** ******************* ******** **** ♦*****♦*<<****** ***»*****#»j si 00050 


SUBROUTINE NS2SUP 
IMPLICIT INTEGER*4 (A-Z) 

COMMON /GLCCOM/ BLANK, CARD! 20 ) , CHKOUT, 


COPFIL, CLASSR, CLASSX, 


CRDSEQ , DATAPE, 


IMAGEX, 


INTEGER * _ 
LOGICAL * 4 
LOGICAL * 1 
EQUIVALENCE 


*NS 1 0006 0 
NS l 000 70 
♦NS1000G0 
NS 1 00090 
NS100100 
NSI001I0 
NS 1 00 1 20 
NS100130 
NS 1 00 1 40 
NS100150 
NS100160 
NS 1 00 1 70 
NS100180 
NS1 00190 
NS 100200 
NS 1 002 1 0 
NS 1 0022 0 
NS 1 002 3 0 
NS100240 
NS 1002 50 
NS l 00260 


CLUSIX, CCNPUT, CPYOUT, CRDRDR, 

DUPLTP, DUPRUN, ERRMSG, FBPNT, 

FILESV, FLDBNO, HDATA, HE AD ( 8 8 ) , I D t 200 » 

IMAKK, KEYSD, MAPTAP, MAXCHA, MAXCLS, 

PAGSIZ, PNCH, POINT, PRESUX, PRNTR, READIN, 

RE S TR T , RUNFIL, RUN T AB ( 10 , 3 ) , 

SOATA, SEPARX, SEPTPX, SPAREI10J, TEMPASI30), 

TP STAT ( 6 ) , TTFLDX, TYPEHR, 

TOP, A RR AY I 1 2 500 ) 

REAL * 8 ARRAY 
REAL * 4 FRQCALI5.30) 

INTEGER * 4 C0MENTIL6), DAT E I 5 ) , HEDK16J, HED2I16), TIMEI5) 

CLANK2.NN 1DI400) 

CHKCUT 
BLANK I 

(DATS AV, 1011)1, (CURRIN.iDI 3) ), ( F RQC AL I 1 ) , I D I 51 

1 (HEDl(l),HEAUi'3)>, (DATE! 1), HEAD! 26) ) , I HED2 1 l ) ,HE AD ( 39 ) ) , NS 1 00270 

2 (TIME! I ) ,HEAD(58) ) , ( COMENT ( 1) , HE AO I 72 ) ) , NS100280 

3 (MAPSAV.TPSTATI 1) ) , NS100290 

4 (SEPSCR,TPSTAT (2) ) • ( D'JP I N , TP ST AT ( 3 ) ) , I DASTAT , TPSTAT { 4) ) , NS 1 00300 

5 (CGPSER.TPSTATIS) ) , I TRAQUT,TPSTAT(6) 1 , NS100310 

6 (BLANK, BLANK2, CLANK1) NS100320 

♦NS100330 

), MTXSIZ, NQFET3, NOPOOL, OFILE. OTAPE, 

FILE, RQTAPE.STKPTR, VECSIZ, WRKSIZ, 


, PQLSTK 


C 

C *NS l 00400 

C ***************** ********************************************** ********NS 1004 10 
C *NS 100420 

C CALL READER AFTER INITIAL MESSAGE PRODUCED *NS100430 

C *NS l 00440 

C ********************************** ********+***♦********************** **NS l 004 50 
C ‘NS100460 

WRITE! TYPE WR , 100) NS100470 

ICO FORMAT!* I COO UNSUPERVISEO ECHO FUNCTION 1PHASE 2)*, NS100460 

* • REQUESTED ( NS2SUP ) * I NS100490 

C _ *NS 1 00500 

CALL NS2RDR NS100510 

C *NS l 00520 

C *************** ******* ************************************ ************* ns 1 0053 0 
C *N S 100540 

C CALL PROCESSING ROUTINE THEN PRINT TERMINATION MESSAGE *NS100550 

C *NS 1 00560 

C ********************************************************** *************NS 1 00570 


0011 

r 

COMMON /NS2C0M/ CSET13 

1 POLNAM I 2,60), RESULT, 

2 CSEL! 30) , POLP TR ( 2 ,60 

0012 

u 

r 

RE AL*4 CSET 

0013 

u 

r 

INTEGER+4 POLNAM 

0014 

u 

r 

INTEGER*2 CSEL, POLPTR 

0015 

w 

C 

L0GICAL*1 POLNMl 

0016 

C 

LOGICAL*! POLNMl 


CALL NS2INT 
C 

WRITE! TYPE WR , 9000 ) 

9000 FORMAT! • l OOOOUNSUPERV I SEO ECHO FUNCTION IPHASE2) COMPLETED 
* * I NS2 SUP) * ) 

RETURN 

ENO 


*NS 1 00580 
NS100590 
*NS l 00600 
NS 1 0061 0 
NS 1 0062 0 
NS 1 0063 0 
*NS 1 00640 
NS100650 
NS100660 
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Purpose : Read functions control card 

System/Language : CMS /FORTRAN 
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MODULE ABSTRACT 

NS2RDR reads and interprets all function control cards for NS2ECH0. Also 
a results tape is readied. 


PURDUE UNIVERSITY 
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1220 Potter Drive 
West Lafayette, Indiana 47906 
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1 . Module Usage 
NS2RDR 

CALL NS2RDR 

This section lists the actions taken when the following control 
cards are read: 

RESULTS — TAPE 


RESULTS— FILE 


RESULTS— INITIALIZE 


INTERMEDIATE— TAPE 


INTERMEDIATE— FILE 


CLASSES 


2 . Internal Description 

NS2RDR uses standard card reader logic in using LARSYS system routines 
CTLWRD, CTLPRM, and IV AL in reading and interpreting the control cards. 

After initializing flags and arrays which will convey control card infor- 
mation, NS2RDR functions in a loop of reading and interpreting control cards 
until DATA or END card is read indicating the end of function control card. 
After control cards have been read some checks are made on the data. If 
the classes are pooled, the pooling is checked for validity. Then the 
results tape is mounted and, if requested, initialized. 


- The variable RQTAPE is set to the 
given tape number. 

- The variable RQFILE is set to the 
given file number. 

- The local flag INITFG is set to TRUE. 
This flag is used by NS2RDR to call 
MMTAPE to initialize a tape. 

- The variable OTAPE is set to the 
given tape number. 

- The variable OFILE is set to the 
given file number. 

- LARSYS System Subroutine POLSCN is 
called to set up the arrays POLNAM, 
POLPTR, POLSTK, and P0LNM1 and compute 
NOPPOL and STKPTR based on the inter- 
pretation of the classes card. 


3. Input Description 

Function control cards for NS2ECH0 are read via LARSYS System routine 
CTLWRD. 
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4 . Output Description 

Control card error messages are written to both the printer (PRNTR) 
and the console (TYPEWR) . A brief list of these follows: 

ERROR IN CLASSES CARD. CORRECT ALL CLASSES CARDS AND START OVER. 

ERROR IN RESULTS CARD (TAPE OR FILE PARAMETER) - TYPE IN CORRECT CARD. 

Syntax error in the TAPE or FILE specification for the intermediate or 
results tape. Standard corrective action. 

A POOL HAS NOT BEEN DEFINED. CORRECT CLASSES CARD. POOL NUMBER IS 

The pool number is written on the next line. The function terminates. 

Pool numbers must be consecutive and start at 1. 

BOTH FILE AND INITIALIZATION OPTION REQUESTED. FILE REQUESTED IGNORED. 
FUNCTION CONTINUES. 

The results tape is initialized. Only file 1 can be initialized. 

UNEXPECTED END OF FILE ON INPUT DATA 

The end of the input deck was reached while reading cards for the function. 
This normally means that the END card was omitted. 


5 . Supplemental Information 

See LARSYS System Manual for card reader requirements. 
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F3JT3ETT IV3TTVEL 20. 7 
f 1LE :$2*0* 


OATF »" 7 7/79 


OJ\*l 

O'XV 

OvM*» 


0011* 
OJf‘5 
0 wilt 

oou/ 

OOUb 

OOfT? 

OUltl 


SUHAOUllNl NS/ROH 
ixptiCfi i\7h;ir»* m-/) 


■NflM /Cl l CUN/ PLAN*. CARIM/01. CHMIUI , ClN>l| l • Cl ASSN, I.IASSX, 
Li us i x , ciinpui. cpvuui . cmukdk, ckostw, caiapi, 

UUPIIP. .JUI'KUN, 1 PR “SI* • I RI'N I • • 

ROUND, MUAIA. Ill ADC II H I • 101/001, INAt.lX, 

“*1XCHA. MAXCI S 


N M VHO. MAPTAP, .. 

*» I 1 ACS I 2 i PNCH. Mill N r . PHI MM, PR.NTK, 11X011. 

t, »rsr«r, HUNT [ l . HUNT AH( 1U, II , 

/ SO A t A ( SI PARK, STPfPA, SPARLIIOl. lf“PASII01. 

« tPSTAtlfcl. I If (.OX . ITPIM1. 

9 im, ARRAY!) 7S0l>) 

RL *1 • b A (KAY 
KFAL * A f IOCALI 10) 

I.NtEufcR * A CQ»lNllHI. DATEI 5) , HLOI I 161 , HE 0/ I 16) . UPEI5I 

IMtuCH * / IILANKZ.NNIDISOO) 

LU-.1CAL • * LHKOUI 

LOGICAL « | HLA.VKI 

EUU IVACE'lCC I CAT SAV « I 0 ( I II ■ I CURK1N. I CO I I . I THUC »l I II . I D l S |)l . 
t I H ; Ul I I I ,M£ AO ( 8 Jl , IOAIEI 1 I.mEAUI? 6M, IHLD2I 11 ,hf ADI 19) I 

/ I I INfcl I I *1)1581 I . I CORF Ml ] liliUm /} I I I 

1 (H \P SAV . fl'STAf I III. 


^aiNniiioi 


£L»NUU /.'IS/C'JN/ CSFR I, 101, HTKSI/, NllTITl. Ml 
I P(>L'lAM/,(.t)l , ■( L ' UL T ■ RUT ILL. KUIAPL,SrKPIK, 
/ C'.blHOI. PULPIRI/.bOl, POlSfK 

■*1 *L** CSC I 
nTEui.ll** '»OLNA* 

r«tiat«*/ Lict , polptr. polstk 

t LG 1 C A L* I pOLN“1 


**5/00010 
•NS/OOOZO 
NS/00010 
•NS/000* o 
••NS/00050 
•NS700060 

ns/ooo/o 

•NS7O0UHO 

Mpnnnvo 

• NS/nmoo 
N 5/110 110 
NS/001/0 
NS/001 10 
NS/OOIAO 

NS/001 SO 
NS7O016O 
NS/001 7 0 
NS/00180 
NS/00190 
NS/00700 
NS/00/ 1 0 
NS/00// 0 
NS/00710 
NS/OO/AO 
NS/OO/SO 
NS/OO/AO 
NS/oo/ro 
, NS/OO/Bo 
NS/00/90 
NS/00100 
.NS/OOMO 
NS/001/0 
Ns/oimo 

• NS ?odho 


007/ 
0071 
002 * 
Od/5 
00/6 
OJ7.7 
OO/rl 
00/9 
00 10 
aou 
on/ 


LOGIC AL* 1 POL NP 1 


INTEGER** ;DCNR60), KEPPIS160), REDOUFT/Ol, 
INTEGER** ^CY-ROIS) .RESCODI 1) 

LOGICAL** IN1TFC 

INTEGER* 2 olNFLG,f ETVCCUOI , 7CTVCJI10I 
0 A 1 A K EYhRL /*CLAS*« * RESU ' * * KITE ' .'DATA* , ’END 1 
" sf *»*FlL|*t , lNIT'/» 


INI »FC * .FALSE. 

sroprG * o 

RESULT * NAPIAP 
RUfAPf 


JTAPC « 0 
UTILE • 0 
OFCIl * 0 


AOFC. . 
.NOPOOL ■ 
SfKPTR = 
STAfAY * 

no iu i* 
CSCLIT » 



FORTRAN |V G UVTL 20.7 


00 18 
flu 19 
00*0 
00*1 


oo*s 

00*6 

00*1 

OOAB 


III FL I VC »ll» = 0 
c o i i i • i , so 
POL STH III 1 '! 
Atpprsrn * o 
PuL W* 1 I I 1 * .THUF . 
PULPTH (1, 1 1 >0 
POLP I R (/ • 1 I =■ 0 
00 IS l»i , 1 
DO IS J*l , 10 
IS CSf III ,JI*-SOOOO. 


LCCULl CUNIKCl CAROS 


OCCOUL ACY^ORO AMO GO ACCORDINGLY 


NS/oorio 

NS/OO 7/0 
NS/00 7 >0 
S/007*0 
S/00 MO 
'200760 
.200770 
S/OO/BO 
S/00790 
S/00800 


NS/00810 

NS/008/0 


• NS. 


assign loo ri covic n£. 

Cl 1,0 ic l«OVfcC, 1100, S*Ol Nf 

IUO CALL C ILPRIilCARO.COL.KT VWRD, S , CODE . S, ERHCOR i N 

iricopi - e «.. s ) spare i /) = i m 

IHCRPCOH .tO. *1 CALI KTNAIN N 

ITU^KCOR .14, 2 I CALL FRPRNU 119, •STOP’ I 


2008/0 
.200810 
$2008*0 
S/008SO 
S/00860 
S/00870 
S/000110 
S/00890 
S/00900 
S/00910 
S/009/0 
S/00910 
5/009*0 
S/009SO 
S200960 
S/00970 
S2 00980 


0050 

0051 
00$ 2 
0053 
005* 


005b 

0057 

OOSU 

0059 


007* 

0.lt5 

Odbfe 

03*7 

0068 


0069 

0070 
00 M 
00 7/ 


111) URl TE I TYPE hR, 91 IOICaro 
9110 FORMA T 1 5* , /O A* I 
\> CALL E RPRN 1 1 ERRriUP, * GO *1 
II IFIE1RC0R . Nf . 31 ERRCOR = 1 
CO TO 100 


LASSES CA.NO 


CALL PCL SC Nl POLNA0, 
• CARC, C/10) 

GC i": 100 


POLPTR, POLSTK , PCLNM1, NO POOL , StRPTR, COL , 


Li SUL T TAPt ESPtClF IC4TI0NS 


TAPE SPfC I F1CATICNS 

VLCS : - l 

CALL IVAUCKO, CCL . RQTAPE, VECSI. C7I5) 
IFIVECSZ . zl. Ol GO TO 115 
GO 1 C 100 
ERRNUM * 1*4 
uT I" 110 

FILE SPECIFICATION 
VCCSl * 1 

CALL IVALlCAhD, CCL, ROFILE, VECS/. C115I 
IflVCCSI .CJ.OI GC IQ 115 
uG Til 100 


NS/O 

NS/O 

NS/O 

NS/O 

NS/O 

NS20 

NS/01 

NS/OI 

■NS/O 

NS/O 

NS/OI 

NS/OI 

•NS20 


0/0 

010 

0*0 

oso 

060 

.070 

080 

090 

100 

110 

I/O 

110 

1*0 

ISO 

160 

170 


NS/O 

NS/O 

nSzo 


IFICOL .£5. 7/1 GC rc 101 

CALL C1LPR“1CAR0, COL, RESCGO, 1, CODE. UUI 
GO TLI 110, 3/0. 3)01, CODE 


■NS/O 

■NS/o: 

ns/o; 

nS/oi 

NS/O 

NS20 

ns/o: 

NS/O 

NS/O 

NS/O 

-NS/O 

NS20 

NS/O 

NS/O 

NS/O 

NS/O 

NS/O 


200 

210 

220 

210 

7*0 

250 

260 

270 

/BO 

290 

100 

310 

I/O 

110 

1*0 

ISO 

160 

170 

180 

>90 

*00 

*10 

*/0 

1*10 

***0 

.*50 

1*60 
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INITIALIZE SPECIE ICAt l OM 


I NT tRMIUI A I E TAPE SPEC I F I C AT I CNS 


1 FICCL .tO. J?) CC TO 100 

CAll C ILPR“ICARD,COL . HE SCOD. ?. CUCT , tl 

GO 101*10, *?U», CODE 

TAPI SPECIFICATIONS 
VtCS.' * I 

C At l IVALICARO, CCt, OTAPL, VTCS/, LAIS 
ir i v. c sz . . o. n » r.o TO *i*» 

GO 111 *00 
IRKNUM » »*U 
GO TO 110 


VFCS/ » I 

LAtl IVAl ItAttO. CCl, 01111, YIC52, L*|-,l 
II I VI C SI .10. (>) GO ID * IS 
GO III *00 


•NSZO * T 0 

• HSZO *-10 

• NS70 *oc 
NS20 500 
NS70 510 

• NS20 5 20 
••NS20 SJO 

• NS20 5*0 
•NS20 55C 

• *1 $ 7 0 SAC 
••NS?n s 7 c 


cur CH I OH PHI si NCI AIJO VALIDITY Of AIL INf lIRMA T ION IN CDHI 


IFISTOPFG .EO. 1) CALt RTMAIM 
IF (NOPOOL .Lt. 01 GO TO 5*0 
00 575 I • 1, NOPCOL 
TfIPOLPTRI l.h .EG. 0 )G0 TO 510 
CONTINUE' 

GO n: s*o 

CALL ERPHNr(J**,»r.O*l 
HKlTF (PRNTH.NSJO) I 
WRI TE I TYPE WR » 95 10 > I 
FORMAT I* *,1*1 

CALL H IMA IN 


CHECK RESUlTj Ci 
INI T WILL HL SUL I 
FILE WRI TTIN ON 
HECOHO 


CARO ENTRICS. ANY FILE t.NIRY WITH RCSULTS 
ILT IN CHANGING THE FILE NUMBER OF THf FIRST 
IN THE INITIALIZED TAPE. NS7INT WRITES THIS 


IF I .NOT. INITFGI GO TO 5*5 
IFIRjF ILF .CW. 01 GO TO 5*7 
HR I It I TYPE wR .95*01 
WRITE IPR-NT 1,05*01 

FORMAT | • 10050 BOTH FILL AND I N I 1 1 AL 1 1 At I C.N OPTION REQUESTED. 
• ‘FILE REQUESTED IGNORED. FUNCTION CONTINUES I NS2RCR > * ) 

HOF IlE » 1 

CALL MPTAPFIROTAPL, T), II 
GO IJ SOO 

CALL MMTAPClROTAPl.RgFILE.II 
IFIRUFVLE .GC. 0>GO TO 900 
ASSIGN 5*0 TO GOVEC 
ERRCUR * J 

ROFILE * IABSIROFILEI 
GO Tu 100 
Kt T URN 
FND 
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MODULE ABSTRACT 


NS2INT carries out required initialization of the rest of the variables 
used by NS2ECH0 and finishes reading the statistics to be used by the 
processor. Also the intermediate tape is readied and the subroutine 
that performs the classification is called. 
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1 . Module Usage 
NS2INT 

CALL NS2INT 

The program NS2INT is called with no arguments. Any variables to be 
used or changed are contained in common blocks GLOCOM and NS2C0M. 


2. Internal Description 

NS2INT performs the following functions: 

1. Uses LARSYS system routine STAT to read the statistics from cards. 

2. Computes space needed for arrays to hold original statistics. 

3. Use LARSYS system routine CLSCHK to check class validity. 

4. The intermediate tape is mounted and positioned via TAPMNT. 

5. With the information read from the intermediate tape it computes 
array bases for reduced statistics. 

6. Reduce statistics using LARSYS system routine REDSAV. 

7. The first three record types are written on the results file. The 
first record has a 0 in the sixth full word to Indicate the absence 
of weights in the file. 

8. Some processor Information is printed out. 

9. With the information from the Intermediate tape the array bases for 
calling N2ECH0 are computed. 

10. N2ECH0 is called for performing the classification when finished 
it writes the needed tape marks and the check record. 


3 . Input Description 

Statistics are read via a call to LARSYS system routine STATS. Infor- 
mation from the ID record of the intermediate tape is obtained via a call 
to LARSYS system routine TOPRD. 


4 . Output Description 

Information concerning serial number, classes, field and channels is 
written on the printer (PRNTR) . The results tape is written using unformatted 
FORTRAN write statement. 
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Supplemental Information 
Not applicable. 



6. Flowchart 


NS2INT 









F UR TRAN IV 
FILE 'IS? INI 


LLVEL 20./ 


t»At,E 0P01 


Ouo* 

0005 

0006 
000 / 
0000 
0:>Ci 
Oo to 


0019 

00?0 

0121 


002 * 
OH?* 
0026 
00 2/ 


SUPR.)UTINF ’IS2INT 
1-i‘LICTT I Ml *»FK •* IA-71 


CO«m:*j /biccur/ blank, c&«r.i/o), chkfiui , Ccptil, classr. 

1 CLUSTX. CONPUT, CPYOUT. CROROR. CROSCO, 0»tA»l . 

2 IJPLTP. OUTRUN. tRRMSC. rO'*NI. 

I niciv, riDHNU, *'Dir«, iiCAOItlbl. lut?oo), imaoex 

I PARK, KEYBO. MAPTAP, «AxCHA, “AxCLS, 

* If R. RfcAl 


PAGSlK RVCHi PCI *11 1 PHE SUX • " PRY T I 

restrTi runfil. runtabi u>,ii , 


fc AD I N . 


. F LOX ■ FYPFWi, 

» TOP. A 4RAYI I2SQ0) 
fif A L • It ARRAY 
REAL • A F :0(,AL I S, 10) 

INI toe R * * COPEN » I lb) , OAT E I * ) i HEUIU'j), HE07I16I. TIP 
INIEStP * i OLANK? 

LOGICAL • ■« CHKOUT 
L CO 1C A l « I BLANK! 

tuUIVALENCL (OATSAV, ll)( l ) ) , ICU-U III. 11)1 t) I , (FKOCALlll.! 

1 (HE31 C I ) .HEADIAM , I 0 A IE ( | | , HEAP I 26 I 1 , IHCU7IU, 

2 IT1HEI II .HCADI5B) I , I COMtATI I I.XCAOI T?> I . 

» (PAPSAV.TPSTAIM )l , 

- I S..PSCR • IP STM 12) I . ICUPIA.TPSTATI J> ). ICASIAT.T 

0 1C IPSFR, TPSTlIlS) I , I TR AUL T , I PS I A I ( 6 I 1 , 

b (Blank. blank? .llanm i 

T*.£ 


CL“"UN /NS2C3*/ CSE1 1 1 . 10 1 ■ '■TXSIZ. NCFETS, NOPCCL , CFILE. OtAPt. 

1 POlNAttl?. 60) . RESULT. RUFlln. RJ I APF , S IKP |R , VECSlZ. KKKSIZ. 

2 CSfcllSO). POL PTRI2.60I. PCJLSTKI60I. PJLMP1I60) 

REAL** CSE f 
INTEGER** RJL’IAP 
IMfcuC R*2 COLL, 

LOGIC A l* I PGLN“I 

LOGICAL*! POLMHI 


•NS700UI 0 
•NS200020 
*152000 10 
• *152000* 0 
••NS700H50 
•.’4S700060 
NS700070 
•NS200010 
JS700Q *0 
•NS2O0IOH 
NS7O0U0 
VS7 00120 
.15200110 
n S 200 i*n 

NS200I50 
•’152001 6U 
-N 5 2 0 0 1 / 0 
M 5200 IDO 
NS? 00 l NO 
NS200200 
1JS200210 
NS200220 
1452002 10 
NS2002*0 
NS2002S0 
NS200260 
.’45200? TO 
I . ‘4S7007HO 
NS?002**0 
NS200 TOO 
I . NS200 1 1 O 
.*4 52 00 120 
NS700 1 to 
•’45200 t*0 
44S700350 
’4 N $200160 
NS2004/0 
•NS2001BO 
NS200190 
445200*00 
•.’45200*10 


PQLPT.R. POLSTk 


INItW *•* \OCNTI60). ... 

INTEGER*? hi IfLG.FEWCCI 101 ,f ET VC 3 IJO l 
LATA EC.i/ • lOS ' / **4“ ’ ’ ’ 


KEPP IS ( 60) • _ REDPUF ( ?0 I . PREFIX 12 I 


RE AU AND SAVE ST Al 1ST ICS 


NUFLUt * 0 
NUFLG5 * 0 
UNIT * CRD (Cl 

IFICARtm .LO. END I CALL ERPHNK75*. ‘STOP* 


.NO SUAfA 
XCLS * NCCLS 
1SIZ = NCFEAT* 


40FEATM )/2 


•NS200510 
•NS200520 
U 52005 10 
NS?005*0 
NS200550 
•NS200560 
••NS200S 70 
•N52005BO 
NS200590 
•NS200600 
• •NS2O051 0 
•NS200620 
NS200610 
NS200**0 
•4 S 2006 50 
NS20Q660 
NS20067Q 
•1S7O06H0 
NS700640 
NS700/00 


FORTRAN IV 
FILE *452 11T 


Ou J2 
001 1 
OJt * 
00 15 
00 16 


G LEVLL 20.7 


00 19 
00*0 
0U*I 
Ull*2 
00 * ) 
00 ** 
00*5 
00*6 
00*7 
00*11 
00*9 

0050 

0051 

0052 


OUSS 
005* 
0055 
00 56 


0061 
0062 
006 3 


CL VAR I 
AVAR I 
CLS10I 
CATBA5 

coat ■ 

IT (COR 
CALL R 
• Hi 


COVAR!*(VARSIZ*HAXCLS * 11/2 
3 A V Ail I . I NOF F AT *H AXCLS* I 1/2 
« CLSIlll ♦ HAXCIS 
TOP -IOATBAS-T )*8 

’ CRPRNfl 171, ’STOP* I 


E -LC. 01 CALL CRPRNTI t7l, •$ 

EOS? i (ARRAY IC0VAR1 I , AAR AY ( A V AR 1 ). AKRA Y ( CL S I 0 I I .KEPPTS. 
NFLC, NUCNT. NOCLS, CSCT, CSEL. FllVtC. VARSIZ. NOFEAII 


CHECK CLASSES VALIOI TY 


n.’icis. nopcm . pih -mr, PCt STK. polnai 


•UUH THE INIfRPf 01 AIT RCSULIS TAPf 


pure • o 

CALL r APHNMUI APL . 0* II E , MOiitl 
CALL I FPRO II1AT API , fiOO . IAkOR, III. 7. 
KUFtl i ^ 101 ’>1 

vfcsiz* nor c r t 
PTASI/*VlCSIZ*ttfftSI/*l 1/2 
Cl; I/O K « 1, ’4EIF t T 3 
FI FVCTlKI - NNIOI 19*?»K1 
L • I 

:c i to <7 = i , io 

L Sb L I K ? I • 0 

I F ( FE T VC 3 ( L I .NE. K2I GO FO ISO 
CS*L (K2I « .’4.’4 1 0 I *0‘ ?»L ) 


CALCULATE UA5C5 OF REDUCED ARRAYS 


VARSZ 3 =* NCFLI KINOFCI 1*11/2 
CCVAR3 * CLVARI 

AVAR! * CQ VAR 1 ♦ I V AS S L 3 *.NOPOCL *1)/? 
SCALAR * AVAR) ♦ I NOF E F 1*N0P00L ♦ l » /2 
CETAR . SC MAR ♦ (NOPOOL ♦ 11/2 
CGNAR = OEIAA ♦ ( NOPOOL •!) /2 
•GRkOS * CCNAR ♦ f NOPOUL • i I /2 
rEOUCE Fm£ STATISTICS 
CALL R I 

• JISJYhyjiii. h'i-II, .-4UPLUL , ’IU 

• PCI ST x. F.1VEC. f E TvCI » VARSZ1) 


PRINT SUPERVISOR I \F 3?**A T I CTJ 


L L G ISE.RL (SERIAL | 

; I IE I FRNT i , 100 <1 

|R«JTI//' NO'JSURCRVISEO ECHO I PHASE? ) IHFCRPA! ION... . 
ll ICIPRNM.inollSERIAL. NO®CCl, NOF E It. ’lOFLOS 
IR-lTISx. 'CLASSIF ICAIIO’. STUCY.... *.(( 

/>X,'v.;. CF CLASSES *,l* 

/5 * . • NC. zr CHARNELS *.l* 

/Y».'*lj. CF FIELDS '.1*1 

I It I PINT 10021 (FEIVC3I I). I ■ 1 • NO r ITS) 

iRpAi i5x, 'C hannels selectec are miit.isiii, 


NS200710 

NS700720 

T.S2O0M0 

NS?007*0 

NS200750 

NS20076O 

NS200770 

NS2OO70O 

•NS200790 

•*NS 20 oano 
•NS200H 1 0 
• NS200P20 
•NS700R 10 
•*NS? 0 OB*O 
•NS700H50 
. NS700B60 
NS200H70 
•NS700HR0 
• •NS200B90 
•N5200900 
•NS7009 | 0 
•NS200970 
• *N5200 v TO 

•NS70Q9A0 
NS700950 
NS?0n-»60 
NS? 00 I/O 
NS70O9HO 
:iS700'«0 
•4S70IOOO 
NS70I0I 0 


NS70: 
NS20I 
4IS20I 
I4S20I 
NS20 1 
NS20 
’4S20 
NS20 
•NS20 
••NS20 
•NS20 
•NS?n 
• NS20 
••NS20 
••4520 
NS?n 
•NS20 
NS20 
N5?0 
•i 570 
.’4 520 
NS20! 
•*IS?0 
•4520 
•NS20 
ns?o 

NS20 
74520 
•NS20 
• **4S20 
•NS20 
•NS20 
•NS20 
••NS20 
•4J52P! 
N520I 
'•520! 
NS20 
45201 
*4520 1 
N520I 
NS201 
NS20I 
• S20 1 


020 
u to 
0*0 
05 0 
060 
070 
00 0 
OlO 
too 

no 
120 
I 30 
1*0 
150 
IbO 
1 70 


2 to 
2*0 
?5Q 
7*0 
2/0 
700 
7 tO 


*70 

*30 

**3 

*50 
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FORTRAN IV 
FILE NS2INI 


0068 
0069 
OU70 
00 71 
01)72 
0075 


i LEVEL 70. r 


DAtt * 1722'* 


.0/4 

0075 


0090 

8831 

009 l 
0U9* 

0095 

0096 
009 7 

0098 

0099 

0100 
0101 
oio; 
0105 
010* 


hRITE riRST RECCJRO On RESULTS 7 APE 


PREFIK2I * 0 

bRI TE ( RE SULT 1 PREFIX, ROT APE, ROf ILL.U. ID. SERI AL . ( 1 0. I > 
REWIND SO* T A 
PREFIXM] * * 

UR] TEIRESULUPREFlX.NO.NOFEri.NOFLCl.NDPCCt.lFETVC M I 

• ^MPOLPTRII j J l, l»l.?l • J* 1, NOPOOL ), (POL SIKCli, IM.NCJ 

RE AO { SO At A <9220 « ENO» 300ICARO 
FORMAT i/QAS) 

uKI IE (RCSUIMPRIFU.CAHO 
I F I CARO I I ! .NE. EUSl 00 TU 270 
CONTINUE 


•NS20I 
• •.NS20. 
•NS201 
•NS201 
•NS2flr 
••NS20 
•NS20! 
NS20I 
NS20 


S20 




l-l .NDFCT.3 JNS20! 
UNfiFCTlU NS20I 
l .NQPQOL I N S ? 0 
L NS20 

■ 1520 
NS20 

r*s?o 

NS20 

NS2(1 

MS/O 
s?o 
s/n 



500 
510 
520 
530 
5*0 
550 
56C 
570 
580 
590 
60 0 
610 
620 
630 
6*0 
650 
bbO 
670 
60 0 
69 0 
/DO 
71 0 


(Nil* I T I * I NUT I I I 
IV* CM2 • 11/2 
INIKS1/ . ||/? 

I V( LSI 2 • NOP 0(11 
I VI CS I / * 7 . 11/ 
INnPUOl *11/2 


NUT C T 3 • I NQf t (!♦/! 


CALL N2ECHUIARRAV (COVARI) . AR RAY I AVAR 1 1 , AMR AY I SL Al AR ) , 

• AMR A VIDE I AR I | ARRAY I CUN AH I , ARM AY I WIIM7 US I , AMR A 1 1 1 A VH AS I , 

• A HR A Y t f C VB AS I « AMRAVtSAVRAt). APR AY 1 III OR AS ) . ARtt A Y I CH »»A N > , 

• ARKAY (OATHS II. OATHS!! 


*15201 t 
NS20IHOO 
NS/OIIIIO 
NS/<> 10/0 
*5201010 
♦*57010*0 
•NS20IM60 
NS70I060 


N52* 


H 70 


HIM THE l AS T FILL ON R! SUL I S TAPE 


FILESV « F lit SV ♦ I 
COPFIL ■ C CPF I L ♦ I 
PREFlAlll * l 
PKEFIII2I * 0 
WRI TE (RE SUL 1 1 PREFIX, 
CALL TOP(F IHCSULTI 
CALL ICPRF (RESULT T 
CALL TCPtF (RtSULT I 
CALL IOP0F (RESULT I 

call tcprf ire Sul t i 

RETURN 

END 


ROTAPE, FUCSV, 


i K, (TO, J-* 1 ,i 


NS7O1H0O 
N520 | H'»0 
•NS/O 1900 
••NS/uj'UO 
•NS20I920 
NS/o i 'M n 
•NS20I960 
••NS701950 
•NS 70 1 960 
NS/O l '» 10 
NS2019H0 
NS/01 990 
NS/O/OOO 
NS/02010 
NS202020 
NS/02010 
NS2020*0 
NS202050 
NS/02060 
NS2O2070 
NS202080 
NS202090 
NS202100 
NS/02110 



LARS Program Abstract 

MODULE IDENTIFICATION 

Module Name : N2ECHO Function Name; NS2ECH0 

Purpose : performs classification for NS2ECH0 and writes results tape. 

System/Language : CMS /FORTRAN 

Author Date : 

Latest Revisor: C. A. Pomalaza Date: 8/21/77 


MODULE ABSTRACT 

N2ECH0 is called by NS2INT to perform the field classification on those 
fields and singular points identified by NS2ECH0 (field extraction phase) 
and wirtten on an intermediate tape. N2ECH0 writes a standard results 
file to tape. 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 
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N2ECHO-2 


1 . Module Usage 


N2ECH0 

CALL N2ECH0 (COVMTX, AVEMTX, SCAPRD, DETCOV, CONST, WORK, MEAN, COV, 
SAVAVE, HOLD, CHISQR, RDATA, DATBS3) 

Input Arguments 

COVMTX R*1 Array of covariances obtained in NS2INT. 

Dimensioned NTXSIZ*NOPOOL where NOPOOL is the 
number of classes contained in the reduced STAT 
DECK and MTXSIZ = VECSIZ * (VECSIZ + l)/2 where 
VECSIZ is the number of channels requested for 
classification (in Phase 1). 

AVEMTX R*4 Array of mean vectors for each class. Dimen- 

sioned VECSIZ*NOPOOL. 

SCAPRD R*4 Array of the scalar product of the mean and 

covariance matrices for each class dimensioned 
NOPOOL. 

DETCOV R*4 Array of the determinant of the covariance matrix 

for each class dimensioned NOPOOL. 

CONST R*4 Array of the constant terms used in the computa- 

tion of the likelihood values dimensioned NOPOOL. 

WORK R*4 Working array for a subroutine cell dimensioned 

VECSIZ * (VECSIZ + 2). 

MEAN R*4 Array used as a holding buffer for a subroutine 

cell. Dimensioned VECSIZ. 

COV R*4 Array used as a holding buffer for a subroutine 

cell. Dimensioned MTXSIZ. 

SAVAVE R*4 Array used for saving the mean from each class 

of the STAT deck. Dimensioned VECSIZ*NOPOOL. 

HOLD R*4 Array used as a holding buffer for subroutine 

call. Dimensioned 7*VECSIZ. 

CHISQR R*4 Array used in a subroutine call. Dimensioned 

NOPOOL. 

RDATA R*4 Area to be used for dynamic allocation. This 

must be the first unused element of ARRAY in the 
calling program. 
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DATBS3 1*4 Number of bytes on ARRAY which one in use 

(i.e. the number of bytes in ARRAY which precede 
RDATA) . 

Output Arguments 
Not applicable. 


2. Internal Description 

1 - N2ECH0 calls to LARSYS system routine SMMULT and SAMINV to invert 

the covariance matrices in COVMTX and produce the determinant in 
DETCOV needed for classification. 

2 - A loop begins where each field statistics matrix (record type 4) 

on the intermediate tape is classified by LARSYS system routine 
SAMCLS. The absolute field number (see tape record description) 
indexes the class. 

3 - The records that associates each pixel with its appropriate field 

are read. The processor classifies each pixel in an homogeneous 
field by looking up its absolute field number. Singular pixels 
are classified by a call to CLASS. The results are written line 
by line following the standard LARSYS format (see RESULTS FILE 
description in LARSYS System Manual). 

4 - The RESULTS FILE is closed and control returns to NS2INT. 


3. Input Description 

N2ECH0 reads field statistics information and pixel data from the 
intermediate tape via call to LARSYS system routines TOPRD, and GADLIN. 


4 . Output Description 

A standard results file is output to tape. See LARSYS System Manual 
for format description. 


5. Supplemental Information 

See non-supervised ECHO intermediate tape description. Also Results 
Tape description on LARSYS System Manual. 
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6. Flowchart 


BEGIN 


/WRITE CLASS COVARIANCE 
'MATRIX AND MEAN VECTORS 
ON RESULTS TAPE i 



Winter- > 

MEDIATE 
.TAPE J 


I READ STATS FOR 
1 FIELD K FROM 
'INTERMEDIATE TAPB 


ISAMPLE CLASSIFY 
FIELD K 
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FORTRAN IV 0 LEVEL 
FILE N2ECHO 


20. r 


N2ECHQ 


DATE » 77229 


1505 2052 


FA GC 0001 


"vl 

Ln 


0001 

0002 

0005 


000 A 

0005 

0006 
000 7 
0008 

0009 

0010 


0011 

0012 
0015 
001 A 

0015 

0016 


0017 

0018 
0019 
00 20 
0021 
0022 
0025 
002* 

0025 

0026 

0027 

0028 
0029 
00 *0 
oo n 

00 52 

0035 

0036 


N2ECM0 


LARS XXXX 

WRITTEN BY C. A. POPALAZA 


SUBROUTINE .\2CCHC(C0VMX, AVEMTX, S CAPR Q, DE TCOV, CONST, WORK, 
t HEAN.COV, SAVAV6 .HOLDtCHiSQR, ROA T A 1 0 A TBS3 ) 


implicit integer** u-zi 


OAT E I 5 1 • HE01 ( 16 ) « HE02II6), TIMEI5) 


COMMON /GLCCOM/ BLANK. CARD(20>, CHKOUT, COPFIL, CLASS*, CLASSX, 

1 CLUSTx, CONPUT, CPYCUT, CRDRDR » CRDSEQ, OATAPE, 

2 OUPIT?. CUPRUN, ERRMSG, FRPNT, 

3 FllESV, FLC3N0, HD AT A * H£A0(88!t 10(2001, IMAGEX, 

A I MARK , KEVOO* MAPTAP, MAXCHA, MAXCLS, 

5 P AUS I t $ PNCH, POINT, PRESUX. PRNTR, REAOIN, 

6 KESTRT, RUNFIL, RUNTA8( 10,3) , 

7 SOATA. SEPARX, SCPTPX, SPAREUO), TEMPAS ( 30 ) , 

8 TPSTAT (6) , TTfLDX, TYPEWR, 

9 TOP, A:uAY(12500) 

REAL • 8 ARRAY 

REAL * A F HOC Al (5,501 
INTEGER * A CCJ V L'NT ( 16 
LOGICAL * A CHKOUT 
INTEGER * 2 HLANK2 
LOGICAL * I CLANK I 

CQUIVALENCt ( 0 A T S A V , 1 0 ( 1 1 1 , ( CURRUN, I 0( 3)1 , ( FROG AL 11) , IOC 51) 1 , 

1 (HCOIU) t HEAD(8U * (OATEl l> ,HEAO( Z6U » ( n£D2( 1 1 , HE ADI 39 1 

2 ITIMU U .HEA0I58) ), ( COME NT C 11 , HE AOt 72 ) ) , 

3 (MAPSAV.TPSTATU) ) , 

6 f SEPSCR, TPSTAT (2 ) ) , ( DUP 1 N, TP ST AT ( 3) } T < OASTAT .TPSTAT l * 1 

5 CUP3ER,TPSTATI5>) r l TRAOUT, T PSTAT I & > > , 

6 (BLANK. BLANK2.BLANK1) 

COMMON /NS2C0M/ CSET(3,30) f MTXSIZ. NOFET 3* NOPOOL, OFILE, OTAPE 

1 J* e fy LT K R0 £ ILE i rotIpe.stkpIr, VECsfz, wrkJiz, 

2 CSEL ( 30) , POLPTR (2,60) , P0LSTKI60) , P0LNM1 ( 60 1 
REAL** CSET 

INTEGER** POLNAM 

INTEGE R*2 CSEL, POLPTR, POLSTK 

LOGICAL*! POLNMl 

LOGICAL*! POLNMl 

DEFINITION OF LOCAL VARIABLES 

REAL * * SPAR( 10) ,AR( 125) 

RFAL* A C SE T 3 ( 3,39) 

INTEGER* a ARK 125) 
tNffc-uEH* 2 A , 0 t AK2 ( 250 I 
INICGFR*2 CSEL 31 50) • BLOCK (6 ) 

LOGICAL* l LR(8), RSI2) 

LOU I VALENCE ( AR ( 1 ) , AR l ( l ) • Aft 2 ( 1 ) ) 
ecu i vale nci turn. I R I I ) ) 

REAL *6 MLANIVECSIZI ,COV( MTXSIZ) 

INTEGER** REJBUFI20) ,IR(2) 

# iNffc^LR*? 3 INF LG* NEGCNE , I N, FETVEC 130), 
t FEfVCJMO) 

COU I VALCNCC (IN,RS(1I) 

INTEGER** INFO! I 7) ,PR£FIX<2) 
iNTEGfc R*2 UUFFEh(2500) 

INTcGtR • 2 SJRUF (22500) 

RE AL ** ROATA( l ) 

INTEGE R*2 TAPCUFI3000) 

REAL** CQVMTX(MTXSIZ,NOPOOl), AVEMTX I VECS I Z, NOPOOL) .WORK (WRKS1Z) • 
6 SC APR C I NO POOL 1 ,CH1 SOR ( NOPOOL ) , OETCOV ( NCPCOL 1 • CONST (NGPQQL ) * 

& HOlO I 7, VECSi 2 ) * SAVAVE I VECSi Z, NOPOOL 1, THRTA&(23*>, PCREJ 


•CCCOOOIO 
GCGOO’JZ 0 
GC000030 
♦GCUOOOA 0 
♦♦GC000050 
*CCQ00060 
GC00C3 T 0 
CCUOOQBO 
•GCU00070 
GCOOOIOO 
•GCOOO 1 l 0 
GCUOO l 20 
GCOOO 1 30 
GCOOO 1* 0 
GCOOO 1 5 0 
GCUOO 160 
UCOOOi 70 
GCC00180 
GCOOOIOO 
GC000200 
GCOOOZIO 
GC000220 
GC0002 30 
GCOOOZAO 
GC000250 
GC000260 
GC000270 
GC00028 0 
1 , GCOOOZSO 
GC00C300 
GC00031 0 
) , GCOOO 320 
GC000330 
GC000J60 
•GC00035 0 


•GC000A50 
GCOOO',60 
•GCCOOA 7 0 
GC0004Q0 
GC000A90 
GC000530 
CC0095I0 
GC000520 
CC000530 
GCC005*0 
GC000550 
•GCU00560 
GCO0057 0 
CC000580 
GC000590 
GCC00600 
CC000610 
GCU00620 
GCU00630 
GC0006AO 
CC000650 
GC000660 
GC000670 
GC000680 
GC000690 
•GCOOO 700 



FORTRAN IV 0 LEVEL 20 . 7 


N2ECHQE 


DATE * 77229 


13052052 


FADE 0002 


FILE N2ECH0 


*r 


0015 

0036 

0037 
0018 
0039 


00*0 
00 A 1 
00*2 
00'* 3 
00** 
00*5 
00*6 


00*7 

00*8 

00*9 

0050 

0051 

0052 

0053 
005* 

0055 

0056 

0057 
0050 


0059 

0060 
0061 
0062 
0063 
006* 

0065 

0066 
0U6 7 
0068 
0069 
00 70 
00 71 

0072 

0073 

00 7* 
0075 


0076 


WRITE COVARIANCE MATRIX ANO MEAN VECTOR ON RESULTS TARE 


REWIND SOATA 
PREFIX (l )=* 

PKEF 1 X (21= I 

ViRt Tt l RESULT) PREFIX, COVMTX,AVEMTX 
NEGONEs-l 


INVERT THE LEADING STAT DECKS ON INTERMEDIATE TAPE 


CC 510 1 = 1 , NOPOOL 
CALL SAMMvtCOVMTXI I , I ) ,NOFET 3,DE TCOV { I ),WORK) 

IF (OCTCOVm.LE.O.) GO TO 520 
„ CC 7 8 3 U = l, VECSU 
788 SAVA VC IN • I ) = AVEMTX(N.I) 

W UL TI COVMJ X < l , l ) .AVEHTXI l, I ) ,SCAPROm,WORK,NOFET3) 
510 CONST ( 1)=-.5*AL0GI (6. 2831853* *N0FET3) *DE TCCV ( l ) } 


5 i 1 
520 
9222 


the code uelow builds a table of threshold values 

CORRESPONDING TU FRACTIONS REJECTED FROM 0.1 PERCENT 
TO 97 PERCENT. PCEREJ IS 1.0 -PERCENT REJECT ION/ 100. 

PCREJ = 1.0 
CO 511 1=1, 23* 

PCKEJ = PCREJ - .001 

I F ( I .GT. ISO) PCREJ ^ PCREJ - 0.009 

CALL ThRESUPCREJ. THRTAB(I), VECSIZ, ERROR) 

IFIERRCR .NE. 0) CALL ERPRNT 1 377, 'STOP'l 
CCN T I NLc 
GU TO 5 '5 

CALL ERPRNTt I60« v GOM 
Ml TE ( PRNTR.9222) 

FORMAT ( »SI NCULAR CLASS*) 

STOP 


CLASSIFY EACH FIELD 


535 CALL TCPRO C 12, 68, ERROR, INFO, 2.WLF) 

CALL T QPRD 112 , 16, ERROR, SPAR, 2 t«LF> 
CLLwfH = S par ( I ) 

LCOU-N T =8 ♦ ( VECSI Z^MTXSIZ)** 

L l * i D ( * ) 

OU 20 J * l * L 1 

CALL TCPRO < 12 vLCOUNT , ERROR *AR ( 1), 2,WLF) 
A = AK2 ( 1) 

R = AR2 < 2) 

0--A.U( 2) 

oo \ i k2=i , v:cs i z 

11 MLANIK2) »A.<(2*K2) 

CC 12 K2=1,HTXS1Z 

12 C0VIK2 )*AR t 2*VECS1Z*K2) 


CALL S.'MCLS(CLASS1,CHI SQR , A1 , CCVMTX, AVEMTX ,SCAPRD*CONST , COV, 
t 0,VECS1Z.hTXSIZ,N0P00L».FALSE.) 


20 


SOBUF ( J)=CLASSl 

continue 


write RECORD 5 ON RESULTS FILE 
PREFIxm«S 


GC000710 
••GCCP0720 
♦GCC 00730 
♦GCC007*0 
♦GC000750 
•♦GCCCC 750 
•GCG00770 
GCC00790 
GC000790 

Gcnooaoo 

GCCOOBI 0 
GCO0OP2O 
•GC000830 

• ♦CCUC0 ; j*0 
•GCCOCoS 0 

GC000560 

♦ CCCOOo 7 0 
•♦CCDCOCPO 

•GCC00B90 
GCCOOVOO 
GCC00910 
GCU00920 
GC000730 
GCCOO ■** 0 
GCU00950 
GCCOO 760 
•GC0009/0 
♦GCQ00930 
GCOC 0*79 0 
GCC01000 
•GC001C10 
GCDO 1020 
GCC01010 
GCOO I 0* 0 
GCC0I050 
GCG01060 
GCOO 1 0 70 
GCGOIOBO 
GC GO I 09 0 

Gcnouoo 

GCOO l 110 
GCC01120 
GC001 130 
♦GCOO 11*0 
••GCCOl ISO 
•GCCO) 160 
•GCOO 1170 

♦ GCOO 1 l fl 0 

• •GC001 190 
♦GCC01200 

GCC01210 
GCOO 1 220 
GCC01230 
GCOO 12*0 
GCDO l 25 0 
GCOO 1260 
GCOO l 270 
GCCOl 230 
GCOO 1290 
GCCOl 100 
GCDO 1 310 
GCCOl 320 
GCDO l 130 
GCCOl 3*0 
MEAN, GCCOl 350 
GCCOl 360 
GCCOl 170 
GCCOl 330 
♦GCGOl 390 

• *GCOO 1*00 
•GCCO 1*10 
•GCOO 1*20 
•GCCO 1*30 

••GCOO 1**0 
•GC00l*50 
GCOO 1*60 
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FORTRAN IV G LEVEL 
FILE N2ECHO 

0077 

0078 

0079 

0080 


N2ECH0C 


DATE » 77229 


1 3052052 


0081 
00d2 
0083 
00 8* 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 
009* 


0095 

0096 

0097 

0098 
0 0 99 
0100 
0101 
0102 
0103 
010 * 

0105 

0106 
0107 


0108 

0109 

01 to 
0111 
0112 
01 1 3 
01 l* 
01 15 
0116 

0117 

0118 

0119 
01 20 
0121 
0122 
0123 
012* 

0125 

0126 

0127 

0128 
0129 
01 30 

0131 

0132 

0133 
013* 
013$ 


PREFIX (2 ) = l 

PTO=l lNFarJ>-INF0<7>M )/ INFO! 91 
LINQM 1.-;F3I5>-1NF0C*)M)/INF0(6> 

WRl Th( RE SUIT) PREFIX, PTC. LI NQ, INFO, ID, <( CSET I I , J ) , I- 1 , 3 7 , Z , 30) 


CHECK FOR CORE OVERFLOW 


9600 

C 

C 

C»**i 

c 

c 

.c 

C»«+4 

c 

6 50 


NSR a PTQ ♦ 6 
NSD = 10(6 ) 

RBASE - .NSR*VECS 1 1 ♦ 1 
LB = !NSD*VECSIZ ♦ 3)/* ♦ RBASE 
IF (M0D<LB,2) .EQ. 1> 18 - 18 ♦ 1 
UASC * 0ATaS3*3 ♦ LB** 

CORE = TOP- BASE 
I F ( CORE .GT. 0) GO TO 650 
CORE = -CO-IE 
CALL ERPRNH 376 , 'GO' ) 

WR1 TE t TYPEWR, 9600) CORE 
WRITE! PANT A, 9600) CORE 
FORMAT UOX, 16, 'BYTES * ) 

CALL RThAIN 


ENTER MAIN LOOP OF CLASSIFICATION 


600 


620 

610 

C 

c 


CO 600 1 » 
CSE» 1(1) « 
Co 600 J = 
CSET3I J,|) 
K = 0 

OU 610 I « 
1 F ( C SE L ( 1 ) 
K = K ♦ 1 
CSEL3IK) * 
OU 620 J * 


1,30 

0 

1.3 


-50000. 


iei 30 


0) GO TO 610 


CSEL(I) 

1.3 

CSEMIJtK) = FRQCAL ( J*2 , 11 

CONTINUE 

CCNTlULfc 


690 


70S 


710 


BLOCK ( 2 ) = l 
BLOCK ( 3 ) - PTO 
BLOCK!*) a i 
PTS ^ PTO 
P I S6 = PTS ♦ 6 
JPTS-P TS/CELWTH 

eur =**jpts 

COUNT = VECS I Z* ( PTS6) ♦* 

DO 6 )2 K Ha 1 ,L l NO 
BLOCK U) = km 

CALL GAOL INtBLOCK! 1 ) , CSEL3 . CSET3, ID. 12. VECSIZ, NSR, 
» RDATA(RBASE) • RDAT A ( 1 ) , ROLL. ERROR) 

CALL TCPAOt 12, BUF .ERROR, BUFFER. 2, WLF> 

RCPTR=1 

6UFPTH=1 

CU T P T K a 3 

TAPQUF ! 1 ) =0 

TAPbUF <2 )=»KM.|nF0(*|-1 

CO 7'jJ K9* l , VECS 1 1 

UU 703 K10M.6 

FCLDiKlOM, K 9 ) * RDATA! PTS ♦!K9-1)*PTS6 ♦ K10) 

00 700 KCO-l.JPTS 

IF IBUFFcRIBUFPTR) . EQ.NEGONE) GO TO 720 
CO 710 K2» l.CELWTH 
BUFN’JN =BUF F ER < 8UFPTR> l ) 

TAPBUF !QUT PTk ) *SQ8UF (8UFNUH) 

OUTP TR* OUT PTR* 1 
ROPTR«RDPTR*CELWTM 
BUFPTR »BUFPTR** 


GCOO 1 *70 
GCGO 1*80 
GCU0K90 
GCOO l 500 
►•GCC01510 
♦GC00I520 
♦GCL01S30 
♦GCCO 15*0 
►•GCGO l 550 
•GCOO 1 560 
GCCO 1570 
GCCO 1 58 0 
GCOO 1590 
GCOO 1600 
GCOO 1610 
GCCO 1620 
GCOO 1630 
GCCO 16*0 
GCOO 1650 
GCOO 1660 
GCOO 1 6 7 0 
GCOO 1 60 0 
GCOO 1690 
GCCO 1 700 

• GCCO I 710 

• GCOO 17 2 0 
►•GCOOl 730 

• GCOO 17*0 

• GCOOL 750 
•GCOOl 760 

►•GCUO 1 770 
•GCOOl 780 
GCOOl 790 
GC001800 
GCUO 1 8 1 0 
GCCO 1 82 0 
GCCO 1830 
GCCO 1 8*0 
GCGO 1 85 0 
GCOO 1 860 
GCCO 18 70 
GCOO 1880 
GCUO 1 990 
GCUC1900 
GCOO l ? I 0 

• GCCO l 92 0 
•GCCO 1930 

GCCO 19*0 
GCUO 1950 
GCOO l 96 0 
GCCO 1970 
GCCO 1980 
GCCO 1990 
GCC02C00 
GCC070 1 0 
GC002020 
GC0O2O3 0 
GCU020* 0 
GCC0205 0 
GC002060 
GCOO 2070 
GCG02080 
GC002C90 
GC002100 
GC0021LQ 
GC D02120 
GCUO? 130 
GCC02 1*0 
GC002150 
GC002160 
GCQ02170 
GC002180 
GC002190 
GCC02200 
GC00221 0 
GC002220 
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FORTRAN IV C LEVEL 
FILE N2ECHO 


20.7 


N2ECH0C 


OATE - 77229 


13052052 


PAGE 0004 


0136 


GO TO 700 

0137 

720 

00 740 K /= l # CC L WTH 

01 38 


CO 7)0 K 1* l . VEC S I Z 

01 39 

730 

HCLOK ,K U *HDATA(ROPTR ♦ (K3~n«PTS6 ) 

0140 

777 

CALL CLASS {hOLO i VECSIZ»1,NQPOOL,SAVAVE.CO i 
C THRTAP1 

0141 


RS121=LR(2 ) 

0142 


RSI 1 ) =LR<6 ) 

0143 


TAPDUF (OUTPrRl*IN 

0144 


OUTPIR*GUTPTfUl 

0145 

740 

ROPIR*RD J TRM 

0146 


eUFPTR=G0FPTR^6 

0147 

700 

cc-irnue 

0148 


PREFIX ll)«>6 

0149 


PREFIX <2)*l 

01 50 


p r = p t s • 2 

0151 


WR! TC ( \ l ) PREFIX, IT APBUFU) ,J«l,PT ) 

0152 

692 

ccntinle 

0153 


TAPbUF ILliO 

0154 


TAPbUF (2) * 0 

01 55 


PREFIX m»7 

0 156 


Will TC<RcSULT)PREFIX,<TAP8UFU>,J-l, PT) 

0157 


PTC =» C 

0158 


PREFIX (I) * 8 

0159 


LINO « 0 

0160 


NRl TC ( RESULT )PREF IX,PTQ,LINQ* INFQ,I0,CSET 

0161 


RETURN 

0162 


END 


GCC022 30 
GCC0224Q 
GCCJ02250 
GCU02260 
GCC02270 
GCC02280 
GCC02290 

Gcno2 ioo 

GCC02310 
GC002 320 
GCG02 i 30 
GCCC2340 
GCG02350 
GC 00 2 36 0 
GC002 370 
GCC02 38 0 
GC002390 
GC002400 
GC002410 
GC002420 
GCC024 30 
GCC02440 
GCCG2450 
GC002460 
GCC02470 
GCC02480 
GC002490 
GC002500 
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Ancillary Data Files Required 
by the 

Nonsupervised ECHO Processor 
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Nonsupervised Field Extraction Processor Disk Files 


The NS1ECH0 processor writes two disk files during execution. The 
first file BUFFER FILE is written as a sequential file by Subroutine NSECHO 
to store field annexation information row by row during the processing. 
BUFFER FILE is read by RDWRTE during the writing of the intermediate tape. 
The second file is also written by NSECHO to store the field statistics. 
This file is used by RDWRTE to compute the data values on the Intermediate 
Results Tape if the 'MAP' options has been selected (see control card 
description for NS1ECH0). 


BUFFER FILE 

This file contains one record of length m bytes for every line processed 
by NSECHO. The length m equals four times the number of cells/line. The 
format of the record is 1*2 with the following structure. 


BUFFER(2i+l) = 


-1 if cell i is not homogeneous. 

The relative field number of cell i if cell i is homogeneous. 


BUFFER ( 2 i+2) = 


-1 if cell i is not homogeneous. 

Absolute statistics record number (in the STAT SCRATCH 
disk file) if cell i is homogeneous. 
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ST AT SCRATCH 


This file contains one record of 2+NC+NC*(NC+l)/2 words for every 
field opened by NSECHO where NC is the number of channels considered by 
NSECHO. Each record has the following structure: 

STAT(l) 1*2 Relative field number. 

STAT(2) 1*2 Absolute sequence number of the field. 

STAT(3) 1*4 Number of pixels in the field used to compute 

these statistics. 

STAT(3+1) R*4 Means for i th channel if field STAT(2). 

STAT(3+NC+1) R*4 Correlation matrix of the field, 

until 

STAT(3+NC+NC*(NC+l)/2) 

The BUFFER FILE is written sequentially using unformatted fortran 10. 

STAT SCRATCH is written using system support routine DEFINE FILE and unfor- 

matted fortran 10. 
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INTERMEDIATE TAPE FILE 


This file is written on tape by Subroutine RDWRTE in the Nonsupervised 
ECHO processor's Field Extraction Phase and becomes the primary input to 
NS2ECH0 (the Classification Phase). The Nonsupervised ECHO Intermediate 
Tape File is composed of six different types of data records. One each of 
record types 1, 2, and 3 are written on each Intermediate Tape File, one 
record type 4 is written for each set of field statistics recorded in the 
STAT SCRATCH file (one covariance matrix and a vector of channel-means for 
each field isolated by NSECHO) , and one type 5 followed by a type 6 for 
every line processed in the Nonsupervised ECHO Field Extraction Phase. 


Record Type 1 (ID Record) 

This is similar to the conventional 800 Byte LARSYS ID record described 
in the LARSYS SYSTEM MANUAL 1 . This records structure is as follows: 


Bytes 

Format 

Size 

Description 

1-4 

1*4 

1 

word 

Intermediate tape number 

5-8 

1*4 

1 

word 

Intermediate tape file number 

9-12 

1*4 

1 

word 

Intermediate run number 

13-16 

1*4 

1 

word 

Number of fields isolated by the field 
extraction phase 

17-20 

1*4 

1 

word 

Number of data channels 

21-24 

1*4 

1 

word 

Number of data samples per channel per 
line 

25-40 

Alpha 

4 

words 

Flightline Identification (16 characters) 

41-44 

1*4 

1 

word 

Month data was taken 

45-48 

1*4 

1 

word 

Day data was taken 

49-52 

1*4 

1 

word 

Year data was taken 

53-56 

Alpha 

1 

word 

Time data was taken 

57-60 

1*4 

1 

word 

Altitude of aircraft 

61-64 

1*4 

1 

word 

Ground heading of aircraft 

65-76 

Alpha 

2 

words 

Date data run was generated (12 characters) 

77-80 

1*4 

1 

word 

Number of lines in this run 
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Bytes 

Format 

Size 

Description 

81-82 

1*2 

% word 

Number of the first channel used by 
NSECHO 

83-84 

1*2 

% word 

Calibration code of the first channel used 
by NSECHO 

85-200 

1*2 


For each channel used repeat information 
of the half words in bytes 81-82 and 
83-84. The remaining bytes are equal to 0 

201-203 

R*4 

1 word 

Lower limit in Micrometers of the first 
spectral band on the original MIST tape 

205-208 

R*4 

1 word 

Upper limit in Micrometers of the first 
spectral band on the original MIST tape 

209-212 

R*4 

1 word 

The suggested value of "CO" calibration 
pulse for the first spectral band 

213-216 

R*4 

1 word 

The suggested value of "Cl" calibration 
pulse for the first spectral band 

217-220 

R*4 

1 word 

The suggested value of "C2" calibration 
pulse for the first spectral band 

221-800 

R*3 


Repeat of words in bytes 200 to 220 for 
the channels on the original MIST tape. 
The remaining bytes are set to 0. 

Record Typ 

e 2 



This record is 17 fullwords long. It describes the area which has 
been processed by the Field Extraction Phase of the Nonsupervised ECHO 
processor to produce the Intermediate Tape File. 

Bytes 

Format 


Description 

1-4 

1*4 


MIST run number of processed area 

5-12 

Alpha 


Field designation on field description card 

13-16 

1*4 


Beginning line number 

17-20 

1*4 


Last line number 

21-24 

1*4 


Line interval 

25-28 

1*4 


First column number 
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Bytes 

Format 

Description 


29-32 

1*4 

Last column number 


33-36 

1*4 

Column interval 


37-68 

Alpha 

Information from columns 51-80 
description card 

on the field 

Record Typ 

e 3 



This record is 33 words long and stores the parameters used in the 
Nonsupervised Field Extraction Phase. 

Bytes 

Format 

Parameters 


1-4 

R*4 

Cell width 


5-8 

R*4 

Variance test threshold 


9-12 

R*4 

Mean test threshold 


13-16 

R*4 

Homogeneity test threshold for 
channel used 

the first 

17-132 

R*4 

Homogeneity test thresholds for 
channels 

the remaining 
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Record Type 4 


This record is the same as on the STAT SCRATCH file. There are m 
records of type 4, where m is the number of fields isolated by the Field 
Extraction Phase of the Nonsupervised ECHO processor. Each record has 
2+NC+NC*(NC+l) /2 words where NC is the number of channels used by NSECHO. 
The structure of each record is: 


Size 

Format 

Contents 

1-2 

1*2 

Relative field number 

3-4 

1*2 

Absolute sequence number of the field 

5-8 

1*4 

Number of pixels in the field with this 
statistics 

9- (NC+2) *4 

R*4 

Mean value for the 1 st , 2 nc *, . . . NC^ 
channel used by NSECHO 

(NC+2) *4+1 - R*4 

(NC+2 ) *4+NC* (NC+1 ) *2 

Correlation matrix of the field 


Record Type 5 I 

This is similar to a standard data line in a LARSYS MIST tape. If 
the MAP option was specified, the original data has been altered so that for 
the pixels identified as falling in an object, the channel mean of the 
object replaces the raw data value. If the MAP option is not active, the 
raw data is copied unaltered to the intermediate tape. 

Each data record will contain one scaft line of data from ID(5) (see 
ID Record) channels. The first halfword (2 bytes) of the record will be 
the line number. The second halfword (2 bytes) will be the roll parameter 
(which is a number indicating relative position of the roll of the aircraft 
for this line of data). If the roll parameter is -32,767, the data for 
the given line does not exist. If the roll parameter has not been calculated, 
it will be set to 32,767. The fifth byte will be the first sample from 
the requested channel. The sixth byte will be the second sample from the 
first requested channel, and so on through ID(6) samples and ID(5) channels. 

A Type 5 record will be ID(5)*ID(6)+4 bytes long. 

All data for each channel is from the field of view of the scanner 
except the last six bytes. The last six are calibration data in the order 
of appearance. 


1. 

c o 

"0" or dark 

level 

2. 

vc o 

Variance of 

c o 

3. 

C 1 

Calibration 

source C^ 

4. 

VC 1 

Variance of 

C 1 
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5 . 


Calibration source 
6. VC2 Variance of C ^ 

where C - Calibration value i and VC^ - calculated variance of calibration 
value i 

During the reformatting process a record may be had due to tape or 
other errors. When this happens, the data roll parameter and calibration 
points will be set to zero. On good data records all data and calibration 
values will be in the range of 0 to 255 (bit form) with no sign included 
in the eight bits. A data value of 0 to 255 means that the data point was 
cut off during the digitization process. Data values then range between 
0 and 255 with 0 indicating low relative irradiance and 255 indicating 
high relative irradiance. 


Record Type 6 

Identical to a record in BUFFER FILE, i.e. it has m words for each line, 
processed, where me is the number of cells per line. The structure of each 
record is: 

Bytes Format Contents 

1-2 1*2 -1 if the cell is singular, otherwise 

relative field number of the cell 

3-4 1*2 -1 if the cell is singular, otherwise 

number of the statistics record of the 
field the cell belongs to 

5-2m 1*2 Similar to the first two halfwords des- 

cribing the nature of the remaining (m-1) 
cells 
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Final Report Distribution List 
NAS9-14970 


NAME 


NUMBER OF COPIES 


NASA/ Johnson Space Center 
Houston, Texas 77058 


ATTN: 

J. 

D. Erickson/SF3 

( 1 ) 

ATTN: 

M. 

C. Trichel/SF3 

( 1 ) 

ATTN: 

L. 

F. Childs/SF 

( 1 ) 

ATTN: 

K. 

J. Demel/SF5 

( 1 ) 

ATTN: 

F. 

Weber/SF5 

( 1 ) 

ATTN: 

G. 

0 . Boatwright/SF3 

( 1 ) 

ATTN: 

K. 

Baker /SF4 

( 1 ) 

ATTN: 

H. 

G. DeVezin, Jr./FM 8 

( 1 ) 

ATTN: 

R. 

P. Heydorn/SF3 

( 1 ) 

ATTN: 

M. 

C. McEwen/SF3 

( 1 ) 

ATTN: 

D. 

H. Hay/ SF12 

( 1 ) 

ATTN: 

D. 

L. Amsbury/SF5 

( 1 ) 

ATTN: 

J. 

G. Garcia/SF3 

( 1 ) 

ATTN: 

F. 

G. Hall/ SF2 

( 1 ) 

ATTN: 

B. 

L. Carroll/C09 

( 1 ) 

ATTN: 

E. 

Laity/SF121 

( 2 ) 

ATTN: 

R. 

Shirkey/JM 6 

(4) 

ATTN: 

J. 

T. Wheeler/AT3 

( 1 ) 

ATTN: 

G. 

E. Graybeal/SF4 

( 2 ) 

ATTN: 

I. 

D. Brovne/SF3 

(5) 


IBM Corporation 
FSD Mail Code 56 
1322 State Park Drive 
Houston, Texas 77058 

ATTN: R. E. Oliver (l) 

Department of Mathematics 
Texas A&M University 
College Station, Texas 77843 

ATTN: L. F. Guseman, Jr. (l) 


ERIM 

P. 0. Box 86 l 8 

Ann Arbor, Michigan 48107 


ATTN: 

R. 

F. 

Nalepka 

( 1 ) 

ATTN: 

W. 

A. 

Malila 

( 1 ) 

ATTN: 

R. 

C. 

Cicone 

( 1 ) 


Kansas State University 
Department of Statistics, Calvin 19 
Statistical Lab 
Manhattan, Kansas 66506 

-ATTN: A. M. Feyerherm 


( 1 ) 



NAME 


NUMBER OF COPIES 


U. S. Department of Interior 
Geological Survey 
GSA Building, Room 5213 
Washington, D. C. 2021+2 

ATTN: Mr. W. A. Fischer 

NASA Wallops 

Wallops Station, Virginia 23337 
ATTN: Mr. James Bettle 

U. S. Department of Interior 
EROS Office 

Washington, D. C. 20242 
ATTN: Dr. Raymond W. Fary 

U. S. Department of Interior 
EROS Office 

Washington, D. C. 2021+2 
ATTN: Mr. William Hemphill 

NAV0CEAN0, Code 7001 
Bay St. Louis, MS 39520 

ATTN: Mr. J. W. Sherman, III 

University of Texas at Dallas 
Box 688 

Richardson, Texas 75080 
ATTN: Dr. Patrick L. Odell 

Department of Mathematics 
University of Houston 
Houston, Texas 77004 

ATTN: Dr. Henry Decell 

U. S. Department of Agriculture 
Statistical Reporting Service 
Room 4833, South Bldg. 

Washington, D. C. 20250 

ATTN: W. H. Wigton 

Goddard Space Flight Center 

National Aeronautics & Space Administration 

Greenbelt , Maryland 20771 

ATTN: Mr. W. Alford, 563 
ATTN: Dr. J. Barker, 923 


I 


( 1 ) 


( 1 ) 


( 1 ) 


( 1 ) 


( 1 ) 


( 1 ) 


( 1 ) 


( 1 ) 


( 1 ) 

( 1 ) 



NAME 


NUMBER OF COPIES 


EROS Data Center 

U. S. Department of Interior 

Sioux Falls, South Dakota 57198 

ATTN: Mr. G. Thorley 

U. S. Department of Agriculture 

Soil & Water Conservation Research Division 

P. 0. Box 267 

Weslaco, Texas 78596 

ATTN: Dr. Craig Wiegand 

U. S. Department of Interior 
USGA National Center 
Mail Stop 115 
Geography Program 
Reston, Virginia 22092 

ATTN: Dr. James R. Anderson 

Director, Remote Sensing Institute 
South Dakota State University 
Agriculture Engineering Building 
Brookings, South Dakota 57006 

ATTN: Mr. Victor I. Myers 

U. S. Department of Agriculture 

Forest Service 

2l+0 W. Prospect Street 

Fort Collins, Colorado 80521 

ATTN: Dr. Richard Driscoll 

University of California 
School of Forestry 
Berkeley, California 9^720 

ATTN: Dr. Robert Colwell 

Environmental Remote Sensing 
Applications Laboratory 
Oregon State University 
Corvallis, Oregon 97331 

ATTN: Dr. Barry J. Schrumpf 

U. S. Department of Interior 
Director, EROS Program 
Washington, D. C. 2021+2 

ATTN: Mr. J. M. Denoyer 
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( 1 ) 


( 1 ) 


( 1 ) 


( 1 ) 


( 1 ) 


( 1 ) 
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NUMBER OF 


Remote Sensing Laboratory- 
129 Mul ford Hall 
University of California 
Berkeley, California 94720 

ATTN: C. M. Hay 

Department of Mathematics 
University of Houston 
Houston, Texas 77004 

ATTN: H. P. Decell, Jr. 

NASA Lyndon B. Johnson Space Center 
Public Affairs Office, Code AP 
Houston, Texas 77058 

National Aeronautics and Space Administration 
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Department of Watershed Sciences 
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ATTN: Mr. J. P. Claybourne/AA-STA (l) 

Texas A&M University 
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